AD-A224  461 


OUC  FILE  COW 

NORSAR  Scientific  Report  No.  1-89/90 

Semiannual  Technical  Summary 

1  April  —  30  September  1989 


Kjeller,  December  1989 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 


DTIC 
OaEcTEn 

w  g  XJ 


UNCLASSIFIED  _ 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-0188 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 
VOT  APPT  TOART  17 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

Scientific  Rep.  1-89/90 


lb  RESTRICTIVE  MARKINGS 

NOT  APPLICABLE 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

APPROVED  FOR  PUBLIC  RELEASE 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

-  Scientific  Rep.  1-89/90 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 

_ _ „  .  „  «f  applicable)  ™ ~ 


NTNF/NORSAR 


6c  ADDRESS  {City,  State,  and  ZIP  Code) 

Post  Box  51 

N-2007  Kjeller,  Norway 


8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION  „  , 

Defence  Advanced 

_  Research.  Pro iects  Aeencv _ 


8c  ADDRESS  (City.  State,  and  ZIP  Code) 

1400  Wilson  Blvd. 


HQ/AFTAC/TTS 


7b  ADDRESS  {City,  State,  and  ZIP  Code) 

Patrick  AFB ,  FL  32925-6001 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

Contract  No.  F08606-89-C-0005 


10  SOURCE  OF  FUNDING  NUMBERS 


11  TITLE  (Include  Security  Classification) 

SEMIANNUAL  T~ ' MNICAL  SUMMARY,  1  APRIL  -  30 

12  PERSONAL  AUTHO 

L.B.  Loughran  (ed.  ) 

SCIENTIFIC  SUMMARY 


16  supplementary  notation 


December  1989 


NOT  APPLICABLE 


COSATI  CODES 


FIELD  GROUP  SUB-GROUP 


18  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

NORSAR,  NORWEGIAN  SEISMIC  ARRAY 


19  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 


(Please  see  reverse  side) 


20  DISTRIBUTION/AVAILABILITY  OF  ABSTRACT 
□  UNCLASSIFIED/UNLIMITED  0  SAME  AS  RPT 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

MR.  LEE  BRIDGES 


□  DTIC  USERS 


z:  abstract  security  class. fication 


iiKmpjKKiiaun 


22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

(407)  494-7765  AFTAC/TTS 


OD  Form  1473,  JUN  86 


Previous  editions  are  obsolete 


SECURITy  CLASSIFICATION  OF  THIS  PAGE 


Abstract 


This  Semiannual  Technical  Summary  describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Nor¬ 
wegian  Regional  Seismic  Array  (NORESS)  and  the  Arctic  Regional  Seismic 
Array  (ARCESS),for  the  period  1  April  -  30  September  1989. 

”  The  NORSAR  Detection  Processing  system  has  been  operated  through¬ 
out  the  period  with  an  average  uptime  of  94.9  as  compared  to  98.2  for  the 
previous  reporting  period.  A  total  of  1908  seismic  events  have  been  reported 
in  the  NORSAR  monthly  seismic  bulletin.  The  performance  of  the  continu¬ 
ous  alarm  system  and  the  automatic  bulletin  transfer  by  telex  to  AFTAC  has 
been  satisfactory.  Processing  of  requests  for  full  NORSAR/NORESS  data  on 
magnetic  tapes  has  progressed  according  to  established  schedules. 

The  satellite  link  for  transmitting  NORESS  data  in  real  time  to  the  U.S. 
has  had  an  average  uptime  of  97.8%  as  compared  to  99.9%  for  the  previous 
period.  On-line  detection  processing  and  data  recording  at  the  NORSAR 
Data  Center  (NDPC)  of  NORESS  and  ARCESS  data  have  been  conducted 
throughout  the  period,  with  an  average  uptime  of  92.1%  for  NORESS  and 
92.2%  for  ARCESS. 

Field  maintenance  activity  has  included  regular  preventive  maintenance 
at  all  array  sites  and  occasional  corrective  actions  when  required.  All  the 
NORSAR  communications  systems  have  been  checked  during  the  period,  and 
although  there  have  been  some  irregularities  in  connection  with  01A,  02C  and 
06C,  the  performance  has  been  satisfactory.  Modcomp  problems  began  at  the 
end  of  August  and  fault  finding  continued  throughout  the  rest  of  the  reporting 
period. 

Various  alternatives  for  a  possible  future  upgrade  of  the  NORSAR  system 
have  been  considered  during  the  period.  The  NORSAR  Detection  Processing 
has  been  transferred  to  a  more  powerful  IBM  computer  which  became  available 
after  moving  NORESS  to  a  SUN-system.  There  have  been  no  changes  in  the 
routine  processing  of  NORSAR  events. 

The  routine  detection  processing  of  NORESS  is  running  satisfactorily  on 
a  dedicated  SUN-3/280.  Event  processing,  which  can  be  performed  on  any 
of  the  SUN  workstations,  is  performed  in  two  steps:  (i)  Phase  estimation, 
and  (ii)  plot  and  epicenter  determination.  The  routine  detection  processing 
of  ARCESS  is  also  running  satisfactorily  on  a  separate  SUN-3/280.  Event 
processing  for  ARCESS  is  identical  to  that  of  NORESS. 
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The  “generalized  beamforming”  method,  which  has  been  developed  at 

NORSAR  and  represents  a  new  approach  to  associate  detected  seismic  phases 
from  a  regional  or  global  network,  has  been  further  tested.  A  data  base  of  77 
regional  and  local  seismic  events  recorded  by  NORESS  and  ARCESS  has  been 
used,  and  the  results  have  been  compared  to  output  of  the  Assess  procedure  in 
the  Intelligent  Monitoring  System  (IMS).  We  have  found  that  the  generalized 
beamforming  method  matches  or  exceed*  the  Assess  results  as  far  as  phase 
grouping  is  concerned  and  that  the  epicenter  determinations  are  also  quite 
accurate.  It  is  concluded  that  the  new  method  would  be  useful  to  supplement 
expert  systems  techniques.  The  recommended  approach  is  to  use  it  as  a  pre¬ 
processor,  and  subsequently  invoke  rule-based  and  script-based  algorithms  in 
order  to  further  refine  the  solutions. 


NORSAR  RMS  Lg  magnitude  measurements  as  published  in  previous 
Semiannual  reports^have  been  compared  tn-yieiTTclata  recently  published  by 
Soviet  scientists  for  Shagan  River  explosions.  The  correspondence  appears  to 
be  excellent,  and  considerably  better  than  the  correspondence  between  world¬ 
wide  mi,  and  these  yield  estimates.  However,  only  a  few  large  events  have  been 
available  for  this  study  and  more  data  will  be  needed  before  firm  conclusions 
may  be  drawn. „  Based  on  these  same  yield  data,  and  analyzing  the  smallest 
explosions  available,  it  appears  that  during  normal  noise  conditions  the  NOR¬ 
SAR  P-wave  detection  threshold  for  fully  coupled  Semipalatinsk  explosions  is 
below  1  kt,  evep  at  the  single  sensor  level.  Additional  significant  improvement 
would  be  obtaihablgjhjmugh  au^y  processing  of  these  data. 


The  stability  of  RMS  Lg  magnitudes  of  Semipalantinsk  explosions  as  a 
function  of  azimuth  and  distance  has  been  further  studied  by  analyzing  data 
from  the  China  Digital  Seismograph  Network  (CDSN).  Data  from  the  stations 
WMQ  and  HIA  have  been  processed  for  a  total  of  12  explosions  and  compared 
to  NORSAR  M(Lg).  The  resulting  scatter  in  the  magnitude  measurements 
is  very  small,  and  indicates  that  single  station  M(Lg)  may  be  measured  with 
an  accuracy  (one  standard  deviation)  of  0.03-0.04  magnitude  units.  This  cor¬ 
roborates  previous  results  comparing  data  from  NORSAR  and  IRIS  stations 
in  the  Soviet  Union.  Furthermore,  Lg  detectability  is  excellent  for  these  two 
CDSN  stations.  Thus  the  station  WMQ,  which  is  situated  at  960  km  from 
Semipalatinsk,  has  an  Lg  detection  threshold  estimated  at  2  magnitude  units 
lower  than  that  of  NORSAR. 

For 

A  new  data  acquisition  system  for  the  FINESA  array  in  Finland  has  been - 

developed  at  NORSAR,  and  installed  in  cooperation  with  the  University  of  tI 
Helsinki.  This  system  enables  real  time  acquisition  and  processing  of  data  from 
the  15  SPZ  seismometers  of  the  array  both  at  Helsinki  and  at  the  NORSAR  x 

data  center.  Initial  on-line  processing  has  been  established  both  in  Finland  _ 

and  Norway,  and  data  quality  appears  very  satisfactory.  The  FINESA  data 
are  currently  being  prepared  for  integration  into  the  Intelligent  Monitoring  , 
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System,  along  with  data  from  NORESS,  ARCESS  and  the  GERESS  system 
under  development.  Initially,  the  beam  deployment  for  FINESA  has  been 
adapted  from  the  NORESS/ ARCESS  deployments,  and  the  real  time  process¬ 
ing  methodology  is  essentially  identical. 

A  project  to  establish  a  mining  explosion  data  base  for  Fennoscandia  and 
adjacent  regions  has  been  initiated.  A  questionnaire  requesting  information 
on  mine  locations,  charge  sizes,  firing  practices  as  well  as  details  of  individual 
explosions  after  1984  has  initially  been  submitted  to  67  organizations  and 
companies  known  to  conduct  mining  activities.  So  far,  responses  have  been 
obtained  from  16  companies,  and  the  available  information  is  being  compiled 
into  a  data  base  conforming  to  the  2.8  CSS  format.  The  data  base  will  be 
continually  updated  and  expanded  as  more  information  is  obtained,  and  the 
aim  is  to  map  mining  activity  in  Fennoscandia  as  well  as  Central  Europe 
and  Western  USSR.  This  knowledge  will  ultimately  be  of  great  value  in  the 
regular  analysis  of  array  data,  and  will  also  be  the  basis  for  addressing  a  range 
of  interesting  and  important  topics  within  seismological  verification  research. 

A  study  has  been  conducted  modelling  the  transverse  components  of  Lg 
waves  and  relating  the  synthetic  seismograms  to  refraction  data  observed  at 
NORSAR.  It  is  shown  that  an  upper  crust  with  a  uniform  distribution  of  mi¬ 
crocracks  aligned  in  the  regional  stress  field  can  account  for  transverse  compo¬ 
nents  of  Lg  waves  induced  by  explosions.  The  preferred  model  to  explain  this 
in  terms  of  upper  crustal  anisotropy  is  a  distribution  of  liquid-filled  cracks  in 
the  upper  10  km  of  the  crust  with  a  crack  density  larger  than  0.02  and  most 
favorably  around  0.05.  This  model  is  consistent  with  previously  published 
observations  of  shear  wave  splitting  in  northern  Scandinavia. 

A  study  addressing  statistical  procedures  for  seismic  signal  detection  and 
estimation  using  small-array  data  has  been  undertaken.  A  generalization  of 
Capon’s  maximum  likelihood  technique  for  detection  and  estimation  of  seismic 
signals  is  introduced.  By  using  multidimensional  autoregressive  approximation 
of  array  seismic  noise.  Capon’s  group  filter  can  be  implemented  for  on-li no 
processing.  Such  autoregressive  adaptation  to  the  current  noise  matrix  power 
spectrum  yields  good  suppression  of  mutually  correlated  array  noise  pn  .esses. 
An  example  is  shown  of  this  technique  as  applied  to  detection  o',  a  small 
underground  explosion  (Semipalatinsk,  December  28,  1988)  recorded  at  the 
ARCESS  array. 
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1  Summary 


This  Semiannual  Technical  Summary  describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Nor¬ 
wegian  Regional  Seismic  Array  (NORESS)  and  the  Arctic  Regional  Seismic 
Array  (ARCESS)  for  the  period  1  April  -  30  September  1989. 

The  NORSAR  Detection  Processing  system  has  been  operated  through¬ 
out  the  period  with  an  average  uptime  of  94.9  as  compared  to  98.2  for  the 
previous  reporting  period.  A  total  of  1908  seismic  events  have  been  reported 
in  the  NORSAR  monthly  seismic  bulletin.  The  performance  of  the  continu¬ 
ous  alarm  system  and  the  automatic  bulletin  transfer  by  telex  to  AFTAC  has 
been  satisfactory.  Processing  of  requests  for  full  NORSAR/NORESS  data  on 
magnetic  tapes  has  progressed  according  to  established  schedules. 

The  satellite  link  for  transmitting  NORESS  data  in  real  time  to  the  U.S. 
has  had  an  average  uptime  of  97.8%  as  compared  to  99.9%  for  the  previous 
period.  On-line  detection  processing  and  data  recording  at  the  NORSAR 
Data  Center  (NDPC)  of  NORESS  and  ARCESS  data  have  been  conducted 
throughout  the  period,  with  an  average  uptime  of  92.1%  for  NORESS  and 
92.2%  for  ARCESS. 

Various  alternatives  for  upgrading  the  NORSAR  system  have  been  con¬ 
sidered  during  the  period.  The  Nanometrics  system  used  at  Yellowknife  and 
Hagfors  in  Sweden  has  features  which  seem  to  make  it  suitable  for  an  up¬ 
grade  of  the  NORSAR  array.  The  NORSAR  Detection  Processing  has  been 
transferred  to  a  more  powerful  IBM  computer  which  became  available  after 
moving  NORESS  to  a  SUN-system.  There  have  been  no  changes  in  the  routine 
processing  of  NORSAR  events. 

The  routine  detection  processing  of  NORESS  is  running  satisfactorily  on 
a  dedicated  SUN-3/280.  Event  processing,  which  can  be  performed  on  any  of 
the  Sun  workstations,  is  performed  in  two  steps:  Phase  estimation,  and  plot 
and  epicenter  determination.  The  routine  detection  processing  of  ARCESS 
is  also  running  satisfactorily  on  a  separate  SUN-3/280.  Event  processing  for 
ARCESS  is  identical  to  that  of  NORESS. 
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The  research  activity  is  summarized  in  Section  7.  Section  7.1  presents  a 
performance  test  of  the  generalized  beamforming  method  applied  to  a  data 
base  of  regional  events.  Section  7.2  discusses  NORSAR  P-wave  detection  and 
yield  estimation  of  selected  Semipalatinsk  explosions.  Section  7.3  presents 
an  analysis  of  data  from  the  Chinese  Digital  Seismograph  Network  for  Soviet 
nuclear  explosions.  Section  7.4  describes  a  new  data  acquisition  system  for 
FINESA.  Establishment  of  a  mining  explosion -data  base  is  presented  in  section 
7.5.  In  section  7.6,  transverse  components  of  explosion-induced  Lg  waves 
and  upper  crustal  anisotropy  are  discussed.  Section  7.7  discusses  statistically 
optimal  event  detection  using  small  array  data. 
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2  Operation  of  all  Systems 


2.1  Detection  Processor  (DP)  Operation 

There  have  been  93  breaks  in  the  otherwise  continuous  operation  of  the  NOR- 
SAR  online  system  within  the  current  6-month  reporting  interval.  The  uptime 
percentage  for  the  period  is  94.9  as  compared  to  98.2  for  the  previous  period. 

Fig.  2.1.1  and  the  accompanying  Table  2.1.1  both  show  the  daily  DP  down¬ 
time  for  the  days  between  1  April  and  30  September  1989.  The  monthly 
recording  times  and  percentages  are  given  in  Table  2.1.2. 


The  breaks  can  be  grouped  as  follows: 

a)  Hardware  failure  16 

b)  Stops  related  to  program  work  or  error  2 

c)  Hardware  maintenance  stops  11 

d)  Power  jumps  and  breaks  0 

e)  TOD  error  correction  20 

f)  Communication  lines  14 


The  total  downtime  for  the  period  was  225  hours  and  46  minutes.  The 
mean-time-between-failures  (MTBF)  was  1.9  days,  as  compared  to  1.8  for  the 
previous  period. 

J.  Torstveit 
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6 

1 

6 

2 

TOD  RETARED  23MS 

260 

0 

0 

12 

24 

MODCOMP  FAILURE 

195 

6 

0 

7 

45 

LINE  FAILURE 

261 

15 

10 

15 

51 

MODCOMP  FAILURE 

207 

15 

7 

17 

15 

CE  MAINTENANCE 

261 

17 

11 

17 

22 

MODCOMP  FAILURE 

209 

7 

31 

8 

5 

MACHINE  ERROR 

262 

5 

15 

5 

50 

MODCOMP  FAILURE 

212 

8 

41 

8 

42 

TOD  RETARED  30MS 

262 

7 

0 

7 

34 

MODCOMP  FAILURE 

212 

8 

49 

8 

54 

LINE  FAILURE 

262 

ii 

54 

12 

37 

CE  SERVICE  MODC. 

212 

8 

58 

9 

4 

LINE  FAILURE 

263 

6 

1 

6 

3 

TOD  RETARED  28MS 

212 

9 

10 

9 

12 

LINE  FAILURE 

263 

15 

53 

16 

4 

MODCOMP  FAILURE 

216 

15 

15 

24 

0 

MODCOMP  FAILURE 

264 

12 

47 

12 

52 

CE  SERVICE  MODC. 

217 

0 

0 

24 

0 

MODCOMP  FAILURE 

265 

19 

0 

24 

0 

MODCOMP  FAILURE 

218 

0 

0 

24 

0 

MODCOMP  FAILURE 

266 

0 

0 

8 

30 

MODCOMP  FAILURE 

219 

0 

0 

9 

0 

MODCOMP  FAILURE 

267 

0 

10 

8 

10 

MODCOMP  FAILURE 

219 

9 

48 

9 

52 

LINE  FAILURE 

267 

10 

36 

12 

0 

MODCOMP  FAILURE 

221 

13 

42 

14 

21 

MODCOMP  FAILURE 

268 

12 

38 

12 

55 

CE  SERVICE  MODC. 

222 

7 

50 

7 

53 

TOD  RETARED  22MS 

269 

2 

22 

4 

8 

MODCOMP  FAILURE 

228 

6 

9 

6 

11 

TOD  RETARED  13MS 

269 

12 

13 

12 

49 

MODCOMP  FAILURE 

229 

7 

44 

8 

15 

CE  MAINTENANCE 

269 

13 

13 

13 

43 

MODCOMP  FAILURE 

230 

8 

55 

9 

3 

CPU  FAILURE 

269 

13 

57 

14 

18 

MODCOMP  FAILURE 

234 

6 

6 

6 

8 

TOD  RETARED  23MS 

269 

15 

5 

16 

40 

CE  SERVICE  MODC. 

237 

3 

41 

7 

4 

MODCOMP  FAILURE 

269 

20 

8 

24 

0 

MODCOMP  FAILURE 

239 

5 

38 

19 

33 

MODCOMP  FAILURE 

270 

0 

0 

6 

17 

MODCOMP  FAILURE 

239 

23 

36 

24 

0 

MODCOMP  FAILURE 

270 

8 

48 

9 

7 

CE  SERVICE  MODC. 

240 

0 

0 

6 

34 

MODCOMP  FAILURE 

270 

20 

7 

20 

59 

MODCOMP  FAILURE 

240 

9 

20 

10 

55 

MODCOMP  FAILURE 

27  1 

8 

34 

19 

30 

CE  SERVICE  MODC. 

240 

19 

2 

24 

0 

MODCOMP  FAILURE 

272 

9 

18 

9 

23 

CE  SERVICE  CPU 

241 

0 

0 

5 

22 

MODCOMP  FAILURE 

273 

23 

16 

23 

59 

MODCOMP  FAILURE 

241 

6 

31 

9 

40 

MODCOMP  FAILURE 

Table  2.1.1  Daily  DP  downtime  in  the  period  1  April  -  30  September  1989. 


r> 


Month 

DP  uptime 
hours 

DP  uptime 
*/. 

No .  of 

DP  breaks 

No.  of  days 
with  breaks 

DP  MTBF* 
(days) 

APR 

715.60 

99.4 

6 

5 

4.3 

MAY 

743.88 

99.9 

7 

6 

3.9 

JUN 

719.75 

99.9 

10 

7 

2.7 

JUL 

741.27 

99.4 

10 

6 

2.8 

AUG 

617.67 

83.0 

21 

16 

1.2 

SEP 

630.10 

87.5 

39 

26 

0.7 

94.9 

93 

66 

1.9 

♦Mean-time-between-failures  =  total  uptime/no.  of  up  intervals. 


Table  2.1.2  Online  system  performance.  1  April  -  30  September  1989. 
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2.2  Array  communications 


General 

Table  2.2.1  reflects  the  performance  of  the  communications  system  through¬ 
out  the  reporting  period. 

The  comm,  adapter  2701  SDAII  failed  in  April,  causing  "garbled-’  data 
between  the  IBM  4341  and  the  Modcomp.  Already  the  next  day  (24  April)  the 
fault  was  found  and  the  systems  were  started  again.  Although  irregularities  in 
connection  with  01A.  02C  and  06C  occurred,  the  performance  of  the  comm, 
systems  may  be  characterized  as  most  satisfactory. 

In  May  02B  power  problems  started  in  connection  with  lightning  in  the 
area.  NTA  carrier  systems  affected  01 B,  03C.  04C  and  06C. 

June  was  characterized  by  many  outages.  Specially  02B  and  06C  were  af¬ 
fected.  02B  lost  power  after  lightning,  and  06C  had  SLEM  problems.  Besides. 
NTA  transmission  system  (in  Oslo)  affected  all  subarrays. 

Also  in  July  most  subarrays  were  affected  ( -01 A .  01B),  02B  by  lightning  (as 
usual).  02C  by  a  broken  cable.  03C,  04C  and  06C  were  affected  by  NTA  cable 
work  (planned). 

Modcomp  problems  started  in  August  (25),  and  02B  had  been  down  since 
28  July  due  to  power  outages.  Estimated  comm,  error  figures  were  calculated 
while  the  Modcomp  was  down. 

In  September  Modcomp  fault  findings  continued.  01A  was  affected  by 
a  broken  comm,  cable  (week  38).  02C  was  affected  week  39.  and  06C  by 
excessive  error  rate  (week  36).  Estimated  average  error  figures  for  each  comm, 
system  were  calculated  during  Modcomp  stops. 

Detailed  summary 

April  (weeks  1 J- 17),  3-30-4-89 

Week  16  01A  and  02C  were  affected  by  comm,  line  irregularities.  06C 
was  down  week  17,  reason  not  stated.  23  April  the  comm,  adapter  IBM  2701 
SDAII  failed  causing  "garbled”  data  transfer  between  the  IBM  4341  and  the 
Modcomp.  24  April  a  bad  card  was  found  and  the  system  started  again. 


May  (weeks  18-22),  1.5-4-6.89 


02B  power  affected  after  lightning  (week  19).  Local  power  company  re¬ 
paired  the  damage  19  May.  22  May  the  subarray  was  visited  by  NMC  staff 
and  the  SLEM  was  reset.  01B,  04C  and  06C  were  also  affected  week  20,  01B  by 
3.6%,  04C  and  06C  approx.  10.7%.  NTA  carrier  equipment  was  the  probable 
reason. 

An  increasing  number  of  resynchronizations  were  observed  over  several 
days.  NTA/Lillestrom  confirmed  carrier  equipment  instability  and  the  prob¬ 
lem  was  solved  (week  21).  25  May  the  same  problem  was  observed,  and  all 
subarrays  were  affected.  NTA  found  a  60-channel  group  in  Oslo  was  faulty. 

I  June  a  coaxial  cable  carrying  03C  and  04C  data  was  damaged.  06C  and 
the  satellite  control  line  between  Eik  earth  station  and  NORESS  were  also 
affected.  NTA/Hamar  switched  03C,04C  comm,  paths  to  a  spare  coax  carrier 
system. 

June  (weeks  23-26),  5.6-2.7.S9 

02B  was  out  of  operation  from  5  June,  and  resumed  operation  14  June  after 
SLEM  reset.  Faulty  transmission  system  in  Oslo  affected  all  comm,  systems 

II  June.  A  faulty  transmission  system  between  Oslo  and  Hamar  affected  02B. 
03C,  04C  and  06C  23  June.  02B  was  down  approx.  20  hours  between  21  and 
22  June,  and  went  down  again  30  June.  It  remained  down  throughout  the 
week  (26).  22  June  we  lost  contact  with  the  subarrav.  29  June  the  subarray 
resumed  operation. 

July  (weeks  27-30),  3-30.7.89 

02B,  which  went  down  30  June  due  to  lightning  in  the  area,  was  turned  on 
5  July  after  SLEM  reset.  28  July  the  subarray  was  affected  again  by  a  power 
outage,  caused  by  lightning,  and  remained  down  throughout  the  month.  02C 
was  affected  4  July  due  to  a  broken  cable.  03C.04C  and  06C  were  affected  5 
and  6  July  in  connection  with  planned  NTA/Hamar  cable  work. 

August  (weeks  31-35),  31.7-3.9.S9 

02B  was  turned  on  2  August  after  having  been  down  since  28  July.  The 
Modcomp  reliability  was  reduced  since  25  August,  with  many  stops.  Per  31 
August  TM-data,  responsible  for  corrective/preventive  maintenance  of  the  ma¬ 
chine.  had  not  been  able  to  find  the  items  causing  the  stops.  The  performance 
of  the  comm,  systems  (excl.  the  Modcomp)  was  in  the  period  most  satisfac¬ 
tory,  and  estimated  average  error  figures  for  the  systems  while  Modcomp  was 
stopped  were  calculated. 


September  (weeks  36-39),  ).  9-1 .10.89 

Modcomp  fault  finding  continued  in  September  in  spite  of  TM-data's  ef¬ 
forts  to  solve  the  problem.  The  next  step  will  be  to  analyze  the  CPU  mi¬ 
crocodes.  As  far  as  comm,  systems  are  concerned,  01A  was  affected  by  a 
broken  comm,  cable  week  38.  02C  was  down  week  39  of  unknown  reason.  Ex¬ 
cessive  error  rate  could  be  observed  week  36  in  connection  with  06C.  During 
the  Modcomp  stop  periods,  estimated  average  error  figures  were  calculated  for 
each  subarray  comm,  system. 


O.A.  Hanson 


Sub- 

Apr  (4) 

May  (5) 

Jun  (4) 

Jul  (4) 

Aug  (5) 

Sep  (4) 

Average 

arrays 

3-30.4 

1. 5-4.6 

5. 6-2.7 

3-30.7 

31.7-3.9 

4.9-1.10 

1/2  Year 

01 A 

*0.5100 

0.030 

*1.370 

0.020 

0.007 

‘°)*0.050 

0.330 

01B 

0.0002 

*0.770 

*1.390 

*0.060 

0.002 

0.001 

0.370 

02B 

0.0040 

‘1*0.004 

2>*N/A 

4>*0.020 

9>*0.003 

0.004 

13  '0 .007 

02C 

*0.4400 

0.030 

*0.930 

5>*0.020 

0.003 

“1*0.002 

0.240 

03C 

0.0040 

0.390 

*6.730 

6'*0.020 

0.007 

0.003 

1.190 

04C 

0.0010 

*2.540 

*2.710 

7)*0  030 

0.004 

0.002 

0.b30 

06C 

0.1600 

*2.950 

3)N/A 

8)*1 .230 

0.002 

1 2 1*0 .060 

,4)0.880 

AVER 

0.1600 

0.960 

2.630 

0.200 

0.004 

0.020 

0.560 

[  Less 

02B.06C 

*See  Section  2.2  regarding  figures  preceded  by  an  asterisk. 

Figures  representing  error  rate  (in  per  cent)  preceded  by  a  number  1),  2),  etc.,  are  related  to 
legend  below. 


11  Power  outages  weeks  19-22 

2)  Power  outages  weeks  23-26 

3)  Outages  weeks  23-26  (faulty  SLEM) 

4>  Average  2  weeks  (28,29),  27,30  N/A 

5-8l  Average  3  weeks  (28,29,30),  27  N/A 
91  Average  4  weeks  (32-35),  31  N/A 


*°1  Average  3  weeks  (36-38),  39  N/A 
111  Average  3  weeks  (36-38),  39  N/A 
121  Average  3  weeks  (36-38),  39  N/A 

13)  Average  5  months,  weeks  23-26  N/A 

14)  Average  5  months,  weeks  23-26  N/A 


Table  2.2.1  Communications  performance.  The  numbers  represent  error  rates 
in  per  cent  based  on  total  transmitted  frames/week  (1  April  -  30  September 
1989). 
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2.3  Event  Detection  operation 

In  Table  2.3.1  some  monthly  statistics  of  the  Detection  and  Event  Processor 
operation  are  given.  The  table  lists  the  total  number  of  detections  (DPX) 
triggered  by  the  on-line  detector,  the  total  number  of  detections  processed  by 
the  automatic  event  processor  (EPX)  and  the  total  number  of  events  accepted 
after  analyst  review  (Teleseismic  phases,  core  phases  and  total). 


Total 

DPX 

Total 

EPX 

Accepted  events 
P-phases  Core  phases 

Sum 

Daily 

Apr  89 

8900 

1203 

302 

62 

364 

12.1 

May  89 

8250 

1217 

291 

92 

383 

12.4 

Jun  89 

7325 

1286 

222 

55 

277 

9.2 

Jul  89 

7350 

1070 

244 

53 

297 

9.6 

Aug  89 

5800 

980 

254 

60 

314 

10.1 

Sep  89 

9450 

1412 

221 

52 

273 

9.1 

1534 

374 

1908 

10.4 

Table  2.3.1.  Detection  and  Event  Processor  statistics,  1  April  -  30  September 
1989. 


B.Kr.  Hokland 
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3  Operation  of  NORESS  and  ARCESS 

3.1  Satellite  transmission  of  NORESS  data  to  the  U.S. 


The  satellite  transmission  of  data  to  the  U.S.  from  the  NORESS  field  instal¬ 
lation  was  interrupted  on  16  occasions,  as  shown  in  Table  3.3.1.  Only  one 
outage  period  was  due  to  failure  of  the  ground  station  equipment. 


25  May 

1649 

to 

1713 

Poser  failure 

due  to 

thunderstorm 

1 

Jim 

1055 

to 

1145 

Control  lines 

down 

7 

Jun 

0202 

to 

0416 

Power  failure 

9 

Jim 

0844 

to 

0945 

Power  failure 

13 

Jun 

2235 

to 

14 

Jun 

0152 

Power  failure 

21 

Jun 

0930 

to 

1225 

Power  failure 

due  to 

thunderstorm 

21 

Jun 

2227 

to 

2229 

Control  lines 

down 

22 

Jun 

2332 

to 

2334 

Control  lines 

down 

23 

Jun 

2046 

to 

24 

Jun 

0810 

Satellite  modem  down 

29 

Jun 

0926 

to 

1059 

Power  failure 

due  to 

thunderstorm 

29 

Jun 

1616 

to 

1834 

Power  failure 

due  to 

thunderstorm 

29 

Jun 

1900 

to 

1941 

Power  failure 

due  to 

thunderstorm 

23 

Jill 

1611 

to 

26 

Jul 

1118 

Power  failure 

due  to 

thunderstorm 

26 

Jul 

1229 

to 

1339 

Repair  work  due  to  thunderstorm 

31 

Jul 

2004 

to 

2151 

Control  lines 

down 

1 

Aug 

0136 

to 

0139 

Control  lines 

down 

Table  3.1.1  Outage  periods  for  NORESS  satellite  transmission  system  April 
-  September  1989. 


The  total  uptime  for  the  NORESS  Earth  Station  for  satellite  transmission 
of  data  to  the  U.S.  was  97.8%  as  compared  to  99.9%  for  the  previous  period. 

The  ground  station  for  continuous  transmission  of  NORESS  data  to  the 
U.S.  was  closed  down  on  1  October  1989.  From  this  date,  NORESS  and 
ARCESS  data  can  be  obtained  in  the  U.S.  at  the  Center  for  Seismic  Studies 
(CSS),  using  the  new  dedicated  satellite  link  established  between  NDPC  and 
CSS. 
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3.2  Recording  of  NORESS  data  at  NDPC,  Kjeller 


Table  3.2.1  lists  the  main  outage  times  and  reasons.  As  can  be  seen,  the  main 
reasons  for  the  outage  are  power  failures  at  the  HUB  and  line  failures. 

The  average  recording  time  was  92.1%  as  compared  to  98.5%  for  the  pre¬ 
vious  period. 


Date 

Time 

Cause 

25  May 

1654-1817 

Power  failure  HUB 

due  to  thunderstorm 

7  Jun 

0208-0549 

Power  failure  HUB 

due  to  thunderstorm 

9  Jun 

0553-0601 

Power  failure  HUB 

due  to  thunderstorm 

9  Jun 

0852-0915 

Power  failure  HUB 

due  to  thunderstorm 

13  Jun 

2243- 

Power  failure  HUB 

due  to  thunderstorm 

14  Jun 

-0608 

Power  failure  HUB 

due  to  thunderstorm 

20  Jun 

0726-0736 

Transmission  line 

failure 

21  Jun 

0931-1447 

Power  failure  HUB 

due  to  thunderstorm 

23  Jun 

0832-0909 

Transmission  line 

failure 

29  Jun 

0531-0542 

Power  failure  HUB 

29  Jun 

0935-1148 

Power  failure  HUB 

29  Jun 

1440-1606 

Power  failure  HUB 

30  Jun 

0706-0720 

Transmission  line 

failure 

3  Jul 

0501-1238 

Power  failure  HUB 

4  Jul 

0855-1225 

Transmission  line 

failure 

7  Jul 

1334- 

Transmission  line 

failure 

8  Jul 

-1223 

Transmission  line 

failure 

13  Jul 

1105-1234 

Transmission  line 

failure 

16  Jul 

0949-1336 

Power  failure  HUB 

19  Jul 

0611-0622 

Transmission  line 

failure 

23  Jul 

1655- 

Power  failure  HUB 

due  to  thunderstorm 

3  Aug 

-1628 

Power  failure  HUB 

due  to  thunderstorm 

18  Aug 

1231-1255 

Hardware  failute  at  NDPC 

3  Sep 

0626-0831 

Power  failure  HUB 

20  Sep 

1255-1333 

Hardware  maintenance  at  NDPC 

20  Sep 

1621-1732 

Transmission  line 

failure 

Table  3.2.1  Interruptions  in  NORESS  recordings  at  NDPC.  April-September 
1989. 

Monthly  uptimes  for  the  Noress  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera- 
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tion)  affecting  this  task  were  as  follows: 


April 

99.6% 

May 

99.8% 

June 

96.8% 

July 

§8.2% 

August 

88.5% 

September 

99.4% 

Fig.  3.2.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  NORESS  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 
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Fig.  3.2.1  NORESS  data  recording  uptime  for  April  (top).  May  (middle)  and 
June  (bottom)  19S9. 
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Fig.  3.2.1  NORESS  data  recording  uptime  for  July  (top),  August,  (middle) 
and  September  (bottom)  1989. 
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3.3  Recording  of  ARCESS  data  at  NDPC,  Kjeller 


The  main  reasons  causing  most  of  the  ARCESS  outage  in  the  period  are: 
Hardware  failure  at  NDPC,  transmission  failure,  and  maintenance  work.  A 
bad  cable  w'as  found  causing  most  of  the  hardware  failures  at  NDPC  in  the 
period. 

The  average  recording  time  was  92.2%  as  compared  to  98.5%  for  the  pre¬ 
vious  period. 


Date 

Time 

Duration  Cause 

18  Apr 

0209-0508 

Hardware  failure  at  NDPC 

18  Apr 

2100- 

Hardware  failure  at  NDPC 

19  Apr 

-0517 

Hardware  failure  at  NDPC 

1  May 

1040-1122 

Hardware  failure  at  NDPC 

3  May 

0600-0725 

Hardware  failure  at  NDPC 

18  May 

0720-1603 

Hardware  work  at  NDPC 

18  May 

1618-1658 

Maintenance  at  NDPC 

30  May 

0904-1334 

Change  of  transmission  satelite 

2  Jun 

1200-1252 

Maintenance  at  HUB 

9  Jun 

0646-0707 

Maintenance  at  HUB 

13  Jun 

0016-0517 

Hardware  failure  at  NDPC 

17  Jun 

2245- 

Transmission  line  failure 

18  Jun 

-0057 

Transmission  line  failure 

19  Jun 

0358-0505 

Hardware  failure  at  NDPC 

6  Jul 

1418- 

Hardware  failure  at  NDPC 

7  Jul 

-0608 

Hardware  failure  at  NDPC 

11  Jul 

0537-0627 

Hardware  failure  at  NDPC 

11  Jul 

0719-0739 

Hardware  failure  at  NDPC 

12  Jul 

2205- 

Transmission  line  failure 

13  J_1 

-0434 

Transmission  line  failure 

13  Jul 

0627-0643 

Hardware  failure  at  NDPC 

15  Jul 

0540- 

Satellite  modem  failure 

21  Jul 

-0141 

Satellite  modem  failure 

28  Jul 

0741-0902 

Hardware  failure  at  NDPC 

4  Aug 

1301-1358 

Hardware  failure  at  NDPC 

10  Aug 

0317-0529 

Hardware  failure  at  NDPC 

16  Aug 

0844-0904 

Hardware  failure  at  NDPC 

16  Aug 

0929-0944 

Hardware  failure  at  NDPC 

16  Aug 

2015-2031 

Hardware  failure  at  NDPC 

16  Aug 

2054- 

Hardware  failure  at  NDPC 

17  Aug 

-0504 

Hardwarw  failure  at  NDPC 
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17  Aug 

0717-0736 

Hardware  failure  at  NDPC 

24  Aug 

0040-01520 

Hardware  failure  at  NDPC 

24  Aug 

1627- 

Hardware  failure  at  NDPC 

25  Aug 

-0516 

Hardware  failure  at  NDPC 

26  Aug 

0936- 

Transmission  line  failure 

27  Aug 

-1909 

Transmission  line  failure 

4  Sep 

1250-1332 

Maintenance  at  HUB 

8  Sep 

1233-1704 

Maintenance  at  NDPC 

12  Sep 

0653- 

Maintenance  at  antenna  at 

NDPC 

13  Sep 

-0638 

Maintenance  at  antenna  at 

NDPC 

13  Sep 

0837-1009 

Maintenance  at  antenna  at 

NDPC 

13  Sep 

1355-1433 

Maintenance  at  antenna  at 

NDPC 

14  Sep 

0720-1540 

Maintenance  at  antenna  at 

NDPC 

15  Sep 

1023-1106 

Adjusting  the  antenna  at  NDPC 

17  Sep 

0324-1136 

Hardware  failure  at  NDPC 

24  Sep 

1036-1200 

Hardware  failure  at  NDPC 

Table  3.3.1  The  main  interruptions  in  ARCESS  recordings  at  NDPC.  April- 
September  1989. 


Monthly  uptimes  for  the  ARCESS  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera¬ 
tion)  affecting  this  task  were  as  follows: 


April 

98.4% 

May 

97.8% 

June 

98.7% 

July 

77.8% 

August 

89.3% 

September 

91.3% 

Fig.  3.3.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  ARCESS  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 

J.  Torstveit 
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Fig.  3.3.1  ARCESS  data  recording  uptime  for  April  (top).  May  (middle)  and 
June  (bottom)  19S9. 
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Fig.  3.3.1  ARCESS  data  recording  uptime  for  July  (top),  August  (middle) 
and  September  (bottom)  1989. 
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3.4  Event  detection  operation 
NORESS  detections 

The  number  of  detections  (phases)  reported  during  day  091,  1989  through 
day  273,1989  was  38803,  giving  an  average  of  246  detections  per  processed  day. 
During  this  period.  25  days  were  not  processed  due  to  NORESS  downtime  (see 
paragraph  2.2). 

Table  3.4.1  shows  daily  and  hourly  distribution  of  detections  for  NORESS. 

Events  automatically  located  by  NORESS 

During  days  091  through  273,  2166  local  and  regional  events  were  located 
by  NORESS,  based  on  automatic  association  of  P-  and  S-tvpe  arrivals.  This 
gives  an  average  of  13.7  events  per  processed  day  ( 158  days  processed).  66  % 
of  these  events  are  within  300  km,  and  87  %  of  these  events  are  within  1000 
km. 

ARCESS  detections 

The  number  of  detections  (phases)  reported  during  day  091,  1989  through 
day  273,1989  was  56375.  giving  an  average  of  348  detections  per  processed 
day.  During  this  period.  21  days  were  not  processed  (see  paragraph  2.3). 

Table  3.4.2  shows  daily  and  hourly  distribution  of  detections  for  ARCESS. 

Events  automatically  located  by  ARCESS 

During  days  091  through  273,  2732  local  and  regional  events  were  located 
by  ARCESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This 
gives  an  average  of  16.9  events  per  processed  day  (162  days  processed).  45  % 
of  these  events  are  within  300  km,  and  78  %  of  these  events  are  within  1000 
km. 

J.  Fyen 
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Table  3.4.1  Daily  and  hourly  distribution  of  NORESS  detections.  For  each 
day  is  shown  the  number  of  detections  within  each  hour  of  the  day,  and  the 
number  of  detections  for  that  day.  The  statistics  at  the  end  of  the  table  give 
total  number  of  detections  for  each  hour  and  the  total  sum  of  detections  during 
the  period.  The  averages  show  number  of  processed  days,  hourly  distribution 
and  average  per  processed  day. 
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Tuesday 

249 

3 

9 

14 

12 

2 

3 

5 

6 

13 

13 

19 

31 

10 

12 

20 

8 

11 

20 

21 

2 

15 

2 

2 

8 

261 

Sep 

06 

Wednesday 

250 

2 

10 

9 

2 

2 

5 

7 

6 

16 

1 

24 

7 

31 

24 

21 

5 

20 

14 

2 

1 

13 

1 

6 

13 

242 

Sep 

07 

Thursday 

251 

3 

9 

16 

11 

2 

3 

13 

13 

26 

19 

13 

21 

17 

18 

14 

15 

6 

8 

15 

23 

10 

8 

3 

15 

301 

Sep 

08 

Friday 

252 

2 

6 

25 

1 

6 

3 

4 

4 

5 

4 

7 

4 

7 

20 

5 

2 

9 

8 

7 

6 

7 

2 

3 

3 

150 

Sep 

09 

Saturday 

253 

2 

1 

3 

9 

3 

7 

2 

5 

5 

5 

4 

2 

5 

8 

9 

18 

10 

14 

2 

19 

5 

13 

7 

10 

168 

Sep 

10 

Sunday 

254 

16 

4 

21 

1 

6 

6 

13 

2 

15 

10 

3 

21 

10 

30 

17 

15 

19 

39 

5 

1 

11 

5 

15 

20 

305 

Sep 

11 

Monday 

255 

15 

7 

16 

5 

3 

6 

6 

15 

29 

8 

9 

25 

9 

8 

22 

5 

18 

16 

5 

5 

14 

0 

11 

2 

259 

Sep 

12 

Tuesday 

256 

12 

2 

16 

7 

10 

8 

7 

12 

20 

15 

11 

15 

12 

35 

23 

15 

12 

9 

2 

4 

12 

8 

11 

4 

282 

Sep 

13 

Wednesday 

257 

6 

4 

11 

2 

4 

13 

4 

3 

18 

10 

9 

15 

25 

16 

8 

6 

11 

7 

13 

17 

3 

1 

6 

7 

219 

Sep 

14 

Thursday 

258 

4 

6 

16 

3 

1 

5 

13 

10 

11 

19 

10 

13 

9 

7 

6 

12 

16 

6 

13 

4 

17 

1 

6 

6 

214 

Sep 

15 

Friday 

259 

3 

4 

34 

15 

9 

13 

8 

5 

2 

4 

10 

16 

4 

22 

42 

5 

5 

11 

6 

3 

1 

0 

3 

13 

238 

Sep 

16 

Saturday 

260 

7 

36 

13 

11 

10 

10 

13 

2 

7 

11 

0 

2 

9 

3 

7 

17 

12 

14 

5 

4 

2 

4 

6 

5 

210 

Sep 

17 

Sunday 

261 

11 

0 

14 

4 

4 

8 

11 

3 

4 

5 

8 

7 

18 

14 

11 

12 

7 

12 

2 

2 

13 

8 

6 

15 

199 

Sep 

18 

Monday 

262 

5 

9 

20 

6 

15 

0 

6 

5 

16 

11 

8 

10 

16 

22 

9 

11 

12 

18 

6 

1 

16 

5 

4 

12 

243 

Sep 

19 

Tuesday 

263 

13 

3 

11 

3 

9 

13 

0 

9 

21 

5 

11 

12 

12 

6 

15 

15 

2 

5 

9 

0 

1 1 

0 

6 

7 

198 

Sep 

20 

Wednesday 

264 

7 

3 

7 

0 

1 

2 

4 

6 

6 

13 

6 

7 

13 

20 

5 

11 

12 

14 

9 

6 

8 

1 

8 

11 

180 

Sep 

21 

Thursday 

265 

4 

4 

15 

5 

10 

1 

7 

7 

13 

19 

11 

17 

11 

10 

6 

3 

5 

10 

11 

12 

18 

4 

0 

10 

213 

Sep 

22 

Friday 

266 

0 

2 

17 

1 

3 

12 

1 

2 

5 

3 

2 

6 

2 

7 

6 

5 

1 

5 

13 

6 

2 

1 

3 

1 

106 

Sep 

23 

Saturday 

Table  3.4.1  (cont.) 


24 


267 

5 

2 

2 

5 

6 

1 

0 

8 

21 

3 

6 

9 

6 

13 

13 

4 

20 

4 

10 

0 

0 

15 

6 

3 

162 

Sep 

24 

Sunday 

268 

12 

5 

11 

11 

2 

4 

3 

17 

0 

7 

21 

10 

8 

9 

12 

13 

6 

27 

16 

3 

7 

3 

2 

3 

212 

Sep 

25 

Monday 

269 

6 

6 

8 

12 

3 

11 

2 

11 

4 

8 

7 

5 

9 

18 

19 

19 

15 

13 

9 

17 

2 

24 

0 

11 

239 

Sep 

26 

Tuesday 

270 

9 

9 

2 

13 

4 

5 

5 

10 

1 

10 

9 

11 

18 

23 

22 

5 

4 

20 

12 

5 

5 

9 

4 

8 

223 

Sep 

27 

Wednesday 

271 

4 

6 

8 

19 

4 

5 

6 

5 

7 

12 

8 

7 

14 

27 

19 

13 

3 

9 

13 

4 

1 

9 

10 

6 

219 

Sep 

28 

Thursday 

272 

19 

11 

2 

10 

6 

9 

0 

3 

3 

10 

7 

14 

17 

20 

12 

2 

4 

15 

4 

8 

12 

7 

1 

3 

199 

Sep 

29 

Friday 

273 

10 

7 

1 

15 

4 

5 

6 

8 

6 

1 

4 

2 

4 

3 

12 

5 

0 

5 

7 

3 

3 

4 

0 

6 

121 

Sep 

30 

Saturday 

NRS 

00 

01 

02 

03 

04 

05 

06 

07 

08 

09 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

Sum 

1552 

1514 

1158 

1274 

1443 

1978 

1877 

1737 

2026 

1368 

1220 

2281 

1655 

2296 

1183 

1225 

1514 

1640 

2037 

1696 

1582 

1583 

1677 

1287 

38803 

Total 

sum 

158 

10 

10 

15 

10 

7 

7 

8 

8 

10 

9 

10 

13 

13 

12 

11 

11 

10 

13 

10 

9 

11 

8 

8 

14 

246 

Total 

average 

126 

9 

8 

13 

8 

5 

5 

6 

7 

9 

8 

10 

12 

14 

12 

10 

10 

9 

12 

10 

8 

10 

6 

7 

14 

221 

Average  workdays 

57 

10 

10 

12 

10 

9 

9 

7 

6 

7 

7 

8 

9 

6 

7 

7 

8 

7 

9 

6 

7 

7 

8 

7 

10 

192 

Average  weekends 

Table  3.4.1  (cont.) 
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FRS 

Day 

91 

. FKX  Hourly  distribution 

00  01  02  03  04  05  06  07 

17  6  8  5  27  6  1  4 

of 

08 

6 

detections 

09  10  11  12 

6  15  10  12 

13 

14 

14 

3 

15 

0 

16 

4 

17 

7 

18 

9 

19 

8 

20 

1 

21 

16 

22 

23 

23 

7 

Sum 

215 

Date 

Apr  01 

Saturday 

92 

3 

11 

6 

3 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

24 

Apr 

02 

Sunday 

93 

0 

0 

0 

0 

0 

2 

4 

3 

6 

2 

2 

2 

6 

7 

12 

12 

1 

10 

8 

12 

8 

7 

12 

14 

130 

Apr 

03 

Monday 

94 

4 

5 

2 

6 

5‘ 

10 

6 

12 

16 

25 

8 

24 

15 

8 

12 

11 

8 

6 

10 

5 

6 

6 

9 

0 

219 

Apr 

04 

Tuesday 

95 

3 

5 

2 

5 

1 

4 

7 

12 

2 

4 

13 

19 

18 

11 

12 

12 

6 

21 

9 

5 

1 

7 

19 

7 

205 

Apr 

05 

Wednesday 

96 

10 

6 

8 

7 

2 

6 

12 

1 

14 

17 

13 

16 

2 

11 

10 

9 

14 

20 

3 

15 

11 

3 

13 

5 

228 

Apr 

06 

Thursday 

97 

13 

0 

7 

5 

2 

10 

2 

18 

14 

5 

8 

17 

20 

16 

5 

7 

7 

14 

5 

6 

10 

13 

21 

3 

228 

Apr 

07 

Friday 

98 

13 

10 

1 

13 

7 

1 

2 

11 

15 

17 

11 

21 

2 

2 

5 

9 

7 

1 

4 

3 

5 

2 

9 

2 

173 

Apr 

08 

Saturday 

99 

5 

1 

6 

4 

10 

11 

1 

4 

1 

10 

7 

13 

17 

3 

4 

6 

5 

6 

10 

6 

10 

2 

22 

3 

167 

Apr 

09 

Sunday 

100 

3 

3 

3 

3 

1 

4 

3 

2 

7 

8 

9 

6 

16 

22 

6 

4 

3 

3 

13 

12 

6 

5 

10 

6 

158 

Apr 

10 

Monday 

101 

4 

7 

4 

5 

12 

9 

2 

23 

7 

35 

9 

9 

9 

13 

12 

15 

1 

1 

-8 

7 

4 

6 

17 

0 

219 

Apr 

11 

Tuesday 

102 

3 

6 

9 

16 

27 

8 

9 

8 

4 

4 

15 

32 

21 

3 

11 

4 

3 

3 

9 

8 

7 

2 

19 

3 

234 

Apr 

12 

Wednesday 

103 

1 

10 

4 

7 

16 

15 

6 

20 

7 

4 

11 

13 

7 

6 

12 

7 

11 

14 

6 

6 

2 

6 

14 

8 

213 

Apr 

13 

Thur  sday 

104 

6 

6 

4 

1 

1 

6 

6 

14 

16 

19 

15 

24 

25 

20 

12 

18 

19 

9 

8 

13 

3 

5 

20 

15 

285 

Apr 

14 

Friday 

105 

4 

9 

5 

2 

2 

2 

12 

3 

22 

9 

23 

29 

2 

7 

8 

6 

1 

4 

8 

j 

10 

24 

7 

5 

207 

Apr 

15 

Saturday 

106 

7 

6 

1 

1 

7 

11 

20 

10 

22 

4 

8 

9 

7 

14 

1 

5 

10 

2 

8 

7 

20 

15 

29 

17 

241 

Apr 

16 

Sunday 

107 

15 

17 

26 

6 

3 

8 

14 

9 

2 

4 

10 

21 

13 

13 

5 

12 

5 

14 

5 

4 

12 

25 

29 

35 

307 

Apr 

17 

Monday 

i  08 

36 

25 

3 

0 

0 

6 

5 

12 

18 

6 

6 

16 

35 

16 

6 

4 

11 

39 

13 

17 

4 

0 

0 

0 

278 

Apr 

18 

Tuesday 

109 

0 

0 

0 

0 

0 

5 

5 

24 

8 

10 

23 

20 

11 

15 

9 

4 

4 

9 

4 

7 

8 

5 

16 

9 

196 

Apr 

19 

Wednesday 

110 

36 

38 

19 

15 

1 

2 

2 

12 

6 

7 

5 

13 

11 

7 

8 

7 

8 

5 

14 

3 

20 

9 

18 

9 

275 

Apr 

20 

Thursday 

111 

5 

0 

8 

3 

6 

10 

0 

6 

9 

26 

29 

44 

19 

12 

10 

3 

2 

2 

4 

6 

8 

6 

14 

25 

257 

Apr 

21 

Friday 

112 

14 

13 

16 

6 

4 

3 

1 

15 

5 

14 

50 

9 

17 

18 

5 

3 

6 

8 

47 

38 

63 

47 

30 

3 

435 

Apr 

22 

Saturday 

113 

12 

9 

18 

14 

2 

7 

8 

8 

7 

3 

4 

4 

4 

4 

1 

4 

2 

10 

4 

12 

5 

7 

25 

2 

176 

Apr 

23 

Sunday 

114 

8 

10 

2 

7 

1 

7 

6 

13 

5 

5 

10 

10 

8 

11 

4 

6 

5 

14 

8 

5 

12 

14 

10 

7 

188 

Apr 

24 

Monday 

115 

38 

55 

41 

23 

3 

8 

7 

5 

9 

10 

9 

7 

7 

7 

14 

4 

12 

12 

9 

4 

9 

8 

10 

11 

327 

Apr 

25 

Tuesday 

116 

4 

41 

48 

13 

2 

1 

3 

14 

10 

15 

7 

10 

11 

12 

10 

8 

3 

11 

6 

2 

8 

6 

11 

2 

258 

Apr 

26 

Wednesday 

117 

5 

9 

6 

2 

8 

7 

4 

9 

13 

12 

15 

19 

17 

14 

4 

9 

8 

2 

5 

11 

7 

10 

10 

5 

211 

Apr 

27 

Thursday 

118 

4 

4 

6 

2 

3 

3 

2 

5 

14 

10 

9 

22 

7 

11 

10 

6 

6 

12 

2 

10 

10 

10 

12 

6 

186 

Apr 

28 

Friday 

119 

4 

9 

3 

8 

2 

5 

10 

6 

13 

7 

10 

19 

3 

7 

7 

10 

2 

1 

7 

3 

8 

9 

6 

3 

162 

Apr 

29 

Saturday 

120 

12 

10 

14 

6 

5 

18 

7 

0 

12 

11 

4 

1 

2 

4 

1 

8 

13 

1 

8 

1 

8 

i 

7 

12 

166 

Apr 

30 

Sunday 

121 

20 

14 

ie 

4 

10 

6 

4 

1 

7 

4 

1 

0 

4 

1 

6 

3 

5 

5 

8 

6 

1 

3 

12 

4 

147 

May 

01 

Monday 

F:rste 

122 

3 

3 

2 

2 

6 

6 

8 

4 

3 

14 

6 

4 

2 

3 

1 

4 

1 

3 

9 

1 

5 

10 

12 

2 

114 

May 

02 

Tuesday 

123 

5 

3 

4 

7 

3 

8 

0 

1 

4 

6 

11 

11 

4 

7 

10 

7 

5 

11 

8 

4 

8 

5 

10 

3 

145 

May 

03 

Wednesday 

124 

5 

4 

5 

6 

6 

11 

5 

6 

11 

9 

11 

19 

9 

10 

7 

5 

4 

1 

12 

5 

5 

2 

17 

6 

181 

May 

04 

Thursday 

Kristi 

125 

3 

3 

8 

3 

3 

12 

7 

8 

5 

9 

16 

9 

12 

7 

15 

5 

4 

18 

19 

4 

10 

5 

16 

3 

204 

May 

05 

Friday 

126 

7 

4 

6 

8 

11 

0 

7 

4 

7 

13 

21 

19 

9 

1 

3 

4 

2 

1 

7 

5 

8 

9 

11 

3 

170 

May 

06 

Saturday 

127 

10 

5 

5 

8 

2 

2 

4 

5 

5 

9 

9 

5 

8 

6 

3 

19 

28 

26 

7 

22 

11 

6 

25 

10 

240 

May 

07 

Sunday 

128 

3 

2 

8 

8 

2 

l 

5 

8 

10 

16 

6 

17 

6 

9 

20 

5 

4 

7 

8 

4 

3 

13 

22 

7 

194 

May 

08 

Monday 

129 

1 

0 

5 

1 

1 

3 

7 

2 

1 

2 

16 

6 

11 

5 

5 

4 

1 

3 

3 

5 

6 

8 

19 

6 

121 

May 

09 

Tuesday 

130 

5 

2 

3 

8 

2 

2 

5 

7 

11 

16 

34 

30 

15 

21 

22 

19 

20 

21 

8 

14 

12 

9 

30 

3 

319 

May 

10 

Wednesday 

131 

5 

13 

5 

22 

8 

32 

21 

15 

16 

25 

36 

82 

11 

26 

22 

30 

19 

33 

10 

16 

24 

7 

24 

10 

512 

May 

11 

Thursday 

132 

17 

7 

5 

8 

3 

10 

15 

14 

29 

16 

18 

47 

23 

19 

27 

12 

18 

8 

6 

12 

13 

16 

24 

9 

376 

May 

12 

Friday 

133 

10 

3 

9 

13 

8 

26 

9 

11 

24 

20 

49 

15 

10 

15 

16 

13 

2 

1 

6 

13 

16 

7 

8 

4 

3oe 

May 

13 

Saturday 

134 

5 

21 

5 

0 

3 

8 

7 

7 

12 

22 

17 

15 

16 

5 

11 

16 

11 

9 

8 

7 

7 

4 

15 

6 

237 

May 

14 

Sunday 

F; rste 

135 

7 

6 

1 

2 

11 

29 

6 

7 

8 

10 

7 

32 

23 

14 

31 

20 

8 

0 

0 

0 

0 

0 

0 

0 

222 

May 

15 

Monday 

Annen 

136 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

36 

38 

8 

20 

11 

7 

14 

4 

8 

7 

14 

18 

1 2 

197 

May 

It 

Tuesday 

137 

1 

9 

10 

12 

3 

10 

7 

15 

22 

15 

11 

32 

19 

17 

23 

21 

31 

25 

18 

2 

8 

3 

16 

6 

336 

May 

17 

Wednesday 

Sytten 

138 

8 

13 

13 

25 

13 

21 

24 

15 

0 

0 

0 

0 

0 

0 

0 

0 

19 

53 

24 

13 

21 

12 

24 

5 

303 

May 

1  8 

Thursday 

139 

7 

5 

15 

12 

2 

12 

24 

15 

19 

18 

22 

33 

23 

23 

24 

5 

17 

9 

9 

19 

12 

2 

19 

3 

349 

May 

1  9 

Fr  iday 

140 

6 

17 

3 

10 

9 

18 

8 

13 

12 

16 

21 

26 

21 

22 

13 

17 

16 

6 

19 

14 

13 

4 

19 

2 

325 

May 

20 

Saturday 

141 

2 

2 

6 

2 

11 

4 

8 

10 

2 

1 

5 

6 

11 

9 

9 

24 

11 

6 

36 

14 

13 

5 

27 

2 

226 

May 

21 

Sunday 

142 

9 

9 

26 

25 

21 

29 

18 

28 

31 

21 

36 

32 

37 

20 

14 

9 

9 

7 

16 

15 

1  1 

Q 

17 

3 

452 

May 

22 

Monday 

143 

4 

8 

9 

10 

5 

14 

13 

22 

20 

36 

36 

57 

48 

40 

32 

31 

24 

15 

19 

17 

14 

12 

t  f 

17 

532 

May 

23 

Tuesday 

144 

10 

7 

12 

12 

2 

43 

30 

31 

27 

20 

27 

29 

31 

25 

16 

20 

19 

24 

15 

13 

8 

13 

21 

6 

461 

May 

24 

Wednesday 

145 

10 

12 

9 

11 

4 

39 

30 

30 

39 

35 

28 

46 

29 

39 

26 

16 

9 

23 

19 

17 

9 

6 

19 

13 

518 

May 

25 

Thursday 

146 

7 

12 

14 

3 

24 

25 

15 

39 

26 

46 

34 

35 

40 

30 

20 

17 

25 

38 

11 

6 

7 

9 

14 

6 

503 

May 

26 

Friday 

Table  3.4.2  Daily  and  hourly  distribution  of  ARCESS  detections.  For  each 
day  is  shown  the  number  of  detections  within  each  hour  of  the  day.  and  the 
number  of  detections  for  that  day.  The  statistics  at  the  end  of  the  table  give 
total  number  of  detections  for  each  hour  and  the  total  sum  of  detections  during 
the  period.  The  averages  show  number  of  processed  days,  hourly  distribution 
and  average  per  processed  day. 
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147 

8 

7 

6 

10 

15 

8 

7 

15 

16 

16 

26 

44 

16 

16 

19 

20 

8 

4 

16 

9 

29 

8 

19 

3 

34  5 

May 

27 

Saturday 

148 

13 

4 

7 

11 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Aug 

27 

Sunday 

2  40 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Aug 

28 

Monday 

2  T  1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Aug 

29 

Tuesday 

242 

0 

0 

0 

0 

0 

0 

0 

0 

30 

38 

35 

39 

26 

25 

30 

43 

19 

25 

10 

7 

0 

0 

0 

0 

327 

Aug 

30 

Wednesday 

243 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Aug 

31 

Thursday 

244 

3 

12 

16 

4 

3 

24 

25 

51 

34 

44 

36 

24 

39 

34 

39 

28 

33 

17 

11 

10 

12 

10 

14 

5 

528 

Sep 

01 

Friday 

245 

3 

8 

5 

2 

5 

12 

19 

21 

26 

53 

37 

31 

22 

3 

42 

22 

11 

11 

9 

6 

8 

5 

13 

1 

375 

Sep 

02 

Saturday 

246 

6 

1 

5 

4 

4 

9 

10 

8 

13 

19 

14 

9 

12 

10 

15 

25 

2 

31 

32 

11 

17 

16 

22 

3 

298 

Sep 

03 

Sunday 

247 

12 

3 

1 

3 

16 

37 

18 

18 

36 

36 

32 

23 

55 

12 

23 

24 

26 

13 

16 

13 

15 

6 

23 

17 

478 

Sep 

04 

Monday 

248 

5 

4 

8 

6 

9 

9 

25 

37 

21 

27 

20 

50 

27 

27 

25 

20 

9 

17 

13 

20 

5 

IP 

13 

2 

416 

Sep 

05 

Tuesday 

249 

21 

1 

15 

7 

13 

21 

26 

14 

14 

0 

0 

33 

1  7 

16 

10 

25 

28 

13 

13 

15 

2 

13 

20 

6 

348 

Sep 

06 

Wednesday 

250 

23 

0 

3 

0 

0 

26 

18 

28 

28 

30 

35 

48 

41 

45 

21 

5 

27 

26 

20 

7 

7 

13 

13 

4 

473 

Sep 

07 

Thursday 

251 

2 

13 

10 

14 

15 

10 

9 

46 

37 

38 

35 

41 

7 

0 

0 

0 

0 

15 

27 

14 

14 

7 

10 

3 

367 

Sep 

08 

Friday 

252 

5 

11 

7 

11 

4 

25 

5 

11 

16 

13 

21 

21 

17 

13 

14 

4 

4 

16 

4 

9 

6 

3 

15 

k 

260 

Sep 

09 

Saturday 

253 

10 

1 

4 

12 

e 

9 

5 

5 

17 

19 

14 

7 

14 

10 

19 

4 

20 

21 

18 

10 

14 

38 

27 

4 

310 

Sep 

10 

Sunday 

254 

i  i 

3 

6 

2 

22 

12 

22 

14 

34 

19 

24 

32 

23 

20 

18 

15 

12 

17 

24 

14 

6 

13 

17 

15 

395 

Sep 

11 

Monday 

255 

1 

3 

9 

16 

20 

24 

8 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

81 

Sep 

12 

Tuesday 

256 

f) 

0 

0 

0 

0 

0 

9 

34 

21 

0 

24 

31 

39 

17 

15 

24 

29 

6 

15 

19 

2 

18 

21 

6 

330 

Sep 

13 

Wednesday 

257 

1 

9 

5 

28 

23 

22 

13 

5 

0 

0 

0 

0 

0 

0 

0 

4 

16 

16 

10 

27 

15 

2 

19 

6 

221 

Sep 

14 

Thursday 

258 

16 

20 

8 

9 

6 

26 

22 

32 

10 

41 

25 

20 

21 

19 

18 

20 

26 

13 

23 

26 

6 

15 

12 

7 

44  1 

Sep 

15 

Friday 

259 

0 

8 

23 

17 

10 

19 

9 

13 

12 

13 

29 

29 

14 

14 

27 

23 

12 

4 

13 

4 

4 

2 

13 

10 

331 

Sep 

16 

Saturday 

260 

g 

35 

6 

1 

0 

0 

0 

0 

0 

0 

0 

3 

16 

15 

21 

10 

9 

13 

27 

21 

20 

6 

12 

3 

226 

Sep 

17 

Sunday 

261 

4 

4 

3 

22 

14 

19 

27 

10 

20 

12 

16 

46 

12 

14 

34 

10 

25 

21 

23 

16 

17 

25 

13 

13 

425 

Sep 

18 

Monday 

262 

3 

13 

4 

11 

7 

13 

20 

19 

25 

25 

10 

45 

15 

25 

9 

18 

23 

23 

9 

13 

4 

20 

21 

5 

388 

Sep 

19 

Tuesday 

263 

7 

3 

2 

14 

3 

25 

25 

13 

9 

17 

18 

30 

19 

17 

12 

9 

9 

12 

12 

10 

4 

14 

14 

7 

305 

Sep 

20 

Wednesday 

264 

12 

7 

6 

10 

12 

15 

7 

26 

26 

19 

40 

30 

27 

27 

21 

26 

13 

12 

16 

6 

9 

8 

11 

3 

339 

Sep 

21 

Thursday 

/.  6  5 

s 

2 

7 

5 

13 

16 

5 

17 

23 

14 

21 

18 

20 

20 

16 

18 

22 

12 

8 

19 

12 

15 

10 

6 

324 

Sep 

22 

Friday 

266 

1  0 

4 

13 

3 

1 

18 

3 

21 

18 

27 

22 

18 

9 

16 

24 

11 

16 

6 

21 

10 

26 

1 

21 

2 

321 

Sep 

23 

Saturday 

Table  3.4.2  (cont.) 


28 


267 

6 

3 

5 

u 

11 

36 

14 

6 

6 

12 

5 

0 

20 

35 

69 

29 

12 

11 

14 

10 

21 

9 

25 

15 

385 

Sep 

24 

Sunday 

268 

6 

5 

15 

35 

76 

15 

20 

19 

20 

10 

16 

26137106 

82 

90 

15 

10 

15 

10 

16 

7 

40 

61 

852 

Sep 

25 

Monday 

269 

97 

41 

10 

15128 

30 

15 

11 

17 

9114 

93 

64111 

13 

28 

31 

8 

18 

12 

11 

4 

14 

14 

908 

Sep 

26 

Tuesday 

270 

4 

10 

2 

14 

19 

14 

20 

25 

7 

13 

12 

24 

35 

5 

5 

12 

10 

14 

0 

11 

8 

6 

6 

14 

290 

Sep 

27 

Wednesday 

271 

7 

14 

3 

6 

6 

11 

17 

22 

9 

19 

12 

25 

30 

23 

47 

59 

79 

10 

10 

11 

12 

6 

15 

14 

467 

Sep 

28 

Thursday 

272 

4 

14 

1 

8 

4 

4 

19 

12 

22 

16 

9 

23 

42 

29 

22 

30 

23 

7 

4 

7 

12 

6 

7 

28 

353 

Sep 

29 

Friday 

273 

6 

12 

4 

7 

10 

6 

3 

9 

18 

7 

19 

18 

22 

24 

10 

19 

4 

8 

6 

4 

7 

6 

6 

12 

247 

Sep 

30 

Saturday 

FRS 

00 

01 

02 

03 

04 

05 

06 

07 

08 

09 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

Sum 

1245 

1502 

2361 

2902 

3344 

4562 

2894 

2407 

2315 

1832 

1594 

1161 

1243 

1210 

1640 

2510 

3032 

3908 

3597 

3027 

2166 

1985 

1591 

2347 

56375 

Total 

sum 

162 

8 

8 

7 

9 

10 

15 

15 

18 

19 

21 

24 

28 

22 

18 

19 

15 

13 

14 

12 

11 

10 

10 

14 

7 

348 

Total 

average 

126 

7 

7 

6 

9 

10 

14 

16 

19 

19 

21 

23 

28 

23 

18 

18 

15 

13 

14 

10 

10 

8 

9 

13 

7 

338 

Average  workdays 

57 

7 

7 

7 

6 

6 

10 

8 

10 

12 

13 

18 

19 

11 

11 

12 

10 

8 

9 

12 

9 

10 

8 

12 

5 

241 

Average  weekends 

Table  3.4.2  (  cont.) 


29 


4  Improvements  and  modifications 


4.1  NORSAR 

Data  acquisition 

No  modification  has  been  made  to  the  JfORSAR  data  acquisition  system. 
The  MODCOMP  system,  which  communicates  with  the  NORSAR  subarrays 
and  the  IBM  acquisition  and  processing  system,  has  been  very  unstable. 

Due  to  the  many  old  parts  in  the  system  data  flow,  we  have  found  no  way 
to  make  partial  replacements,  like  say  replacing  the  MODCOMP.  The  only 
rational  upgrade/replacement  seems  to  be  to  replace  all  systems  from  A/D 
converters  (today's  SLEM)  to  central  data  recording. 

An  upgrade  of  the  NORSAR  system  should  include  digital  data  transmis¬ 
sion  from  seismometer  to  subarray  vault,  and  the  resulting  dynamic  range 
should  be  comparable  to  that  of  NORESS/GERESS  (at  least  21  bit  dynamic 
range).  The  upgrade  should  include  broadband  seismometers  and  local  record¬ 
ing  at  the  subarray  vault. 

Depending  on  which  system  is  chosen,  cables  between  subarrav  vaults  and 
seismometers  may  have  to  be  replaced,  due  to  power  requirements  and/or  data 
transmission  technique. 

To  use  the  NORESS/ARCESS  technology,  seven  such  systems  would  be 
required.  Moreover,  there  may  be  a  problem  with  cabling  and  power.  Another 
problem  may  be  time  synchronization  for  array  beamforming. 

The  Teledyne  GERESS  system  has  not  yet  been  tested,  but  may  also  have 
too  large  power  requirements  to  use  existing  cables. 

The  new  FINESA  system  would  make  some  improvement  with  less  cost. 
However,  data  would  still  be  limited  to  the  dynamic  range  of  the  analog  signals 
on  the  buried  cables. 

The  Nanometrics  system  used  at  Yellowknife  and  Hagfors  in  Sweden  has 
features  that  seem  to  fit  very  well  for  an  upgrade  of  the  NORSAR  array.  With 
acceptable  dynamic  range  and  16  bit  resolution,  the  system  may  use  existing 
cabling  to  meet  all  other  requirements.  This  system  has  local  recording. 

Detection  processing 

The  NORSAR  detection  processor  has  been  running  satisfactorily  on  the 
IBM  during  this  reporting  period.  Processing  has  been  moved  to  a  more 
powerful  IBM  computer,  that  became  available  after  moving  NORESS  to  a 
SUN-system . 
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Detection  statistics  are  given  in  chapter  2. 

Event  processing 

There  are  no  changes  in  the  routine  processing  of  NORSAR  events,  us¬ 
ing  the  IBM  system.  Full  event  processing  on  the  SUN  system  is  not  yet 
implemented. 


4.2  NORESS 
Detection  processing 

The  routine  detection  processing  of  NORESS  is  running  satisfactory  on  a 
SUN-3/280.  From  day  091:00  through  day  103:08,  the  RONAPP  beam  set  of 
17  coherent  beams  and  3  incoherent  beams  was  used.  The  detection  processor 
recipe  for  this  beam  set  was  described  in  NORSAR  Sci.  Rep.  No.  2-88/89. 
The  beam  set  is  given  in  Table  4.2.1. 

From  1989.  day  103:09,  we  have  used  a  recipe  with  66  coherent  and  10 
incoherent  beams,  which  nearly  doubled  the  number  of  detections  per  day.  The 
beam  list  is  shown  in  Table  4.2.2.  This  processing  includes  STA/LTA  detection 
processing  only,  as  the  computer  intensive  f-k  analysis  for  each  detection  is 
performed  offline  as  part  of  the  event  processing.  This  detector  gives  two  files 
per  day,  e.g.,  NRS89091.DPX,  indicating  detection  processor  report  (.DPX) 
for  NORESS  (NRS),  day  091,  1989,  and  NRS89091.REP,  giving  reports  on 
the  occurrence  of  spikes  and  otherwise  data  of  bad  quality. 

We  have  tested  the  SUN-3/280  array  data  acquisition  system  with  vari¬ 
able  processing  load  and  found  that  STA/LTA  processing  may  be  performed 
without  disturbing  acquisition  or  data  serving.  However,  the  added  load  of  f-k 
analysis  for  this  high  number  of  detections  is  not  advisable. 

Event  processing.  Phase  estimation 

The  NORESS  event  processing  is  performed  in  two  steps.  In  step  1, 
all  detections  are  subjected  to  broadband  f-k,  polarization,  and  onset  anal¬ 
ysis,  called  phase  estimation.  This  analysis  results  in  one  file  per  day,  e.g., 
NRS89091.FKX,  indicating  f-k  analysis  report  (.FKX)  for  NORESS  (NRS), 
day  091,  1989.  The  report  files  give  onset  time,  detection  beam,  SNR,  ve¬ 
locity,  preliminary  phase  name,  azimuth,  coherency,  quality  of  f-k,  frequency, 
amplitude,  detection  STA  and  polarization  results. 

This  processing  may  be  performed  on  any  of  the  SUN  workstations. 
Event  processing.  Plot  and  Epicenter  determination 

In  step  two  of  the  automatic  event  processing,  all  phases  from  one  array  are 
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analyzed  for  association  to  events,  as  described  in  NORSAR  Sci.  Rep.  2-88/89. 
The  location  analysis  results  in  one  file  per  day.  e.g.,  NRS89091.EPX.  indi¬ 
cating  event  analysis  report  (.EPX)  for  XORESS  (MRS),  day  091,  1989.  All 
‘interesting'  phases  are  plotted  on  an  Imagen  postscript  laser  printer.  These 
plots  also  give  event  location,  origin  time  and  magnitude.  Where  a  minimum 
of  two  phases  (at  least  one  P  and  one  S)  have  been  associated,  we  have  located 
the  event  with  the  TTAZLOC  procedure.  _ 


4.3  ARCESS 

Detection  processing 

The  routine  detection  processing  of  ARCESS  is  running  satisfactory  on  a 
SUN-3/280.  From  day  091:00  through  day  130:09.  the  RONAPP  beam  set  of 
17  coherent  beams  and  3  incoherent  beams  was  used.  The  detection  processor 
recipe  for  this  beam  set  was  described  in  NORSAR  Sci.  Rep.  No.  2-88/89. 
The  beam  set  is  given  in  Table  4.3.1.  From  day  1989-130:09,  we  have  used 
a  recipe  with  66  coherent  and  10  incoherent  beams,  which  nearly  doubled 
the  number  of  detections  per  day.  The  beam  list  is  shown  in  Table  4.3.2. 
This  processing  include  STA/LTA  detection  processing  only,  as  the  computer 
intensive  f-k  analysis  for  each  detection  is  performed  offline  as  part  of  the 
event  processing.  This  detector  give  two  file  per  day.  E.g..  FRS89091.DPX, 
indicating  detection  processor  report  (.DPX)  for  ARCESS  (FRS.ARC),  day- 
091,  1989,  and  FRS89091.REP,  giving  reports  on  the  occurrence  of  spikes  and 
otherwise  data  of  bad  quality 

Event  processing.  Phase  estimation 

The  ARCESS  event  processing  is  identical  to  that  of  NORESS.  This  pro¬ 
cessing  results  in  one  file  per  day,  e.g..  FRS89091.FKX.  indicating  f-k  analysis 
report  (,FKX)  for  ARCESS  (FRS.ARC),  day  091,  1989. 

Event  processing.  Plot  and  Epicenter  determination 

In  step  two  of  the  automatic  event  processing,  all  phases  from  one  array  are 
analyzed  for  association  to  events,  as  described  in  NORSAR  Sci.  Rep.  No.  2- 
88/89.  The  location  analysis  results  in  one  file  per  day.  e.g.,  FRS89091.EPX. 
indicating  event  analysis  report  (  EPX)  for  ARCESS  (FRS.ARC),  day  091. 
1989.  All  ‘interesting’  phases  are  plotted  on  an  Imagen  postscript  laser  printer. 
These  plots  also  give  event  location,  origin  time  and  magnitude.  Where  a 
minimum  of  two  phases  (at  least  one  P  and  one  S)  have  been  associated,  we 
have  located  the  event  with  the  TTAZLOC  procedure. 
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J.  Fyen 


BEAM  Velocity  Azimuth 

Filter 

band 

Threshold 

Configuration 

NR01 

99999.9 

0.0 

1.0 

- 

3.0 

4.0 

17 

CD 

NR02 

99999.9 

0.0 

1.5 

- 

3.5 

4.0 

17 

CD 

NR03 

99999.9 

0.0 

2.0 

- 

4.0 

4.0 

22 

BCD 

NR04 

99999.9 

0.0 

2.5 

- 

4.5 

4.0 

13 

BC 

NR05 

99999.9 

0.0 

3.0 

- 

5.0 

4.0 

13 

BC 

NR06 

99999.9 

0.0 

4 . 0 

- 

8.0 

5.0 

9 

AB 

NR07 

99999.9 

0.0 

8. 

- 

16.0 

5.0 

9 

AB 

NR08 

14.3 

0.0 

2.0 

- 

4.0 

4.0 

22 

BCD 

NR09 

14.3 

90.0 

2.0 

- 

4.0 

4.0 

22 

BCD 

NR10 

14.3 

180.0 

2.0 

- 

4.0 

4.0 

22 

BCD 

NR11 

14.3 

15.0 

2.5 

- 

4.5 

4.0 

22 

BCD 

NR12 

34.3 

75.0 

2.5 

- 

4.5 

4.0 

22 

BCD 

NR13 

14.3 

135.0 

2.5 

- 

4.5 

4.0 

22 

BCD 

NR14 

14.3 

25.0 

3.0 

- 

5.0 

4.0 

22 

BCD 

NR15 

14.3 

75.0 

3.0 

- 

5.0 

4.0 

22 

BCD 

NR16 

14.3 

125.0 

3.0 

- 

5.0 

4.0 

22 

BCD 

NR17 

99999.9 

0.0 

2.0 

- 

4.0 

4.0 

13 

BC 

NI18 

99999.9 

0.0 

1.0 

- 

2.0 

2.5 

17 

C 

NI19 

99999.9 

0.0 

2.0 

- 

3.0 

2.5 

22 

C 

NI20 

99999.9 

0.0 

2.0 

- 

4.0 

2.1 

13 

D 

IN18 

99999.9 

0.0 

1.0 

- 

2.0 

2.5 

17 

C 

I N 1 9 

99999.9 

0.0 

2.0 

- 

3.0 

2.5 

22 

C 

TN20 

99999.9 

0.0 

2.0 

- 

4.0 

2.1 

13 

D 

Table  4.2.1  NORESS  ‘RONAP'  beam  table.  The  table  shows  the  name 
of  beam,  velocity  (km/sec),  azimuth  (degrees),  filter  band  (Hz).  STA/LTA 
threshold,  and  configuration.  The  configuration  is  described  with  number  of 
sensors  and  A.B.C  or  D-rings.  Center  AO  is  always  within  the  configuration. 
N 1 18  -  NI20  are  incoherent  beams  using  vertical  components.  IN  18  -  IN20 
are  identical  to  NI18  -  NI20.  This  beam  set  was  used  up  to  1989.  day  103:08. 
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BEAM  Velocity  Azimuth  Filter  band  Threshold  Configuration 


N011 

99999.9 

0.0 

0.5 

N021 

99999.9 

0.0 

1.0 

N031 

99999.9 

0.0 

1.5 

N032 

11.0 

30.0 

1.5 

N033 

11.0 

90.0 

1.5 

N034 

11.0 

150.0 

1.5 

N035 

11.0 

210.0 

1.5 

N036 

11.0 

270.0 

1.5 

N037 

11.0 

330.0 

1.5 

N038 

15.8 

80.0 

1.5 

N039 

10.0 

30.0 

1.5 

NO  41 

99999.9 

0.0 

2.0 

N042 

10.2 

30.0 

2.0 

N04  3 

10.2 

90.0 

2.0 

N044 

10.2 

150.0 

2.0 

N045 

10.2 

210.0 

2.0 

N046 

10.2 

270.0 

2.0 

N047 

10.2 

330.0 

2.0 

N04  8 

15.8 

80.0 

2.0 

N049 

10.0 

30.0 

2.0 

N051 

99999.9 

0.0 

2.5 

N052 

8.90 

30.0 

2.5 

N053 

8.90 

90.0 

2.5 

N054 

8.90 

150.0 

2.5 

N055 

8.90 

210.0 

2.5 

N056 

8.90 

270.0 

2.5 

N057 

8.90 

330.0 

2.5 

N058 

15.8 

80.0 

2.5 

N059 

10.0 

30.0 

2.5 

N061 

99999.9 

0.0 

3.0 

N062 

10.5 

30.0 

3.0 

N063 

10.5 

90.0 

3.0 

N064 

10.5 

150.0 

3.0 

N065 

10.5 

210.0 

3.0 

N066 

10.5 

270.0 

3.0 

N067 

10.5 

330.0 

3.0 

N068 

15.8 

80.0 

3.0 

N069 

10.0 

30.0 

3.0 

N071 

99999.9 

0.0 

3.5 

N072 

11.1 

30.0 

3.5 

N073 

11.1 

90.0 

3.5 

N074 

11.1 

150.0 

3.5 

N075 

11.1 

210.0 

3.5 

1.5 

4.0 

10 

D 

- 

3.0 

4.0 

17 

CD 

— 

3.5 

4.0 

17 

CD 

3.5 

4.0 

17 

CD 

— 

3.5 

4.0 

17 

CD 

- 

3.5 

4.0 

17 

CD 

- 

3.5 

4.0 

17 

CD 

- 

3.5 

4.0 

17 

CD 

- 

3.5 

4.0 

17 

CD 

3.5 

3.5 

17 

CD 

~ 

3.5 

3.5 

17 

CD 

4.0 

4.0 

17 

CD 

4.0 

4.0 

17 

CD 

— 

4.0 

4.0 

17 

CD 

— 

4.0 

4.0 

17 

CD 

— 

4.0 

4.0 

i 

CD 

— 

4.0 

4.0 

17 

CD 

— 

4.0 

4.0 

17 

CD 

- 

4.0 

3.5 

17 

CD 

- 

4.0 

3.5 

17 

CD 

4.5 

4.0 

22 

BCD 

- 

4.5 

4.0 

22 

BCD 

- 

4.5 

4.0 

22 

BCD 

- 

4.5 

4.0 

22 

BCD 

4.5 

4.0 

22 

BCD 

- 

4.5 

4.0 

22 

BCD 

- 

4.5 

4.0 

22 

BCD 

- 

4.5 

3.5 

22 

BCD 

4.5 

3.5 

22 

BCD 

- 

5.0 

4.0 

22 

BCD 

— 

5.0 

4.0 

22 

BCD 

— 

5.0 

4.0 

22 

BCD 

- 

5.0 

4.0 

22 

BCD 

— 

5.0 

4.0 

22 

BCD 

— 

5.0 

4.0 

22 

BCD 

- 

5.0 

4.0 

22 

BCD 

5.0 

3.5 

22 

BCD 

- 

5.0 

3.5 

22 

BCD 

5.5 

4.0 

13 

BC 

- 

5.5 

4.0 

13 

BC 

- 

5.5 

4.0 

13 

BC 

- 

5.5 

4.0 

13 

BC 

5.5 

4.0 

13 

BC 

Table  4.2.2  NORESS  beam  table.  The  table  shows  the  name  of  beam,  ve¬ 
locity  (km/sec),  azimuth  (degrees),  filter  band  (Hz),  STA/LTA  threshold,  and 
configuration.  The  configuration  is  described  with  number  of  sensors  and 
A,B,C  or  D-rings.  Center  AO  is  always  within  the  configuration.  NH01  - 
NH04  are  incoherent  beams  using  horizontal  components  of  the  A0,C2,C4  and 
CT  three-component  stations.  NV01  -  NV04  are  incoherent  beams  using  ver¬ 
tical  components.  This  beam  set  has  been  in  operation  since  1989.  day  103:08. 
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N076 
N077 
N081 
N082 
N083 
N084 
N085 
N086 
N087 
N091 
N092 
N093 
NO  94 

N096 
N097 
N101 
N102 
N103 
N104 
N105 
N106 
N107  9.9 

NHOl  99999.9 

NH02  99999.9 

NH03  99999.9 

NH04  99999.9 

NVOl  99999.9 

NV02  99999.9 

NV03  99999.9 

NV04  99999.9 

NV05  99999.9 

NV06  99999.9 


5.5 

5.5 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
8.0 
10. 0 

5.0  -  10.0 
5.0  --  10 . 0 
5.0  -  10.0 
5.0  -  10.0 
5.0  -  10.0 
5.0  -  10.0 
8.0  -  16.0 
8.0  -  16.0 
8.0  -  16.0 
8.0  -  16.0 
8.0  -  16.0 
8.0  -  16.0 
8.0  -  16.0 
2.0  -  4.0 

3.5  -  5.5 

5.0  -  10.0 
8.0  -  16.0 
0.5  -  1.5 

0.0  1.0  ~  2.0 

0.0  1.5  -  2.5 

0.0  2.0  -  3.0 

0.0  2.0  -  4.0 

0.0  3.5  -  5.5 


4.0 

13 

BC 

4.0 

13 

BC 

4.0 

13 

BC 

4.0 

13 

BC 

4.0 

13 

BC 

4.0 

13 

BC 

4.0 

13 

BC 

4.0 

13 

BC 

4.0 

13 

BC 

4.5 

13 

BC 

4.5 

13 

BC 

4.5 

13 

BC 

4.5 

13 

BC 

4.5 

1 

DC 

4.5 

13 

BC 

4.5 

13 

BC 

4.5 

9 

AB 

4.5 

9 

AB 

4.5 

9 

AB 

4.5 

9 

AB 

4.5 

9 

AB 

4.5 

9 

AB 

4.5 

9 

AB 

2.4 

22 

AC_ne 

2.4 

13 

AC_ne 

2.4 

13 

AC__ne 

2.5 

9 

AC  ne 

2.5 

17 

D 

2.5 

17 

C 

2.5 

17 

C 

2.5 

22 

C 

2.4 

22 

C 

2.4 

13 

C 

11.1 

11.1 

99999.9 

9.5 

9.5 

9.5 

9.5 

9.5 

9.5 

99999.9 

10.5 

10.5 

10.5 

10.5 

10.5 

10.5 

99999.9 

9.9 

9.9 

9.9 

9.9 

9.9 


270.0 

330.0 

0.0 

30.0 

90.0 

150.0 

210.0 

270.0 

330.0 

0.0 

30.0 

90.0 

150.0 

210.0 

270.0 

330.0 

0.0 

30.0 

90.0 

150.0 

210.0 

270.0 

330.0 

0.0 

0.0 

0.0 

0.0 

0.0 


3.5 

3.5 

4.0 


0 

0 

0 

0 

,0 

,0 

.0 


Table  4.2.2  (cont.) 
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BEAM  Velocity  Azimuth  Filter  band  Threshold  Configuration 


FR01 

99999.9 

0.0 

1.0 

FRO  2 

99999.9 

0.0 

1.5 

FRO  3 

99999.9 

0.0 

2.0 

FR04 

99999.9 

0.0 

2.5 

FRO  5 

99999.9 

0.0 

3.0 

FR06 

99999.9 

0.0 

4.0 

FR07 

99999.9 

0.0 

8.0 

FRO  8 

14.3 

0.0 

2.0 

FRO  9 

14.3 

90.0 

2.0 

FRIO 

14.3 

180.0 

2.0 

FR11 

14.3 

15.0 

2.5 

FR12 

14.3 

75.0 

2.5 

FR13 

14.3 

135.0 

2.5 

FR14 

14.3 

25.0 

3.0 

FR15 

14.3 

75.0 

3.0 

FR16 

14.3 

125.0 

3.0 

FR17 

99999.9 

0.0 

2.0 

IF18 

99999.9 

0.0 

1.0 

IF19 

99999.9 

0.0 

2.0 

IF20 

99999.9 

0.0 

2.0 

FI18 

99999.9 

0.0 

1.0 

FI19 

99999.9 

0.0 

2.0 

FI20 

99999.9 

0.0 

2.0 

3.0 

4.0 

17 

CD 

3.5 

4.0 

17 

CD 

4.0 

4.0 

22 

BCD 

4.5 

4.0 

13 

BC 

5.0 

4.0 

13 

BC 

8.0 

5.0 

9 

AB 

16.0 

5.0 

9 

AB 

4.0 

4.0 

22 

BCD 

4.0 

4.0 

22 

BCD 

4.0 

4.0 

22 

BCD 

4.5 

4.0 

22 

BCD 

4.5 

4.0 

22 

BCD 

4.5 

4.0 

22 

BCD 

5.0 

4.0 

22 

BCD 

5.0 

4.0 

22 

BCD 

5.0 

4.0 

22 

BCD 

4.0 

4.0 

13 

BC 

2.0 

2.5 

17 

C 

3.0 

2.5 

22 

C 

4.0 

2.1 

13 

D 

2.0 

2.5 

17 

C 

3.0 

2.5 

22 

C 

4.0 

2.1 

13 

D 

Table  4.3.1  ARCESS  ‘RONAP’  beam  table.  The  table  shows  the  name 
of  beam,  velocity  (km/sec),  azimuth  (degrees),  filter  band  (Hz).  STA/LTA 
threshold,  and  configuration.  The  configuration  is  described  with  number  of 
sensors  and  A.B,C  or  D-rings.  Center  AO  is  always  within  the  configuration. 
FI  18  -  FI20  are  incoherent  beams  usine  vertical  components.  IF  1 8  -  TF'20  are 
identical  to  FI18  -  F120.  This  beam  set  was  used  up  to  1989.  day  103:09. 
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Filter  band  Threshold  Configuration 


BEAM  Velocity  Azimuth 


F011 

99999.9 

0.0 

F021 

99999.9 

0.0 

F031 

99999.9 

0.0 

F032 

11.0 

30.0 

F033 

11.0 

90.0 

F034 

11.0 

150.0 

F035 

11.0 

210.0 

F036 

11.0 

270.0 

F037 

11.0 

330.0 

F038 

15.8 

80.0 

F039 

10.0 

30.0 

F04 1 

99999.9 

0.0 

F042 

10.2 

30.0 

F043 

10.2 

90.0 

F044 

10.2 

150.0 

F04  5 

10.2 

210.0 

F046 

10.2 

270.0 

F047 

10.2 

330.0 

F04  8 

15.8 

80.0 

F049 

10.0 

30.0 

F051 

99999.9 

0.0 

F052 

8.90 

30.0 

F053 

8.90 

90.0 

F054 

8.90 

150.0 

F055 

8.90 

210.0 

F056 

8.90 

270.0 

F057 

8.90 

330.0 

F058 

15.8 

80.0 

F059 

10.0 

30.0 

F061 

99999.9 

0.0 

F062 

10.5 

30.0 

F063 

10.5 

90.0 

F064 

10.5 

150.0 

F065 

10.5 

210.0 

F066 

10.5 

270.0 

F067 

10.5 

330.0 

F068 

15.8 

80.0 

F069 

10.0 

30.0 

F071 

99999.9 

0.0 

F072 

11.1 

30.0 

F073 

11.1 

90.0 

F074 

11.1 

150.0 

F075 

11.1 

210.0 

0.5  - 

1.5 

4.0 

1.0  - 

3.0 

4.0 

1.5  - 

3.5 

4.0 

1.5  - 

3.5 

4.0 

1.5  - 

3.5 

4.0 

1.5  -  . 

3.5 

4.0 

1.5  - 

3.5 

4.0 

1.5  - 

3.5 

4.0 

1.5-  - 

3.5 

4.0 

1.5  - 

3.5 

3.5 

1.5  - 

3.5 

3.5 

2.0  - 

4.0 

4.0 

2.0  - 

4.0 

4.0 

2.0  - 

4.0 

4.0 

2.0  - 

4.0 

4.0 

2.0  - 

4.0 

4.0 

2.0  - 

4.0 

4.0 

2.0  - 

4.0 

4.0 

2.0  - 

4.0 

3.5 

2.0  - 

4.0 

3.5 

2.5  - 

4.5 

4.0 

2.5  - 

4.5 

4.0 

2.5  - 

4.5 

4.0 

2.5  - 

4.5 

4.0 

2.5  - 

4.5 

4.0 

2.5  - 

4.5 

4.0 

2.5  - 

4.5 

4.0 

2.5  - 

4.5 

3.5 

2.5  - 

4.5 

3.5 

3.0  - 

5.0 

4.0 

3.0  - 

5.0 

4.0 

3.0  - 

5.0 

4.0 

3.0  - 

5.0 

4.0 

3.0  - 

5.0 

4.0 

3.0  - 

5.0 

4.0 

3.0  - 

5.0 

4.0 

3.0  - 

5.0 

3.5 

3.0  - 

5.0 

3.5 

3.5  - 

5.5 

4.0 

3.5  - 

5.5 

4.0 

3.5  - 

5.5 

4.0 

3.5  - 

5.5 

4.0 

3.5  - 

5.5 

4.0 

10 

D 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

17 

CD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

22 

BCD 

13 

BC 

13 

BC 

13 

BC 

13 

BC 

13 

BC 

Table  4.3.2  ARCESS  beam  table.  The  table  shows  the  name  of  beam,  velocity 
(km/sec),  azimuth  (degrees),  filter  band  (Hz),  STA/LTA  threshold,  and  con¬ 
figuration.  The  configuration  is  described  with  number  of  sensors  and  A.B.C 
or  D-rings.  Center  AO  is  always  within  the  configuration.  FH01  -  FH04 
are  incoherent  beams  using  horizontal  components  of  the  A0.C2.C4  and  C7 
three-component  stations.  FV01  -  FV04  are  incoherent  beams  using  vertical 
components.  This  beam  set  has  been  in  operation  since  1989.  day  103:09. 
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F076 

11.1 

270.0 

3.5  - 

5.5 

4.0 

13 

BC 

F077 

11.1 

330.0 

3.5  - 

5.5 

4.0 

13 

BC 

F081 

99999.9 

0.0 

4.0  - 

8.0 

4.0 

13 

BC 

F082 

9.50 

30.0 

4.0  - 

8.0 

4.0 

13 

BC 

F083 

9.50 

90.0 

4.0  - 

8.0 

4.0 

13 

BC 

F084 

9.50 

150.0 

4.0  - 

8.0 

4.0 

13 

BC 

F085 

9.50 

210.0 

4.0  - 

8.0 

4.0 

13 

BC 

F086 

9.50 

270.0 

4.0  - 

8.0 

4.0 

13 

BC 

F087 

9.50 

330.0 

4.0  - 

8.0 

4.0 

13 

BC 

F091 

99999.9 

0.0 

5.0  - 

10.0 

4.5 

13 

BC 

F092 

10.5 

30.0 

5.0  - 

-10.0 

4.5 

13 

BC 

F093 

10.5 

90.0 

5.0  - 

10.0 

4.5 

13 

BC 

F094 

10.5 

150.0 

5.0  - 

10.0 

4.5 

13 

BC 

F095 

10.5 

210.0 

5.0  - 

10.0 

4.5 

13 

BC 

F096 

10.5 

270.0 

b .  0  - 

10.0 

4.5 

13 

BC 

F097 

10.5 

330.0 

5.0  - 

10.0 

4.5 

13 

BC 

F101 

99999.9 

0.0 

8.0  - 

16.0 

4.5 

9 

AB 

F102 

9.9 

30.0 

8.0  - 

16.0 

4.5 

9 

AB 

F103 

9.9 

90.0 

8.0  - 

16.0 

4.5 

9 

AB 

F104 

9.9 

150.0 

8.0  - 

16.0 

4.5 

9 

AB 

F105 

9.9 

210.0 

8.0  - 

16.0 

4.5 

9 

AB 

F106 

9.9 

270.0 

8.0  - 

16.0 

4.5 

9 

AB 

F107 

9.9 

330.0 

8.0  - 

16.0 

4.5 

9 

AB 

FH01 

99999.9 

0.0 

2.0  - 

4.0 

2.4 

22 

AC_ne 

FH02 

99999.9 

0.0 

3.5  - 

5.5 

2.4 

13 

AC_ne 

FH03 

99999.9 

0.0 

5.0  - 

10.0 

2.4 

13 

AC_ne 

FH04 

99999.9 

0.0 

8.0  - 

16.0 

2.5 

9 

AC  ne 

FV01 

99999.9 

0.0 

0.5  - 

1.5 

2.5 

17 

D 

FV02 

99999.9 

0.0 

1.0  - 

2.0 

2.5 

17 

C 

FV03 

99999.9 

0.0 

1.5  - 

2.5 

2.5 

17 

C 

FV04 

99999.9 

0.0 

2.0  - 

3.0 

2.5 

22 

C 

FV05 

99999.9 

0.0 

2.0  - 

4.0 

2.4 

22 

C 

FV06 

99999.9 

0.0 

3.5  - 

5.5 

2.4 

13 

C 

Table  4.3.2  (  cont.) 
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5  Maintenance  Activities 


5.1  Activities  in  the  field  and  at  the  Maintenance  Center 

This  section  summarizes  the  activities  in  the  field,  at  the  Maintenance  Center 
(NMC)  Hamar  and  NDPC  activities  related  to  monitoring  and  control  of  the 
NORSAR,  NORESS  and  ARCESS  arrays. 

Activities  comprise  preventive/corrective  maintenance  in  connection  with 
all  the  NORSAR  subarrays.  NORESS  and  ARCESS.  Other  activities  are  re¬ 
lated  to  testing  the  NORSAR  communications  systems. 

Cables  have  been  pointed  out  (01A.02C)  and  spliced  (04C.02B,0lB,02C). 
SLE.Ms  have  been  replaced/reset  (02B.06C).  LP  seismometer  and  related 
equipment  have  been  adjusted/rcplaced  (02B.0lB.03C). 

NORESS 

Power  and  fiber  cables  have  been  cut  and  spliced,  buried  and  rerouted 
in  three  new  trenches  (NORESS,  C3).  Hub  station  power  distribution  and 
power  supply  at  site  D3  have  been  replaced  after  lightning.  Electronics  on 
remote  sites  A1,A3,B3,B5  were  repaired.  Satellite  modem  was  replaced  and 
broadband  instrument  KS36000  installed  after  repair. 

ARCESS 

A  new  NORSAT  B  antenna  was  erected.  Ditches  and  roads  were  sown. 
Fiber  optical  links  at  all  sites  were  adjusted. 

Details  are  presented  in  Table  5.1. 
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Subarray/ 

Task 

Date 

area 

NDPC 

Daily  check  of  SP/LP  data  and  comm,  systems. 
Check  including  NORESS  and  ARCESS  data. 
Weekly  calibrations  of  SP/LP  instruments. 
Adjustment  of  MP/FP  parameters  when  outside 
tolerances  _ 

April 

01A 

In  connection  with  digging,  SP05  pointed  out 
to  a  farmer. 

19  May 

02  B 

SLEM  reset  after  power  outage.  Cable  work  SP02. 

22  May 

04C 

Cables  splicing  SP04. 

24  May 

NORESS: 

Power  and  fiber  cables  to  site  C3  cut  and  spliced. 
Afterwards  buried  in  a  new  trench 

May 

NDPC 

Daily  routine  check  of  NORSAR.  NORESS  and 

May 

ARCESS  data.  Weekly  calibration  of  SP/LP 

instruments. 

NORSAR: 

02C 

Cable  splicing  SP05. 

19.20  June 

02  B 

Cable  splicing  SP02. 

21 .22  June 

06C 

Check  of  EW  LP  seismometer. 

23  June 

04C 

Cable  splicing  SP04. 

26  June 

02B 

Replaced  MP  motor  (RCD)  EW  LP  seismometer. 

27  June 

02B 

Cable  fault  SP04. 

28  June 

06C 

Digital  Unit  replaced  (SLEM). 

29  June 

01B 

Cable  splicing  SP02. 

30  June 

ARCESS: 

Ditches  and  roads  sown. 

1-12  June 

Fiber  opitical  links  at  all  sites  adjusted. 

11.12  June 

NORESS: 

Satellite  transmitter  modem.  Hub  and  remote 
site  (C3.D5.D2.D3)  were  still  down  due  to  lack 
of  spare  parts  (from  Sandia). 

J  une 

NORSAR: 
04  C 

Visit  to  SP02.03  and  00  in  connection  with 

6  July 

01 A 

repair  after  lightning. 

Cable  work  SP02. 

24.26  July 

01 A 

Cable  work  SP04. 

27.28  July 

02C 

Cable  splicing  SP05. 

25  July 

Table  5.1  Activities  in  the  field  and  the  .NORSAR  Maintenance  Center,  in¬ 
cluding  NDPC  activities  related  to  the  NORSAR  array.  1  April  -  AO  September 
1989. 
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Subarrav/  Task  Date 

area 


ARCESS: 

NORESS: 

NDPC: 

Array  visited  by  P.VV.  Larsen.  Satellite 
modem  was  replaced. 

Hub  station  power  distribution  replaced. 

Power  supply  site  D3  replaced  after  lightning. 

HF  supply  faulty 

Lightning  put  the  array  out  of  operation  and 
it  remained  so  throughout  the  month. 

Daily  routine  check  of  NORSAR.  NORESS  and 
ARCESS  data. 

18-20  July 

4  July 

5  July 

17  July 

24  July 

July 

NORSAR: 

06C 

Replaced  SLEM  power  supply. 

1  Aug 

06C 

Cable  splicing  carried  out  SP03. 

17  Aug 

02C 

Located  faulty  cable  SP05. 

31  Aug 

NDPC: 

Daily  check  of  NORSAR.  NORESS  and  ARCESS 
data. 

Aug 

Calibration  every  week  of  NORSAR  LP  instru¬ 
ments.  Measurements  of  Mass  Pos.  and  Free 
Period,  adjusted  when  outside  tolerances. 

Aug 

ARCESS: 

Installed  broadband  instrument  KS36000  which 
has  been  in  USA  for  repair.  In  addition,  a 
new  NORSAT  B  antenna  erected. 

1-7  Sep 

NORESS: 

Electronics  on  remote  sites  A1.A3.B3.B5  repaired. 

12-15  Sep 

NORSAR: 

and  cable  D6  spliced. 

&  20  Sep 

01 A 

Repaired  damaged  cable  SP02. 

18-19  Sep 

01B 

LP  instruments  adjusted. 

22  Sep 

02C 

Repaired  JB  box  SP05. 

11  Sep 

03  C 

SP/LP  adjustments  carried  out. 

25  Sep 

04  C 

Test  Gen.  and  Protection  Card  on  SP04  replaced. 

21  Sep 

02B 

(Telemetry)  Station  2  was  repaired. 

26,27  Sep 

06C 

Broadband  filter  removed  from  ch.  05  and 
experiment  terminated. 

30  Sep 

NDPC: 

Daily  check  of  NORSAR.  NORESS  and  ARCESS 
data.  In  spite  of  a  failing  Modcomp  we  were 
able  to  calibrate  and  check  the  NORSAR  array. 

Sep 

Table  5.1  (cont.) 
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5.2  Array  status 

Broadband  filter  experiment  06C  SP  ch  05  was  terminated  by  the  end  of 
September  1989. 

As  of  30  Sep  1989  the  following  NORSAR  channels  deviated  from  toler¬ 
ances. 


01A  01  8  Hz  filter 
02  8  Hz  filter 
04  30  dB  attenuator 


O.A.  Hansen 
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7  Summary  of  Technical  Reports  /  Papers  Pub¬ 
lished 

7.1  A  performance  test  of  the  generalized  beamforming  method 
applied  to  a  data  base  of  regional  events 

Introduction 

In  the  real  time  processing  of  seismic  network  data,  one  of  the  most  important 
aspects  is  to  properly  associate  phases  detected  at  individual  stations.  For  a 
network  of  regional  arrays,  the  current  state-of-the-art  is  represented  by  the 
IAS  system,  where  the  Assess  procedure  is  used  for  phase  association  and 
location  estimation  (Bache.  1987). 

The  generalized  beamforming  method  introduced  by  Ringdal  and  Kvasrna 
(1989)  is  a  different  approach  that  has  shown  considerable  promise  as  an  ef¬ 
ficient  way  to  obtain  phase  grouping  and  initial  epicenter  estimates.  In  this 
paper,  the  method  is  applied  to  a  data  base  of  77  regional  and  local  events  and 
compared  to  IAS  results.  It  is  found  that  the  generalized  beamforming  is  able 
to  match  the  good  results  bv  Assess  as  far  as  phase  grouping  is  concerned, 
and  that  the  epicenter  determinations  are  also  quite  accurate. 

It  is  concluded  from  this  study  that  the  new  method  could  be  useful  as  a 
supplement  to  the  expert  systems  techniques  in  future  extensions  of  the  IAS 
system.  While  some  further  testing  would  be  required,  it  would  appear  that 
the  method  could  be  particularly  valuable  as  a  pre-processor  in  this  regard.  In 
this  way,  the  results  from  the  generalized  beamforming  could  be  taken  as  input 
to  a  procedure  where  the  rule-based  and  script-based  algorithms  implemented 
in  the  IAS  could  be  used  to  further  refine  the  solutions. 


Data  analysis 

The  data  set  of  event  recordings  from  N’ORESS  and  ARCESS  were  first  pro¬ 
cessed  by  the  IAS  detector  program,  called  Sigpro,  to  generate  detection  lists 
for  each  array.  The  detection  list  was  then  processed  by  the  IAS  rule-based 
expert  system,  called  Assess,  to  define  event  groups  and  event  locations. 

Subsequently  the  automatic  Assess  solutions  were  reviewed  by  an  analyst, 
and  if  necessary,  phases  were  retimed  or  new  phases  were  picked  to  get  more 
precise  locations.  When  false  events  occurred,  these  were  invalidated  by  the 
analyst. 

As  we  consider  the  analyst-reviewed  event  definitions  and  locations  to  be 
correct,  we  will  use  those  as  a  basis  for  evaluating  the  performance  of  the  gen- 
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eraiized  beamforming  method.  For  reference,  we  will  also  present  the  results 
produced  bv  Assess. 

1)  Identification  of  event  groups 


After  analyst  review  of  the  Assess  solutions.  76  events  were  observed. 
Phase  detections  belonging  to  a  common  event,  will  be  called  an  event  group. 
72  of  the  76  event  groups  were  validated  or  relocated  after  analyst  review.  The 
remaining  4  events  can  be  divided  into  the  following  classes: 

a)  One  multiple  event  that  was  observed  as  a  separate  event  group, 
but  with  wrong  phase  id’s.  This  event  was  not  relocated 

by  the  analyst. 

b)  One  triple  event  that  was  observed  as  one  separate  event  group, 
but  with  wrong  phase  id's.  These  events  were  not  relocated 

by  the  analyst. 

c)  One  multiple  event  that  was  observed  by  analyst,  but  was  not  given 
any  separate  event  group  by  Assess. 

Both  the  generalized  beamforming  method  and  Assess  in  a  few  cases  gave 
automatic  locations  that  deviated  considerably  from  the  locations  found  after 
analyst  review.  But  as  long  as  a  real  event  candidate  is  identified,  the  analyst 
can  later  make  the  necessary  adjustments  of  phase  id’s  and  arrival  times.  We 
will  therefore  consider  the  event  in  la)  as  identified  by  Assess,  one  of  the  two 
events  in  lb)  as  identified,  whereas  the  second  event  in  lb)  and  the  event  in 
lc)  is  missed. 

The  generalized  beamforming  method  found  71  of  the  72  events  that  were 
validated  or  relocated  after  analyst  replacement.  The  events  in  la)  and  lb) 
were  missed,  whereas  the  event  in  lc)  was  identified.  One  new  event  was 
identified,  so  we  set  the  number  of  events  observed  by  analyst  to  77. 

If  we  discard  the  multiple  events  in  1  a )- 1  c )  from  the  statistics,  the  gener¬ 
alized  beamforming  method  and  Assess  remain  with  one  missing  event  each, 
both  of  these  being  small  interfering  events  observed  on  one  array  mly. 


Analyst 

Generalized 

Assess 

Events  identified 

77 

73 

74 

Excluding  multiples 

73 

72 

72 
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2)  False  events 


An  important  factor  in  real-time  operation  of  seismic  networks  is  the  num¬ 
ber  of  false  event  groups  created  bv  the  automatic  system.  If  an  analyst  is 
going  to  review  the  automatic  output,  he  will  want  the  number  of  false  events 
to  be  at  a  reasonable  level.  By  split  two-arrav  events  we  mean  that  a  single 
event  was  detected  on  two  arrays,  but  the  automatic  algorithm  created  an  ex¬ 
tra  event  that  was  false.  By  split  one-array  event  we  mean  that  the  event  was 
detected  on  one  array  only,  but  the  automatic  algorithm  created  an  extra  event 
that  was  false.  Implausible  associations  were  decided  on  the  basis  of  detection 
list  information  like  amplitudes  and  frequency  content  of  the  different  phases, 
knowledge  of  propagation  characteristics  and/or  by  inspecting  the  actual  data. 


Generalized 

Assess 

Split  two-array  events 

1 

» 

Split  one-array  events 

0 

4 

Implausible  associations 

7 

4 

Number  of  false  events 

8 

15 

3)  Phase  identification 


The  generalized  beamforming  method  takes  advantage  of  the  travel-time 
pattern  of  the  detected  phases  when  event  groups  and  phase  identifications 
are  defined.  For  events  with  detections  on  one  array  only,  we  know  that  these 
patterns  in  many  cases  do  not  contain  sufficient  information  for  defining  the 
individual  phase  identifications.  In  those  situations  a  rule-based  approach  has 
to  be  invoked. 

For  39  events  in  t he  sample  data  base  with  detections  on  two  arrays,  we 
have  made  some  statistics  for  both  the  generalized  beamforming  method  and 
Assess  that  reflect  the  ability  to  assign  correct  phase  id's  to  the  different 
detections.  The  reference  will  be  the  phase  id’s  selected  by  analyst. 

The  columns  of  the  following  tables  need  some  explanation. 


correct  -  gives  the  number  ot  correctly  identified  detections. 

•■arly  gives  the  number  of  cases  where  a  detection  preceding  the 
correct  one  is  selected. 

late  -  gives  the  number  of  cases  where  a  detection  following  the 

correct  one  is  selected. 
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missed  -  gives  the  number  of  cases  were  a  detection  is  identified  by  the  analyst, 
but  missed  by  the  automatic  method. 

#  phases  -  means  number  of  detections  identified  as  defining  phases  by  the  analyst. 
#  phases=;correct+early+late+ missed 
The  number  of  events  analyzed  is  39.  as  indicated  earlier, 
false  -  gives  number  of  cases  where  the  automatic  method  picked  a  defining 
phase  that  was  not  hlentifiecLby  the  analyst. 


For  both  the  generalized  beamforming  method  and  Assess,  we  give  tables 
for  the  phase  identification  performance  on  Noress  and  Arcess  detections.  By 
adding  the  two  tables,  we  also  indicate  the  overall  performance. 

Generalized  beamforming 

NORESS 

I  correct  early  late  missed  #  phases  false 


rs-pn 
..rs-pg 
nrs-sn 
nrs-lg 


sum 


31 

0 

23 

24 


78 


0  5 

0  0 

0  6 

1  4 


1 


15 


37 

0 

29 

30 


96 


ARCESS 


correct 

early 

late 

missed 

#  phases 

false 

arc-pn 

34 

0 

4 

1 

39 

0 

arC'Pg 

3 

0 

0 

0 

3 

11 

arc-sn 

15 

1 

12 

2 

30 

3 

arc-lg 

13 

2 

10 

1 

26 

3 

sum 

65 

3 

26 

4 

98 

17 

NORESS  and  ARCESS 


correct 

early 

late 

missed 

#  phases 

false 

sum-pn  j 

65 

0 

9 

2 

76 

0 

sum-pg 

3 

0 

0 

0 

3 

12 

sum-sn 

38 

1 

18 

2 

59 

3 

sum-lg 

37 

3 

14 

2 

56 

5 

Assess 


NORESS 


correct 

early 

late 

missed 

#  phases 

false 

nrs-pn 

0 

1 

6 

37 

0 

nrs-pg 

0 

0 

0 

0 

0 

nrs-sn 

0 

0 

-  8 

29 

1 

nrs-lg 

1 

2 

30 

2 

sum 

71 

1 

3 

21 

96 

3 

ARCESS 


correct 

early 

late 

missed 

#  phases 

false 

arc-pn 

34 

0 

3 

•> 

39 

0 

arc-pg 

1 

0 

0 

2 

3 

0 

arc-sn 

12 

1 

3 

14 

30 

1 

arc-lg 

14 

2 

4 

6 

26 

4 

sum 

61 

3 

10 

24 

98 

5 

NORESS  and  ARCESS 


correct 

early 

late 

missed  ! 

1  #  phases 

false 

sum-pn 

64 

0 

4 

8 

76 

0 

sum-pg 

1 

0 

0 

2 

3 

0 

sum-sn 

33 

1 

3 

22 

59 

2 

sum-lg 

34 

3 

6 

13 

56 

6 

sum 

132 

4 

13 

45 

194 

8 

Within  the  data  base  there  were  9  one-array  events  with  three  or  more 
defining  phases.  That  all  of  these  events  were  observed  on  ARCESS.  reflect 
the  fact  that  very  few  events  within  the  data  base  had  epicenters  close  to 
NORESS.  In  cases  with  one-array  events  defined  by  three  or  more  phases, 
there  is  a  possibility  that  the  generalized  beamforming  method  can  utilize  the 
travel-time  pattern  to  come  up  with  correct  phase  id's. 

The  following  two  tables  indicate  the  performance  of  the  two  methods  on 
this  data  set. 
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Generalized  beamforming 


correct 

early 

late 

missed 

#  phases 

false 

arc-pn 

6 

0 

3 

0 

9 

0 

arc-pg 

2 

1 

3 

5 

arc-sn 

7 

1 

8 

i 

arc-lg 

5 

1 

6 

i 

sum 

20 

0 

5 

1 

26 

7 

Assess 


correct 

early 

late 

missed 

#  phases 

false 

arc-pn 

8 

1 

9 

0 

arc-pg 

1 

2 

3 

0 

arc-sn 

5 

3 

8 

0 

arc-lg 

i 

1 

2 

2 

6 

1 

sum 

15 

2 

2 

7 

26 

1 

4)  Initial  location  estimates 


Although  the  generalized  beamforming  method  gave  reasonable  phase  id’s 
for  the  9  one-array  events  considered  above,  we  can  not  give  very  much  con¬ 
fidence  to  the  initial  event  locations.  Within  the  predefined  azimuth  accep¬ 
tance  limits,  all  points  with  the  same  radius  from  the  observation  point  will  be 
equally  good  location  candidates.  As  we  have  not  run  any  location  program 
on  the  phases  identified  by  the  generalized  beamforming  method,  the  location 
estimates  for  these  events  will  not  be  very  accurate.  But  for  completeness,  we 
present  the  results. 

For  the  epicenter  estimates  obtained  by  the  generalized  beamforming  method 
and  Assess,  we  calculated  the  deviations  from  the  epicenters  computed  after 
analyst  review.  These  deviations  were  grouped  into  different  intervals,  as 
shown  below. 

Deviation  from  analyst  reviewed  epicenter  solutions  for  all  one-array  events 
with  three  or  more  defining  phases. 
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Deviation 
in  km 

Number  of  events 
Generalized 

within  distance  interval 
Assess 

d  <  1 0 

0 

3 

10<d<25 

0 

2 

25<d<50 

4 

1 

50<d<  100 

2 

2 

100<d<200 

3 

1 

200<d 

0 

0 

events 

9 

9 

In  the  following  table  we  consider  the  39  events  that  were  detected  on  two 
arrays. 


Deviation  from  analyst-reviewed  epicenter  solutions  for  all 

two-array  events 


Deviation 
in  km 

Number  of  events 
Generalized 

within  distance  interval 
Assess 

d<  10 

3 

8 

10<d<25 

8 

5 

25<d<50 

15 

12 

50<d<100 

7 

7 

100<d<200 

4 

2 

200<d 

2 

5 

events 

39 

39 

Discussion 

We  have  in  the  preceding  section  shown  that  when  processing  data  from 
two  arrays,  the  generalized  beamforming  method  identified  73  out  of  77  event 
groups.  When  discarding  multiple  events  from  the  same  location,  only  one 
small  interfering  one-array  event  was  missed. 

Some  false  event  groups  were  created,  and  this  requires  some  attention.  In 
the  cases  were  two  phases  from  two  arrays  are  erroneously  associated  together, 
a  check  on  the  amplitude  (magnitude)  of  the  two  phases  would  remove  some 
of  the  event  candidates. 

One  of  the  main  problems  in  automatic  identification  of  seismic  phases,  is 
the  relatively  large  number  of  detections  in  the  coda  of  the  seismic  signals.  If 
the  travel-time  functions  of  the  different  phases  do  not  fit  with  what  we  acually 
observe,  there  is  a  possibility  that  the  generalized  beamforming  method  will 
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assign  a  phase  id  to  one  of  the  coda  detections.  For  the  39  two-array  events, 
we  find  that  this  happened  in  -41  cases. 


As  an  additional  test,  we  wanted  to  check  if  we  could  improve  the  per¬ 
formance  by  introducing  some  simple  rules  to  the  procedure.  They  were  as 
follows: 


a)  If  a  phase  is  assigned  the  id  Pn.  and  the  event  group  contains  a 
preceding  Pn  candidate  (within  some  given  time  interval),  the 
preceding  phase  is  given  the  id  Pn. 

b)  If  a  phase  is  assigned  the  id  Sn.  and  the  event  group  contains  a 
preceding  Sn  candidate  (within  some  given  time  interval),  the 
preceding  phase  is  given  the  id  Sn. 

c)  If  a  phase  is  assigned  the  id  Lg,  and  the  event  group  contains  an 
Lg  candidate  with  a  larger  STA  value  (within  some  given  time 
interval),  that  phase  is  given  the  id  Lg. 


The  time  intervals  were  for  all  3  phases  set  to  10  seconds,  and  we  did  of 
course  also  check  how  these  rules  affected  the  previously  correct  phase  assign¬ 
ments.  The  results  are  shown  in  the  tables  below. 


Generalized  beamforming 


NORESS 


correct 

early 

late 

missed 

#  phases 

false 

nrs-pn 

36 

0 

0 

1 

37 

0 

nrs-pg 

0 

0 

0 

0 

0 

1 

nrs-sn 

27 

0 

2 

0 

29 

0 

nrs-lg 

26 

0 

3 

1 

30 

2 

sum 

89 

0 

5 

2 

96 

3 

ARCESS 


correct 

early 

late 

missed 

#  phases 

false 

arc-pn 

36 

0 

2 

1 

39 

0 

arc-pg 

3 

0 

0 

0 

3 

11 

arc-sn 

23 

1 

4 

2 

30 

3 

arc-lg 

14 

2 

9 

1 

26 

3 

sum 

76 

3 

15 

4 

98 

17 
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NORESS  and  ARCESS 


correct 

early 

late 

missed 

#  phases 

false 

sum-pn 

72 

0 

2 

2 

0 

sum-pg 

3 

0 

0 

0 

3 

12 

sum-sn 

50 

1 

6 

'2 

59 

3 

sum-lg 

40 

2 

12 

2 

56 

5 

sum 

165 

3 

20 

6 

194 

20 

What  we  found  was  that  these  simple  rules  improved  the  performance  on 
Pn  and  Sn  considerably,  but  the  Lg-rule  did  not  provide  almost  any  improve¬ 
ment.  Another  approach  that  could  possibly  be  used  to  better  identify  Lg  is 
to  compute  an  event  location  based  on  the  identified  phases  and  thereafter 
reselect  Lg  from  the  detection  that  best  fitted  the  predicted  Lg  arrival  time. 

But  the  fundamental  problem  is  that  for  many  regions,  we  do  not  have 
the  correct  tiavel-time  functions,  so  to  further  improve  the  performance  of  the 
generalized  beamforming  method,  we  need  to  introduce  regionalized  travel¬ 
time  functions  that  reflect  the  automatic  picks. 

For  the  9  one-array  events  with  three  or  more  defining  phases.  5  out  of  26 
phase  id’s  were  assigned  to  coda  detections.  By  introducing  the  same  rules  as 
outlined  for  the  two-array  events,  this  number  was  reduced  to  2.  as  shown  in 
the  table  below. 


Generalized  beamforming 


correct 

early 

late 

missed 

#  phases 

false 

arc-pn 

9 

0 

9 

arc-pg 

2 

0 

1 

3 

5 

arc-sn 

8 

0 

8 

1 

arc-lg 

4 

0 

1 

1 

6 

1 

sum 

23 

2 

1 

26 

1 

For  the  two-array  data  set.  the  analyst  only  identified  3  Pg  phases,  all 
on  ARCESS.  The  generalized  beamforming  method  correctly  identified  these 
phases,  but  in  addition.  12  Pn  coda  detections  were  assigned  the  id  Pg.  Some 
of  these  false  Pg  identifications  may  be  removed  by  giving  this  phase  some 
stronger  constraints  on  apparent  velocity  and/or  polarization  attributes. 

The  generalized  beamforming  produced  initial  locations  that  were  consis¬ 
tent  with  the  relatively  good  performance  on  phase  identification.  The  natural 
next  step  would  be  to  run  a  location  program  on  the  identified  phases,  and 
then  compare  to  the  location  obtained  after  analyst  review. 
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In  this  study  we  have  presented  results  on  the  performance  of  the  general¬ 
ized  beamforming  method  on  a  network  of  two  arrays.  The  next  step  will  be  to 
get  detection  data  (real  or  synthetic)  for  a  larger  network  of  arrays  and  single 
stations,  and  test  the  method  on  this  data  set.  Since  the  computer  programs 
for  this  method  is  almost  fully  parameterized,  only  minor  modifications  are 
needed  for  such  an  experiment. 
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7.2  NORSAR  P-wave  detection  and  yield  estimation  of  se¬ 
lected  Semipalatinsk  explosions 

In  the  most  recent  NORSAR  Semiannual  Technical  Summary  (NORSAR  Sci¬ 
entific  Report  No.  2-88/89)  results  are  presented  showing  that  the  Lg  phase 
provides  very  stable  measurements  of  magnitudes  of  Semipalatinsk  under¬ 
ground  nuclear  explosions,  thereby  showing-considerable  promise  with  regard 
to  explosion  yield  estimation  (Ringdal  and  Marshall,  1989).  A  possibility  to 
compare  Lg  and  P  magnitudes  to  published  yields  for  Semipalatinsk  explosions 
has  now  emerged  with  the  recent  publication  by  Soviet  scientists  quoting  yield 
estimates  for  a  number  of  such  explosions  (Bocharov  et  al.  1989.  Vergino. 
1989).  This  paper  presents  preliminary  results  from  a  study  comparing  NOR¬ 
SAR  measurements  to  these  yield  data. 

Yield  estimation  based  on  NORSAR  Lg  data 

We  have  made  an  effort  to  conduct  detailed  analysis  of  available  NORSAR 
Lg  data  for  the  largest  explosions  in  the  data  set  given  by  Bocharov  et  al 
(1989)  (Table  7.2.1).  The  results  are  as  follows: 

The  Shagan  River  event  of  11/02/72:  This  event  was  included  in  the  data 
base  of  Ringdal  and  Marshall  (1989),  with  an  M(Lg)  value  of  6.118. 

The  Shagan  River  event  of  12/10/72:  This  event  was  likewise  in  the  above 
data  base,  with  M(Lg)  =  6.116.  However,  a  factor  not  compensated  for  in 
this  estimate  was  the  occurrence  of  a  smaller  (Mi  =  5.6)  explosion  at  Dege- 
len  Mountains  with  origin  time  a  few  seconds  before  the  Shagan  River  event 
(Fig.  7.2.1).  The  quoted  M(Lg)  value  includes  a  contribution  from  Lg  waves 
from  this  smaller  event,  and  is  therefore  slightly  high.  Making  the  reasonable 
assumption  that  the  Lg  signals  from  the  two  events  add  up  incoherently,  i.e., 
that  the  signal  power  adds  up,  it  is  possible  to  compensate  for  this  interference 
in  the  same  way  as  the  noise  compensation  is  performed.  To  do  this,  we  need 
an  estimate  of  M(Lg)  for  the  smaller  event,  and  a  range  of  5. 5-5. 7  seems  rea¬ 
sonable.  The  corresponding  bias  values  are  found  to  range  from  0.013  to  0.030 
units,  with  a  mean  of  0.021.  Thus  a  revised  M(Lg)  estimate  of  6.116-0.021  = 
6.095  is  obtained,  with  an  uncertainty  of  about  ±  0.01  units  resulting  from 
this  compensation  procedure. 

The  Shagan  River  event  of  11/30/69:  We  have  been  able  to  recover  initial 
NORSAR  data  from  this  explosion.  At  the  time,  only  13  NORSAR  SP  chan¬ 
nels  (out  of  132)  were  in  operation,  and  5  of  these  were  among  the  42  channels 
used  in  the  previous  analysis.  Based  on  these  5  channels,  we  have  estimated 
M(Lg)  =  6.043  ±  0.021  for  this  event. 

The  Degelen  Mountains  event  of  4/25/71:  We  have  processed  this  event,  using 
exactly  the  same  procedure  as  for  the  Shagan  River  explosions,  with  the  same 
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time  windows,  filter  setting  and  channel  corrections.  The  resulting  estimate 
is  M(Lg)  =  5.862  ±  0.013,  based  on  42  NORSAR  channels. 

Fig.  7.2.2  shows  a  plot  of  M(Lg)  versus  published  yields  (Bocharov  et  al, 
1989)  for  these  four  explosions.  We  note  that  the  correspondence  is  excellent, 
although  there  are  of  course  far  too  few  data  points  to  draw  any  firm  con¬ 
clusion.  Nevertheless,  it  is  interesting  tojiote  that  the  Degelen  event  appears 
to  give  consistent  observations  with  the  Shagan  River  explosions,  which  also 
would  be  expected  from  the  earlier  observed  stability  of  the  M(Lg)  estimates. 

For  comparison,  Fig.  7.2.3  shows  a  plot  of  world-wide  mj  versus  these 
yield  values  for  the  same  four  events.  The  scatter  is  far  greater  than  for  the 
Lg  data.  Looking  at  the  differences  m^-log  Y,  it  is  interesting  to  note  that  the 
two  explosions  in  the  SW  and  NE  Shagan  River  area  show  a  relative  bias  of 
0.13  units,  which  is  very  close  to  the  m(,-M(Lg)  anomaly  of  0.15  magnitude 
units  found  by  Ringdal  and  Marshall  (1989)  between  these  two  subregions. 

Detection  of  low-yield  explosions 

The  cited  Soviet  publication  also  lists  a  number  of  low-yield  nuclear  explo¬ 
sions  conducted  at  Semipalatinsk.  We  have  analyzed  available  NORSAR  data 
for  these  explosions,  for  the  purpose  of  assessing  the  detection  thresholds  of 
the  NORSAR/NORESS  system  in  terms  of  explosion  yields.  Fig.  7.2.4  shows 
NORSAR  single  instrument  P-wave  recordings  of  the  three  smallest  explosions 
with  NORSAR  data  available:  2  Sep  72  (2  kt),  28  Mar  72  (6  kt)  and  16  Aug  72 
(8  kt).  We  have  chosen  to  display  the  instrument  06C02,  which  is  co-located 
with  the  center  seismometer  of  the  NORESS  array  (Note  that  NORESS  data 
have  only  been  available  since  1984).  The  traces  have  been  filtered  in  the  band 
2. 0-4.0  Hz,  and  the  signal-to-noise  ratios  (STA/LTA)  are  given  for  each  trace. 
Note  that  an  STA/LTA  treshold  of  4.0  is  generally  sufficient  to  confidently 
declare  a  signal  detection. 

It  is  clear  from  this  figure  that  the  NORSAR  detection  threshold,  at  the 
single  seismometer  level,  is  well  below  the  yields  of  the  three  events  shown. 
By  a  reasonable  extrapolation,  it  can  be  inferred  that,  given  similar  coupling 
conditions  and  noise  levels  as  in  these  three  examples,  Semipalatinsk  explosions 
of  yields  as  low  as  0.5  kt  would  be  expected  to  produce  detectable  signals  at 
NORSAR,  even  at  the  single  sensor  level. 

While  the  available  data  are  not  sufficient  to  draw  any  firm  conclusions 
regarding  the  NORSAR/NORESS  array  thresholds,  some  preliminary  obser¬ 
vations  may  nevertheless  be  made.  Thus,  it  has  been  established  in  earlier 
studies  that  the  SNR  gains  of  NORESS  in  the  filter  band  shown  are  at  least 
14  dB  (Kvaerna,  1989),  which  corresponds  to  a  factor  of  5  in  signal  ampli¬ 
tude.  This  indicates  a  detection  threshold  of  0.1  kt  or  lower  for  the  NORESS 
array  for  fully  coupled  Semipalatinsk  explosions,  assuming  normal  noise  con- 
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ditions.  Such  an  assertion  is  also  consistent  with  earlier  studies  indicating  that 
NORESS  has  a  detection  threshold,  in  terms  of  mt,  of  about  3.0  or  lower  for 
the  Semipalatinsk  region. 

Conclusions 

Studies  of  NORSAR  recordings  of  Semipalatinsk  explosions  for  which  So¬ 
viet  yield  estimates  are  available  confirm  that  NORSAR  Lg  magnitudes  have 
a  potential  for  providing  very  accurate  yield  estimates  in  the  magnitude  range 
where  there  is  sufficient  signal-to- noise  ratio.  Notably,  the  P-Lg  bias  earlier  ob¬ 
served  between  the  NE  and  SW  parts  of  the  Shagan  River  test  site  is  reflected 
in  a  similar  bias  between  m;,  and  log  (Yield). 

Furthermore,  the  excellent  signal  detectability  of  the  NORSAR/NORESS 
system  has  been  confirmed.  The  detection  threshold  for  fully  coupled  Semi¬ 
palatinsk  explosions,  assuming  normal  noise  conditions,  appears  to  be  below  1 
kt,  even  at  the  single  sensor  level.  Additional  significant  improvement  would 
be  obtainable  through  array  processing.  Clearly,  this  detection  level  will  not 
be  achieved  in  cases  where  the  noise  level  is  abnormally  high  (e.g.,  in  the  coda 
of  a  large  earthquake)  or  if  coupling  conditions  are  not  optimal  (e.g.,  in  case 
of  full  or  partial  cavity  decoupling).  It  must  also  be  noted  that  the  event  iden¬ 
tification  threshold  is  necessarily  higher  than  the  signal  detection  threshold, 
and  that  there  are  currently  no  seismic  methods  available  to  distinguish  reli¬ 
ably  between  small  underground  nuclear  explosions  and  chemical  explosions 
of  similar  yields. 


F.  Ringdal 


References 


Bocharov,  V.S.,  S.A.  Zelentsov  and  V.I.  Mikhailov  (1989):  Characteristics 
of  96  underground  nuclear  explosions  at  the  Semipalatinsk  test  facility, 
Atomic  Energy,  67,  3,  210-214  (in  Russian). 

Kvaerna,  T.  (1989):  On  exploitation  of  small-aperture  NORESS  type  arrays 
for  enhanced  P-wave  detectability.  Bull.  Seism.  Soc.  Am..  79,  888-900. 

Ringdal,  F.,  and  P.D.  Marshall  (1989):  Yield  determination  of  Soviet  un¬ 
derground  nuclear  explosions  at  the  Shagan  River  Test  Site.  NORSAR 
Semiannual  Technical  Summary,  1  Oct  1988  -  31  Mar  1989 ,  NORSAR 
Sci.  Rep.  2-88/89. 

Vergino,  E.S.  ( 1989):  Soviet  Test  Yields,  EOS,  Nov  "  1989,  1511-1525. 


56 


Date 

Region 

mb 

NORSAR 

M(Lg) 

Yield 

kt 

11/30/69 

Shagan  (TZ) 

6.02 

6.043 

125 

04/25/71 

Degelen 

5.94 

5.862 

90 

11/02/72 

Shagan  (SW) 

6.16 

6.118 

165 

12/10/72 

Shagan  (NE) 

5.96 

6.095*) 

140 

*)  Adjusted  for  interfering  explosion 


Table  7.2.1  Parameters  of  four  explosions  analyzed  in  this  study.  Yield  esti 
mates  are  from  Bocharov  et  al  (1989). 
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Shagan  River  12/10/72 
M(Lg)~6.116  (uncorrected) 


Fig.  7.2.1  NORSAR  recordings  (instrument  01A00)  of  the  double  explosion 
on  Dec  12,  1972.  Note  that  the  Lg  phase  is  a  superposition  of  signals  from 
both  sources,  although  the  contribution  from  the  smaller  explosion  (Degelen) 
is  relatively  modest.  By  estimating  a  confidence  interval  of  the  M(Lg)  for  the 
Degelen  event,  it  is  possible  to  compensate  for  the  effects  on  the  M(Lg)  for 
the  Shagan  River  event,  as  described  in  the  text. 
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MAGNITUDE  -  YIELD 

SEM I  PA  LA  TINSK  EXPLOSIONS 

S=  0.90  /=  4.136  D=  0.028  N=  4 


Log  Yield  (Bocharov  et  al) 


Fig.  7.2.2  Comparison  of  M(Lg)  and  log(yield),  of  the  four  explosions  in  Table 
7.2.1,  with  the  yield  values  taken  from  Bocharov  et  al  (1989).  The  slope  of  the 
straight  lines  has  been  restricted  to  0.9.  Note  the  excellent  correspondence. 
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MAGNITUDE  -  YIELD 

SEMI  PA  LATINS  K  EXPLOSIONS 
S=  C.90  /=  4.127  D=  0.069  N=  4 


Fig.  7.2.4  NORSAR  single  seismometer  recordings  (instrument  06C02)  of 
three  low-yield  Semipalatinsk  nuclear  explosions.  The  traces  have  been  fil¬ 
tered  in  the  band  2. 0-4.0  Hz.  and  the  signal-to-noise  ratio  is  shown  for  each 
trace.  Note  the  very  high  SN'R  observed  even  at  the  single  sensor  level.  Sig¬ 
nificant  additional  SNR  gain  would  be  obtained  through  array  processing  of 
the  NORSAR  data. 
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7.3  Analysis  of  data  from  the  China  Digital  Seismograph 
Network  (CDSN)  for  Soviet  nuclear  explosions 

This  paper  is  a  follow-up  to  earlier  work  (Ringdal  and  Marshall.  1989:  Hansen 
et  al.  1989)  aimed  at  evaluating  the  stability  of  seismic  Lg  magnitudes  for 
yield  estimation  purposes.  In  particular,  these  efforts  have  involved  analyz¬ 
ing  available  Lg  data  from  Soviet  nuclear  explosions  at  the  Shagan  River. 
Semipalatinsk  test  site,  and  conducting  comparative  analyses  of  Lg  and  P 
recordings  at  various  seismograph  stations. 

Hansen  et  al  ( 1989 )  analyzed  data  recorded  at  four  digital  stations  installed 
by  IRIS  in  the  Soviet  Union,  and  found  an  excellent  correspondence  between 
Lg  measurements  at  these  stations  and  the  NORSAR  M(Lg)  estimates  pub¬ 
lished  by  Ringdal  and  Marshall  ( 1989).  Furthermore,  they  noted  the  very  high 
Lg  signal-to-noise  ratio  observed  at  the  IRIS  stations,  in  particular  ARl’  and 
GAR.  and  concluded  that  reliable  Lg  measurements  at  these  stations  would  be 
possible  for  explosions  as  small  as  m*,  =  4.0.  assuming  normal  noise  conditions. 

In  this  paper,  we  extend  the  analysis  to  data  from  the  China  Digital  Seis¬ 
mograph  Network  (CDSN),  which  is  operated  by  the  USGS  in  cooperation 
with  the  State  Seismological  Bureau,  Beijing.  Two  of  the  CDSN  stations. 
VVMQ  in  Urumqi  and  HIA  in  Hailar.  have  particularly  good  Lg  propagation 
paths  from  Semipalatinsk,  and  we  have  based  our  analysis  on  data  from  these 
two  stations. 

Fig.  7.3.1  shows  the  locations  of  the  two  stations  in  relation  to  the  test 
site,  as  well  as  locations  of  the  NORSAR  and  the  IRIS  stations.  WMQ  has  an 
epicentral  distance  to  Shagan  River  of  960  km.  whereas  HIA  is  at  a  distance  of 
about  3000  km.  Both  stations  show  excellent  Lg  recordings  of  Semipalatinsk 
explosions,  as  illustrated  by  the  examples  in  Figs.  7.3.2  and  7.3.3. 

In  our  analysis  of  WMQ  and  HIA  Lg  recordings,  we  have  employed  the 
exact  same  procedure  as  described  for  IRIS  data  by  Hansen  ct  al  (1989). 
and  the  details  will  not  be  repeated  here.  Data  from  a  total  of  12  Shagan 
River  explosions,  dating  back  to  1987,  were  provided  to  us  for  this  analysis  by 
the  Center  for  Seismic  Studies.  Table  7.3.1  lists  these  events  along  with  the 
estimated  parameters. 

Fig.  7.3.1  shows  a  comparison  of  WMQ  and  NORSAR  log  RMS  (Lg)  esti¬ 
mates  for  these  12  events.  The  slope  of  the  plot  has  been  restricted  to  1.00. 
and  the  standard  deviation  of  the  differences  between  the  two  stations  is  only 
0.031  units.  This  is  essentially  the  same  scatter  found  earlier  by  Hansen  et  al 
( 1989)  when  comparing  data  from  NORSAR  and  the  Soviet  station  ARC.  and 
confirms  the  excellent  stability  of  the  RMS  Lg  estimates. 

Fig.  7.3.'>  shows  a  comparison  of  HIA  and  NORSAR  log  RMS  (Lg)  esti- 
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mates.  In  this  case,  the  slope  of  the  least  squares  linear  relationship  ( 1.48)  is 
significantly  different  from  unity,  and  we  note  that  a  similar  observation  was 
also  made  bv  Ringdal  and  Marshall  (1989)  when  comparing  NORSAR  and 
Grafenberg  Lg.  We  will  not  go  into  any.  detail  discussing  possible  underlying 
physical  reasons  for  this  variability  in  slopes.  For  our  purpose,  the  important 
point  is  to  note  that  the  scatter  of  the  relationship  is  still  very  small:  the 
standard  deviation  in  the  y-axis  direction  being  0.041  units.  The  "orthogo¬ 
nal”  standard  deviation  relative  to  the  straight  line  fit  is  0.023,  which  in  fact 
compares  very  closely  to  the  orthogonal  standard  deviation  of  0.024  which  can 
be  inferred  from  the  WMQ-NORSAR  data  shown  in  Fig.  7.3.4. 

In  Fig.  7.3.6  we  plot  the  II IA  versus  VV'MQ  log  RMS  (Lg)  values,  and  again 
observed  that  the  least-squares  slope  ( 1.36)  is  significantly  different  from  unity. 
Once  more,  the  scatter  is  very  small,  with  an  orthogonal  standard  deviation 
of  0.028  units. 

Fig.  7.3. 7  is  a  plot  comparing  WMQ  Lg  data  with  maximum  likelihood  ISC 
m*  estimates.  Compared  to  the  previous  figures,  this  plot  shows  a  somewhat 
greater  standard  deviation  of  0.060  measured  in  the  y-axis  direction.  This 
scatter  is  still  quite  small,  but  it  must  be  noted  that  only  one  event  from 
the  northeast  part  of  Shagan  is  in  the  data  base.  Thus,  we  cannot  assess 
whether  the  M(Lg)  versus  mj,  bias  earlier  found  for  this  subregion  (Ringdal 
and  Marshall,  1989)  is  also  present  when  measuring  Lg  at  WMQ. 

Finally,  Fig.  7.3.8  compares  the  signal-to-noise  ratios  (defined  as  Lg  signal 
to  pre-P  noise)  for  stations  at  various  distances,  using  5  large  explosions.  It 
is  noteworthy  that  WMQ  shows  the  best  SNR  for  all  the  events.  The  figure 
suggests  that  WMQ  would  be  able  to  give  Lg  measurements  for  events  two 
magnitude  units  smaller  than  the  NORSAR  threshold  of  approximately  5.5. 
I  nfortunately.  there  were  no  low  magnitude  events  in  our  data  base,  so  we 
have  not  been  able  to  confirm  this  hypothesis. 

In  conclusion,  our  studies  confirm  that  Lg  magnitude  estimates  of  Semi- 
palatinsk  explosions  are  remarkably  consistent  between  stations  widely  dis¬ 
tributed  in  epicentral  distance  and  azimuth.  It  thus  appears  that  a  single 
station  with  good  signal-to-noise  ratio  can  provide  M(Lg)  measurements  with 
an  accuracy  (one  standard  deviation)  of  about  0.03-0.04  magnitude  units. 
Thus,  the  Lg  phase  shows  considerable  promise  for  use  in  yield  determination, 
although  more  data  will  be  needed  before  the  accuracy  of  Lg-estimated  yields 
can  be  firmly  established. 


R.A.  Hansen 
F.  Ringdal 
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No. 

Date 

tnh 

NAO  Lg 

WMQ  Lg 

HIA  Lg 

1 

87171 

6.03 

3.012 

3.851 

2.189 

2 

87214 

5.83 

2.911 

3.693 

2.072 

3 

87319 

5.98 

3.014 

3.870 

2.298 

4 

87347 

6.06 

3.133 

3.907 

2.352 

5 

87361 

6.00 

3.086 

3.851 

2.339 

6 

88044 

5.97 

33)82 

3.911 

- 

7 

88094 

5.99 

3.103 

3.925 

2.307 

8 

88125 

6.09 

3.084 

3.958 

- 

9 

88258 

6.03 

3.014 

3.827 

2.224 

10 

88317 

5.20 

2.307 

3.104 

11 

88352 

o 

00 

lO 

2.846 

3.636 

1.947 

12 

89043 

5.90 

2.836 

3.619 

1.921 

Table  7.3.1  Magnitudes  (mi,)  and  log  RMS  Lg  values  at  NORSAR,  YVMQ 
and  HIA  for  12  explosions  analyzed  in  this  study. 
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Fig.  7.3.1  Map  indicating  the  locations  of  the  Shagan  River  lest  Site,  the 
IRIS  stations  in  the  USSR,  the  NORSAR  array  in  Norway  and  the  stations 
WMQ  and  HIA  in  China. 
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Fig.  7.3.2  Example  of  recordings  from  a  Soviet  nuclear  explosion  (3  April 
1988)  at  the  station  YVMQ.  For  each  of  the  three  components  we  show  the 
unfiltered  trace  (bottom),  the  filtered  trace  (0. 6-3.0  Hz)  and  the  1'20-second 
window  RMS  measure  (top)  as  a  function  of  time. 
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Fig.  7.3.3  Kxamplcof  wordings  from  the  Soviet  JVF  explosion  I  M  Sep  1!),SS) 
at  the  station  HIA.  The  three  traces  for  each  component  are  as  on  Fig.  7.3.2. 


Log  (RMS  Lg)  at  WMQ 


Z  Component  Lg  RMS  Comparison 


Log  (RMS  Lg)  at  NAO 


Fig.  7.3.4  Comparison  of  log  RMS  Lg  measurements  obtained  at  WMQ  and 
NORSAR.  The  standard  deviation  of  the  differences  is  0.034  in  the  y-direction 
and  0.024  orthogonal  to  the  line.  The  dotted  lines  correspond  to  plus  or  minus 
two  standard  deviations. 
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Z  Component  Lg  RMS  Comparison 


Log  (RMS  Lg)  at  WMQ 


Fig.  7.3.0  Comparison  of  log  RMS  Lg  measurements  at  HIA  and  WMQ.  The 
slope  of  the  line  is  1  and  the  standard  deviation  of  the  misfit  of  the  line 
to  the  data  is  0.047  in  the  y-direction  and  0.028  orthogonal  to  the  line.  The 
dotted  lines  correspond  to  plus  or  minus  two  standard  deviations. 


71 


Z  Component  of  Log  (RMS  Lg)  at  WMQ 


Fig.  7.3. 7  C  comparison  of  log  RMS  Lg  at  WMQ  to  world-wide  nif,  magnitude. 
Standard  deviation  is  0.060  units  in  the  y-direction  and  0.042  orthogonal  to  the 
line.  The  dotted  lines  correspond  to  plus  or  minus  two  standard  deviations. 


Signal/Noise  for  RMS  Lg 
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Fig.  7.3.8  Graph  showing  the  variation  of  the  signal-to-noise  ratios  (log  RMS 
minus  log  RMS  noise)  from  the  four  IRIS  stations,  the  NORSAR  array  and  the 
CDSN  stations  WMQ  and  HIA.  Epicentral  distance  to  the  test  site  is  plotted 
along  the  horizontal  axis. 
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7.4  A  new  data  acquisition  system  for  FINESA 

Hie  FINESA  array  was  installed  in  Finland  in  1985  and  comprised  in  its 
original  version  10  SPZ  seismometers.  Results  from  data  analysis  using  the 
FINESA  array  in  its  initial  configuration  are  given  in  Ringdal  et  cil  (1987). 
In  1988.  the  geometry  of  the  array  was  expanded  by  the  addition  of  five  <-le- 
'ents.  and  the  FINESA  array  geometry  currently  comprises  15  vertical  only 
seismometers  within  an  aperture  of  2  km.  as  shown  in  Fig.  7.4.1.  An  evalu- 
aion  of  the  performance  of  the  expanded  array  is  offered  by  Mykkc-ltveit  1 1  al 
(  1988). 

I  ntil  this  year  ( 1989).  t  he  only  way  of  acquiring  FINESA  data  has  been  by 
on-site  recording  on  magnetic  tapes,  and  all  analysis  of  FINESA  data  has  been 
strictly  off-line.  With  the  development  of  the  Intelligent  Array  processing  Sys¬ 
tem  (IAS).  it  has  become  an  objective,  however,  to  provide  the  capability  of 
transmitting  data  irom  FINESA  in  real  time  to  the  NORSAR  data  processing 
center.  At  the  same  time,  it  has  been  an  objective  to  enable  real  time  trans¬ 
mission  of  F  INESA  data  to  Helsinki,  for  the  independent  analysis  of  this  data 
by  the  staff  of  the  Institute  of  Seismology  at  the  University  of  Helsinki.  A  new 
data  acquisition  system,  meeting  the  above  objectives,  has  been  developed 
at  NORSAR  and  put  into  operation.  Short  descriptions  of  the  field  system, 
the  communications  interface  used,  tiie  data  transmission  lines  that  have  been 
established,  and  the  SUN-based  data  acquisition  and  processing  systems  in¬ 
stalled  at  NORSAR  and  Helsinki,  are  the  cubjects  of  this  contribution. 

The  field  system 

The  field  unit  is  a  VME-based  system  manufactured  by  Force  Comput¬ 
ers.  The  software  was  developed  in  Norway  by  the  company  Data  Respons 
A/S.  under  the  guidance  and  supervision  of  NORSAR  personnel.  A  detailed 
description  ot  this  field  system  e  provided  in  a  manual  prepared  by  Data  Re¬ 
spons  A/S.  entitled  NORSAR  Field.  System  Description  Manual.  This  manual 
is  available  upon  request. 

The  main  functions  of  the  field  unit  are: 


•  A/D  conversion  at  a  sampling  rate  of  It)  11/.  of  up  to  42  channels  (15 
channels  currently  in  use) 

•  Synchronized  sampling  using  a  Kinemetrics  clock  feeding  the  system 
with  (JMT  time 

•  Data  buffering,  reformatting  and  adding  of  header  information 

•  SDI.C  communication. 


74 


The  entire  FINESA  acquisition  system  is  shown  schematically  in  Fig.  7.4.2. 
The  units  located  at  the  FINESA  array  site  in  Sysma  are  shown  in  the  upper 
left  corner:  The  field  unit  (Force  Computer  with  A/D  cards),  Kincmetrirs 
clock,  an  RS232  to  V.35  converter,  and  a  modem. 

Communications  interface  and  data  transmission  lines 

NORSAR  has  recently  developed  a  communications  system  named  BI’SC- 
VME  (Buffered  5DLC  Communications  system  for  SUN  VME-based  systems). 
Our  major  requirements  for  a  new  communications  interface  were: 


•  Low  cost 

•  VME-based 

•  SDLC  communication 

•  No  hardware  developments 

•  Easy  and  low  cost  maintenance. 

The  BUSC-VME  system  meets  these  requirements  and  is  easily  interfaced 
with  a  SUN  VME-based  system.  The  BUSC-VME  system  is  based  on  com¬ 
mercially  available  cards  for  installation  directly  into  a  SUN  system. 

The  BUSC-VME  system  is  based  on  two  VME  cards  (CPU  and  memory) 
installed  in  the  SUN  VME  bus.  There  is  no  requirement  for  outside  units. 
The  SDLC  stream  is  plugged  directly  into  the  CPU  board.  The  CPLT  board  is 
MC68020-based  and  handles  two  high  speed  SDLC  links  (array  and  HF-data 
in  case  of  NORESS:  array  data  only  in  the  current  1  i.NESA  implementation). 
Data  received  by  the  CPU  card  are  stored  in  a  cyclic  buffer  on  a  memory 
board.  A  2  Mb  board  has  buffer  capacity  for  10  minutes  of  array  data.  The 
buffer  size  may  easily  be  changed.  The  transfer  speed  between  Sc  .'.  and  the 
BUSC-VME  system  is  only  limited  by  the  SUN  VME  bus  and  the  activity  on 
the  SUN. 

The  software  driver  is  very  simple  since  the  BUSC-VME  memory  is  mapped 
into  the  SUN  VME  bus  and  there  is  no  need  for  interrupt  routines  to  com¬ 
municate  with  the  memory  board.  The  driver  simply  reads  data  out  of  the 
buffer  like  a  program  fetches  data  from  the  regular  SUN  memory  space.  All 
handshaking  is  controlled  on  a  VME  hardware  level  and  has  no  influence  on 
the  SUN  itself. 

Two  spare  cards  are  the  only  requirements  for  backup  and  an  acceptable 
maintenance  scheme  for  the  communications  system. 


The  BUSC-VME  interface  has  been  installed  in  the  SUN  computers  at 
N'ORSAR  and  Helsinki  to  enable  reception  of  FINESA  data  at  these  locations. 

High  speed  (64  Kbit/s)  digital  land  lines  have  been  installed  between  the 
FINESA  array  site  and  Helsinki  and  between  Helsinki  and  N'ORSAR.  As 
shown  in  Fig.  7.4.2.  the  data  stream  is  split  in  Helsinki,  and  forwarded  to 
both  NORSAR  and  the  Seismological  Institute  of  the  University  of  Helsinki. 

Data  acquisition  and  processing 

Data  from  FINESA  are  received  continuously  in  real  time  both  at  N'OR¬ 
SAR  and  in  Helsinki,  as  shown  in  Fig.  7.4.2.  The  data  acquisition  and  process¬ 
ing  systems  are  based  on  SUN  3/2xx  computers.  Both  systems  are  equipped 
with  local  disk,  and  the  SUN  system  at  NORSAR  also  has  an  Exabyte  tape 
drive  for  data  archiving. 

The  data  received  are  written  on  to  a  cyclic  diskloop.  and  only  available 
disk  space  limits  the  size  of  this  loop.  For  the  NORSAR  installation,  the 
disk  drive  installed  permits  a  loop  containing  approximately  70  hours  of  data, 
which  matches  the  length  of  the  disk  loops  implemented  for  the  NORESS  and 
ARCESS  data  acquisition  systems. 

The  detection  processing  (DP)  and  event  processor  (EP)  program  packages 
described  by  Fyen  (1989)  have  been  installed  on  the  acquisition  computers 
both  at  NORSAR  and  in  Helsinki.  The  beam  deployment  chosen  was  designed 
from  the  NORESS/ARCESS  experience  and  adapted  to  the  FINESA  array 
geometry.  An  example  of  an  event  automatically  detected  and  located  (after 
an  association  of  a  P  and  an  Lg  phase)  by  the  DP/EP  packages  is  shown 
in  Figs.  7. 4. 3-7. 4. 5.  The  event  is  a  small  mining  explosion  in  Estonia.  The 
format  of  the  plots  and  the  information  contained  therein  are  as  explained  in 
Fyen  ( 1989). 

In  the  future,  data  from  FINESA  will  be  processed  jointly  with  NORESS 
and  ARCESS  data  (as  well  as  GERESS  data,  when  available)  in  the  IAS. 


R.  Paulsen 
J.  Fyen 
P.W.  Larsen 

S.  Mykkeltveit 
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Fig.  7.4.1  (  configuration  of  th<>  HNFSA  array  in  Sysnia.  Finland,  after  the 
expansion  of  the  array  geometry  in  IPSX. 
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Fig.  7.4.4  Expanded  display  of  the  P- phase  for  the  event  of  Fig.  7.-J.3.  \T06dpx 
is  the  detecting  beam  that  gives  the  best  SNR  for  this  phase.  This  is  an  inco¬ 
herent  beam  with  infinite  velocity,  but  plotted  here  as  a  coherent  beam,  using 
the  same  steering  parameters.  YI06  (shown  both  filtered  and  unfiltered)  is  the 
coherent  beam  using  ti«*»  same  configuration,  and  steered  according  to  the  f-k 
estimate  (here:  phase  velocity  6.50  km/s:  azimuth  165.90°). 


Fig.  7.4.5  Same  as  Fig.  7.4.4.  but  corresponding  to  the  Lg  phase. 


7.5  Establishment  of  a  mining  explosion  data  base 

On  a  daily  basis,  numerous  local  and  regional  events  are  detected  and  lo¬ 
cated  by  the  signal  analysis  procedures  applied  to  the  high-frequency  arrays 
NORESS  and  ARCESS.  The  majority  of  these  events  are  related  to  mining 
activity  in  the  area  monitored  by  these  arrays.  Specific  knowledge  pertaining 
to  mining  explosions  such  as  mine  locations’,  charge  sizes,  firing  practices,  etc., 
will  aid  in  the  regular  analysis  of  the  array  data.  In  addition,  such  information 
will  be  the  basis  for  adequately  addressing  a  range  of  interesting  and  important 
topics  within  seismoiogical  verification  research. 

With  the  purpose  of  obtaining  such  information,  a  survey  of  mining  activity 
in  the  areas  monitored  by  NORESS  and  ARCESS  was  initiated  during  the 
reporting  period.  The  work  consists  of  four  steps  including  gathering  of  back 
ground  material,  acquiring  information  from  the  mines,  defining  relational 
tables  and  selecting  time  series  from  our  data  archives. 

Collection  of  background  information 

The  authorization  of  mining  activity  in  Scandinavia  is  decentralized  to  lo¬ 
cal  pCice:  however,  there  is  no  official  regulation  requiring  recording  of  blasting 
parameters,  etc.  The  initial  work  towards  establishing  a  mining  explosion  data 
base  therefore  was  aimed  at  obtaining  an  overview  of  activities  from  branch 
organizations  and  companies  relating  to  mining  activity.  The  survey  of  ac¬ 
tivity  resulted  in  a  list  of  24  Swedish  and  43  Norwegian  companies  some  of 
which  may  operate  at  several  locations.  These  companies  represent  presum¬ 
ably  continuous  operation  of  mining  activity  underground  or  at  the  surface, 
or  simply  rock  quarries  for  the  purpose  of  producing  building  materials.  At 
present,  no  effort  has  beer  made  to  map  temporary  mining  activity  relating 
to  various  construction  projects  or  military  activity.  However,  data  regarding 
larger  naval  detonations  are  available  via  a  reporting  schedule  to  the  Univer¬ 
sity  of  Bergen.  A  list  of  explosions  from  the  largest  open  pit  mine  in  southern 
Norway  (Titania)  was  also  obtained  as  a  part  of  a  reporting  program  to  the 
University  of  Bergen. 

Submission  of  a  questionnaire  to  actual  mines  and  quarries 

Common  practice  regarding  filing  of  information  about  past  explosions 
varies  considerably  from  one  mine  to  another,  as  the  mining  companies  are 
not  required  to  provide  reports  on  their  explosions.  A  questionnaire  was  sub¬ 
mitted  to  the  67  companies  identified,  asking  for  a  number  of  details  about 
blasting  practice  and  explosions  detonated  since  1984.  Out  of  the  67  companies 
approached,  answers  were  obtained  from  16.  corresponding  to  an  answering 
per  cent  of  about  25.  The  quality  of  the  responses  varied  greatly,  apparently 
due  to  the  lack  of  details  in  their  files.  The  sparse  response  may  also  be  at- 
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tributed  to  the  fact  that  the  workload  on  the  part  of  the  companies  was  not 
associated  with  any  compensation  by  N'ORSAR.  A  map  showing  the  location 
of  mines  presently  contributing  to  the  explosion  data  base  is  given  in  Fig. 
7.5.1.  The  mines  from  which  data  are  acquired  include  four  of  the  main  open 
pit  mining  sites  in  Norway,  located  in  southern  Norway,  western  Norway,  mid 
Norway  and  northern  Norway,  respectively. 

Details  of  explosion  locations  at  the  mine  in  mid- Norway  (Rana  Gruber) 
are  shown  in  Fig  7.5.2,  and  for  the  mine  in  northern  Norway  (Sydvaranger) 
located  close  to  the  Russian  border  in  Fig.  7.5.3.  The  latter  mine  performs 
blasting  at  five  different  locations  within  an  area  of  about  10  km  across.  The 
distances  to  ARCESS  and  NORESS  from  these  four  mines  are  shown  in  Fig. 
7.5.4.  The  amounts  of  explosives  detonated  vary  from  less  than  10,000  kg 
to  more  than  200.000  kg,  thus  providing  a  good  range  of  charge  sizes  and 
distances  for  consideration  in  future  research. 

Establishing  mining  explosion  relations  according  to  CSS  data  base 
structure 

The  mines  responding  to  the  questionnaire  do  provide  a  basis  for  defin¬ 
ing  elements  of  a  mining  explosion  data  base  for  future  research  and  master 
event  selections.  Two  relations  are  proposed  in  conformance  with  the  2.8  CSS 
data  base  structure  for  the  mining  explosions.  The  relational  table  MINFO 
contains  information  specific  for  one  particular  mine  location,  but  common  to 
many  explosions  in  that  mine.  The  description  of  individual  explosions  is  con¬ 
tained  in  the  relation  MINEX.  Details  of  the  proposed  relations  are  provided 
in  Appendix  A. 

The  information  contained  in  the  returned  questionnaiies  has  been  edited 
to  a  computer  file  and  subsequently  processed  and  rewritten  to  a  file  in  confor¬ 
mance  with  CSS  external  file  format.  The  processing  includes  conversion  from 
local  time  to  EPOC-time  and  conversion  from  local  coordinates  provided  by 
the  Norwegian  Mapping  Authority  (NGO)  to  geodetic  latitude  and  longitude. 
Subroutines  needed  for  this  transformation  were  acquired  and  implemented  at 
NORSAR  as  part  of  this  project. 

Examples  of  the  external  files  of  the  relational  tables  MINFO  and  MINEX 
are  presented  in  Tables  7.5.1  and  7.5.2.  The  description  of  ripple  firing  is 
given  in  terms  of  one  tuple  for  each  individual  detonation  in  the  ripple.  For 
many  ripple-fired  explosions  reported  by  the  mines,  detailed  descriptions  of 
the  cap  delay  times  and  charge  distribution  were  not  provided.  In  cases  where 
the  charge  configuration  parameter  ‘chacon’  is  set  to  'obtainable',  details  of 
ripples  may  be  obtained  upon  specific  request  and  entered  into  the  tables  at 
a  later  time. 


Selection  of  master  recordings  of  mine  explosions 

The  work  towards  entering  typical  mine  explosion  recordings  at  our  a  rays 
into  master  files  is  at  a  preliminary  stage,  and  must  await  the  consideration 
of  a  large  number  of  recordings  supported  by  data  processing  results  in  terms 
of  energy  distribution  in  time,  frequency  and  space.  Crucial  issues  for  further 
consideration  ar  the  attenuation  of  the  seismic  waves  from  various  types  of 
mining  explosions  and  the  characteristics  of  ripple  firing  in  time-frequency 
estimates  (spectograms)  of  such  recordings;  see,  e.g.,  Baumgardt  and  Ziegler 
(1988),  Hedlin  et  al  (1989),  Richards  et  al  (1989).  Some  examples  of  mine 
recordings  are  shown  in  Figs.  7.5.5  to  7.5.7. 

Plans  for  further  work 

From  the  information  gathered  from  the  mining  industry  so  far.  we  are 
in  a  good  position  to  judge  which  of  the  mines  that  have  not  yet  returned 
the  questionnaire  are  of  interest  to  us.  We  will  make  an  effort  to  contact 
these  mines  again.  Many  of  the  mines  in  this  category  are  located  in  Sweden, 
and  we  expect  to  obtain  some  help  and  advice  from  Swedish  seismologists  in 
establishing  contact  with  the  mines  in  question.  After  this  second  round  of 
information  gathering,  we  hope  to  be  able  to  compile  a  fairly  comprehensive 
data  base  of  mining  explosions  in  Norway  and  Sweden. 

The  network  of  regional  arrays  providing  real  time  data  to  the  NORSAR 
Data  Processing  Center  is  currently  being  expanded.  Very  soon,  data  from 
the  FINESA  array  in  Finland  and  the  GERESS  array  in  the  Federal  Republic 
of  Germany  will  become  available.  It  will  then  be  necessary  to  expand  the 
survey  of  mining  activity  described  in  this  report  to  new  regions,  in  particular 
Finland,  the  western  USSR,  Poland,  Czechoslovakia  and  the  FRG.  Again,  we 
expect  to  be  able  to  obtain  help  and  information  from  colleagues  in  these 
countries. 


A.  Dahle 
A.  Alsaker 
S.  Mykkeltveit 
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ine  No. 

Name 

Lat,  deg 

Lon,  deg 

Mine  Type 

1 

Olivin 

62.04000 

5.52333 

open  pit 

3 

Rana  Gruber 

66.42063 

4.67853 

open  pit 

4 

Svdvaranger 

69.65189 

30.02533 

open  pit 

5 

Sulitjelma 

67.14769 

16.07007 

shaft 

6 

Fz.f.Bryggja 

61.93000 

5.45000 

shaft 

7 

Grong 

64.87000 

13.88000 

shaft 

8 

Bleikvassli 

65.93000 

13.88000 

shaft 

9 

Fosdalen 

64.07000 

11.20000 

shaft 

13 

Steens  Kalkverk 

60.74000 

11.02000 

open  pit 

14 

Fz.f.Sandvika 

59.90000 

10.50000 

shaft/pit 

15 

Titania 

58.30000 

6.40000 

open  pit 

16 

Fz.f.V'interbro 

59.75000 

10.77000 

open  pit 

101 

Vammala 

61.33000 

23.03000 

102 

Enonkoski 

62.04000 

28.77000 

103 

Taipalsaari 

61.18000 

28.04000 

104 

Pyhaesalmi 

63.66000 

26.05000 

105 

Hitura 

63.85000 

25.05000 

106 

Vihanti 

64.41000 

25.15000 

107 

Elijaervi 

65.78000 

24.70000 

108 

Saattopora 

67.78999 

24.43000 

109 

Faarby 

60.10000 

22.88000 

110 

Mustio 

60.1/000 

23.84000 

111 

Sipoo 

60.25000 

25.39000 

112 

Tytyri 

60.27000 

24.07000 

113 

Parainen 

60.30000 

22.29000 

114 

Ihalainen 

61.03000 

28.18000 

115 

Vampula 

61.05000 

22.64000 

116 

Otama 

61.81000 

21.74000 

117 

Ruokojaervi 

61.94000 

29.03000 

118 

Ankele 

62.07000 

27.41000 

119 

Ryytimaa 

63.15000 

24.02000 

120 

Kalkkimaa 

65.90000 

24.47000 

121 

Kemiaa 

60.14000 

22.59000 

122 

Horsmanaho 

62.82000 

29.25000 

123 

Lipasvaara 

63.02000 

29.23000 

124 

Siilinjaervi 

63.12000 

27.74000 

125 

Nilsiae 

63.16000 

27.99000 

126 

Lahnaslarnpi 

64.12000 

28.06000 

.Mine  No. 

Name 

Lat.  deg 

Lon,  deg 

Mine  Type 

201 

Dala  Kalk  AB 

60.62000 

15.10000 

open  pit 

202 

Dannemora 

60.12000 

17.52000 

203 

Garpenberg  Nord 

60.20000 

16.13000 

204 

Garpenberg 

60.18500 

16.11500 

205 

Zinkgruvan 

58.80000 

15.10000 

206 

Glanshammar 

59.30000- 

15.40000 

207 

Grangesberg 

60.10000 

15.00000 

208 

Yxsjoberg 

60.00000 

14.80000 

209 

Gasgruvan 

59.70000 

14.10000 

210 

Svardsjo 

60.70000 

15.90000 

Table  7.5.1  Example  of  attributes  contained  in  the  MINFO  relation.  Mines 
1-100  are  Norwegian,  mines  101-200  Finnish,  and  mines  201-300  are  Swedish. 


Epoch 

Mine  No. 

Charges  Charge  No. 

Delay, ms 

Configur. 

610812300.000 

1 

9 

1 

0.000 

ripple 

610812300.000 

1 

9 

2 

45.000 

ripple 

610812300.000 

1 

9 

3 

90.000 

ripple 

610812300.000 

1 

9 

4 

135.000 

ripple 

610812300.000 

l 

9 

5 

180.000 

ripple 

610812300.000 

1 

9 

6 

225.000 

ripple 

610812300.000 

1 

9 

7 

270.000 

ripple 

610812300.000 

1 

9 

8 

315.000 

ripple 

610812300.000 

1 

9 

9 

360.000 

ripple 

612021900.000 

1 

8 

1 

0.000 

ripple 

612021900.000 

1 

8 

2 

45.000 

ripple 

612021900.000 

1 

8 

3 

90.000 

ripple 

612021900.000 

1 

8 

4 

135.000 

ripple 

612021900.000 

1 

8 

5 

180.000 

ripple 

612021900.000 

1 

8 

6 

225.000 

ripple 

612021900.000 

1 

8 

7 

270.000 

ripple 

612021900.000 

1 

8 

8 

315.000 

ripple 

507291300.000 

3 

22 

0 

0.000 

obtainable 

508768800.000 

3 

24 

0 

0.000 

obtainable 

509991000.000 

3 

21 

0 

0.000 

obtai  liable 

516042900.000 

3 

26 

0 

0.000 

obtainable 

518026200.000 

3 

27 

0 

0.000 

obtainable 

526588200.000 

3 

21 

0 

0.000 

obtainable 

578831400.000 

3 

14 

0 

0.000 

obtainable 

581738400.000 

3 

20 

0 

0.000 

obtai  liable 

584957700.000 

3 

34 

0 

0.000 

obtainable 

592677600.000 

3 

10 

0 

0.000 

obtainable 

595688100.000 

3 

16 

0 

0.000 

obtainable 

604710000.000 

3 

17 

0 

0.000 

obtainable 

609668100.000 

3 

26 

0 

0.000 

obtainable 

.532384500.000 

4 
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0 

0.000 

obtainable 

613826100.000 

5 

0 

0 

0.000 

obtainable 

602370000.000 

6 

17 

1 

0.000 

ripple 

602370000.000 

6 

17 

2 

100.000 

ripple 

602370000.000 

6 

17 

3 

200.000 

ripple 

602370000.000 

6 

17 

4 

300.000 

ripple 

602370000.000 

6 

17 

5 

400.000 

ripple 

602370000.000 

6 

17 

6 

450.000 

ripple 

602370000.000 

6 

17 

7 

950.000 

ripple 

602370000.000 

6 

17 

8 

1450.000 

ripple 
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Epoch 

Mine  No. 

Charges 

Charge  No. 

Delay, ms 

Configi 

602370000.000 

6 

17 

9 

1950.000 

ripple 

602370000.000 

6 

17 

10 

2450.000 

ripple 

602370000.000 

6 

17 

11 

2950.000 

ripple 

602370000.000 

6 

17 

12 

3450.000 

ripple 

602370000.000 

6 

17 

13 

3950.000 

ripple 

602370000.000 

6 

17  - 

14 

4450.000 

ripple 

602370000.000 

6 

17 

15 

4950.000 

ripple 

602370000.000 

6 

17 

16 

5450.000 

ripple 

602370000.000 

6 

17 

17 

5950.000 

ripple 

599612400.000 

9 

1 

1 

0.000 

single 

Table  7.5.2  Example  showing  some  of  the  attributes  of  the  MINEX  relation. 
The  middle  section  of  the  table  contains  data  of  ripple  fired  charges  where  the 
details  of  ripple  firing  is  not  specified,  however,  the  data  may  be  obtainable 
from  the  mining  company. 
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Fig.  7.5.1  Mines  and  quarries  in  Scandinavia  presently  registered  in  the  min¬ 
ing  explosion  data  base.  The  most  important  mines  in  Norway  are  shown  as 
Titania.  Olivin,  Rana  Gruber  and  Sydvaranger. 
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Fig.  7.5.4  Distance  from  ARCESS/N’ORESS  to  four  of  the  main  mines  in 
Norway.  Location  of  the  mines  is  shown  in  Fig.  7.5.1 
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Fig.  7.5.5  NORESS  single  channel  recording  of  ripple  fired  explosion  at  Sand- 
vika  (near  Oslo),  shown  in  nine  filter  bands.  Charge  was  1.615  tons  over  68 
delays  in  the  time  interval  0-578  ms. 
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Fig.  7.5.7  ARCESS  recording  of  ripple  fired  explosion  at  Rana  Gruber. 
Charge  was  180  tons  over  34  delays  in  the  time  interval  0-1530  ms. 
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APPENDIX  A 


Note  on  Data  Base  Structure 
Proposal  for  Mine  Explosion  Relations 


The  present  note  proposes  two  relations  in  conformance  with  CSS  Database 
Structure,  Version  2.8.  The  relations  and  attributes  are  those  proposed  by 
NORSAR  based  on  a  sparse  return  of  questionnaires  submitted  to  67  mines 
and  quarries.  The  answering  per  cent  is  around  25.  and  it  may  prove  necessary 
to  include  more  attributes  during  future  research  and  contact  with  the  mining 
companies. 

The  relations  proposed  are  denoted  minex  and  minfo.  Minfi}  contains 
the  basic  information  about  each  mine  not  pertaining  to  specific  explosions. 
The  relation  with  minex  is  secured  via  the  relational  key  minid.  the  unique 
identification  number  of  each  mine.  Further  attributes  contained  in  minfo  are 
defined  below. 

The  minex  relation  contains  the  relational  evid  key  and  other  attributes 
proposed  for  mine  explosions  (see  below).  The  minex  relation  describes  each 
detonated  charge  in  ripple  firing  by  one  tuple.  The  number  of  tuples  pertaining 
to  a  ripple  fired  explosion  thus  equals  the  number  of  charges  in  the  ripple. 
The  charge  configuration  attribute  chacon  characterizes  the  charge  as  being 
'single’,  Tipple’,  'unknown-  or  'obtainable1.  In  the  case  of  'unknown',  and 
rdur  given,  the  charge  is  a  ripple  but  the  details  of  the  ripple  configuration 
is  only  known  as  the  time  from  first  to  last  subcharge  in  the  ripple.  In  the 
case  of  'obtainable’,  the  ripple  description  may  be  requested  from  the  mining 
company. 
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minex 


attribute 

field 

storage 

external 

name 

no. 

type 

format 

evid 

1 

i4 

i8 

time 

2 

f  8 

f  15.3 

date 

3 

i4 

i8 

minid 

4 

i4 

i8 

depth 

5 

f  4 

f  9 .4 

elev 

6 

f  4 

f  9 .4 

lat 

7 

f  4 

f  10.5 

Ion 

8 

f  4 

f  10.5 

extyp 

9 

cl5 

al5 

grade 

10 

f4 

f9 .4 

nex 

11 

i4 

i8 

subnex 

CNJ 

*-< 

i4 

i8 

delt 

13 

f  4 

f9.3 

rdur 

14 

f  4 

f  9 . 3 

totcha 

15 

f  4 

f  9 . 1 

charge 

16 

f  4 

f  9 . 1 

chacon 

17 

clO 

alO 

auth 

18 

c  15 

al5 

character 

attribute 

positions 

description 

1-8 

event  id 

10-24 

epoch  time  of  explosion 

26-33 

shot  date  (julian) 

35-42 

mine  id 

44-52 

shot  depth  (km) 

54-62 

surface  elevation  (km) 

64-73 

latiw^de  (geodetic) 

75-84 

longitude  (geodetic) 

86-100 

type  of  explosive 

102-110 

strength  relative  tnt 

112-119 

tot  no.  charges  in 
ripple 

121-128 

actual  charge  number 

130-138 

delay  rel.  first 
charge  (ms) 

140-148 

ripple  duration  (ms) 

150-158 

total  charge  (kg) 

160-168 

indiv.  ripple 
charge  (kg) 

170-179 

charge  config. 
description 

181-195 

author 
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minfo 


attribute 

field 

storage 

external 

^character 

attribute 

name 

no . 

type 

format 

positions 

descri-ption 

minid 

1 

i4 

i8 

1-8 

mine  id 

minam 

2 

cl5 

alS 

10-24 

name  of  mine 

lat 

3 

f  4 

fl0.5 

26-35 

latitude  (geodetic) 

Ion 

4 

f  4 

f  10.5 

37-46 

longitude (geodetic) 

elev 

5 

f  4 

f  9 . 5 

48-56 

surface  elev.  (km) 

prodpt 

6 

f  4 

f  9 . 5 

58-66 

production  depth  (km) 

mintyp 

7 

cl5 

al5 

68-82 

mine  type 

prodct 

8 

c20 

a20 

84-103 

product 

geolog 

9 

c30 

a30 

105-134 

bedrock  geology 

f irprc 

10 

c40 

a40 

136-145 

firing  practice 

auth 

11 

CIS 

al5 

177-191 

author 
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Glossary  of  attributes 


minid  Relations  :  miner  minfo 

Mine  id;  each  mine  is  assigned  a  positive  unique  number  which 
identifies  it  in  the  data  base.  (i4)  No  null. 


extyp  Relation  ;  minex 

Explosive  type;  type  of  explosive  used  in  a  charge. (cl?)  Mixed 
case . 


grade  Relation  :  minex 

Explosive  strength;  messured  relative  to  TNT=1.  Slurry  may  be  of 
grade  0.9,  the  unit  weight  of  slurry  is  equivalent  to  0.4  units  of 
TNT.  used  to  normalize  charge  strength  to  common  reference 
strength.  (f4)  Null  =  -1. 


nex  Relation  :  minex 

Total  number  of  charges  in  a  ripple  firing.  (i4)  Null  =  -1. 


subnex  Relation  :  minex 

Charge  number  in  ripple  firing;  sequencial  unique  number  identify¬ 
ing  the  acctual  charge  for  the  current  entry,  values  1  -  nex. 
subnex=0  means  obtainable  charge  distribution  in  a  ripple  firing. 
(i4)  Null  =  -1. 

delt  Relation  :  minex 

Delay  time  in  ripple  firing;  delay  for  current  charge  (entry) 
relative  to  first  charge  measured  in  miliseconds.  (f4)  Null  =  -1. 


rdur  Relation  :  minex 

Duration  of  ripple  firing  when  chacon  is  ’unqnown’ ; 

relative  to  first  charge  measured  in  miliseconds.  (f4)  Null  =  -1. 
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totcha  Relation  :  minex 

Total  charge;  sum  of  nex  individual  charges  (charge)  measured  in 
kg.  For  single  charges  totcha=charge .  (f4)  Null  =  -1. 


charge  Relation  :  minex  - 

Individual  charge  size  in  a  ripple  firing  measured  in  kg.  (f4) 
Null  =  -1. 


cuacon  Relation  :  minex 

Charge  configuration  attaining  values  ‘single’,  ‘ripple’, 
‘unknown’ ,  ‘obtainable’  (clO)  lower  case 

minam  Relation  :  minfo 

Name  of  mine;  (clS)  Mixed  case 

prodpt  Relation  :  minfo 

Production  depth  below  surface  in  km.  (f4)  Null  =  -1. 


mintyp  Relation  :  minfo 

Type  of  mine  ;  Pit  or  shaft  (cl5)  Case  =  lower 

prodct  Relation  :  minfo 

Type  of  product  of  the  mine;  Ore  (type),  rock  etc.  (c20) 
geolog  Relation  :  minfo 

Type  of  geological  environment;  sedimentary,  crystalline, 
metamorphic,  specific  rock  types  etc.  (c30) 


firprc  Relation  :  minfo 

Firing  practice  information;  Time  of  day,  regularity,  etc.  (c40) 


102 


7.6  Transverse  components  of  explosion-induced  Lg  waves 
and  upper  crustal  anisotropy 

Introduction 

Vertical  components  of  Lg  waves  have  been  successfully  modelled  by  many 
authors  (e.g.,  Bouchon,  1982;  Herrmanir&  Kijko,  1983),  but  specific  problems 
occur  when  one  attempts  to  model  also  their  horizontal  components.  Syn¬ 
thesis  of  Lg  waves  produced  by  explosions  and  propagating  in  an  isotropic 
vertically  layered  crust  predicts  displacement  on  the  vertical  and  radial  com¬ 
ponents  only.  Lg  waves  produced  by  explosions,  however,  commonly  show 
transverse  components  as  large  or  even  larger  than  the  vertical  and  radial 
ones  both  in  Eurasia  (Mykkeltveit  &  Husebye,  1981)  and  in  the  United  States 
(Blandford,  1981).  Examples  of  Lg  recordings  at  NORSAR  3-component  short 
period  station  01A  are  given  in  Fig.  7.6.1  for  3  explosions  fired  at  sea  during 
the  CANOBE  and  FENNOLORA  refraction  experiments  (Cassell  et  al,  1983; 
Ansorge,  1981).  The  data  have  been  band-pass  filtered  between  0.5  and  2 
Hz.  Note  in  particular  that  the  Lg  waves  arrive  on  the  transverse  components 
as  early  as  on  the  vertical  or  radial  ones  and  cannot  therefore  be  considered 
only  as  coda  waves.  These  recordings  are  representative  of  a  larger  collection 
of  NORSAR  explosion  recordings  where  Lg  wave  transverse  components  are 
systematically  present  and  large.  In  the  present  case,  the  explosions  were  fired 
close  to  the  sea-bed,  73  km  from  the  coast  for  the  CANOBE  example,  and 
close  to  the  coast  for  the  FENNOLORA  recordings.  The  near-source  structure 
certainly  generates  SV  waves,  but  less  likely  large  SH  waves.  Propagation  in 
a  complex  structure  is  required  to  produce  these  transverse  components. 

The  heterogeneous  nature  of  the  crust  is  more  commonly  invoked  than 
anisotropy  in  order  to  explain  complex  polarization  of  short  period  waves, 
though  anisotropy  has  been  observed  in  various  regions  of  the  world,  and  its 
effect  on  polarization  of  body  waves  and  surface  waves  is  well-documented 
(see  Crampin,  1977,  for  a  review).  We  plan  here  to  contribute  to  a  better 
understanding  of  Lg  wave  propagation  by  presenting  synthetic  seismograms  of 
explosion-induced  3-component  Lg  waves  propagating  in  models  with  apparent 
upper  crustal  anisotropy. 

Anisotropy  in  the  upper  crust 

Shear-wave  splitting,  the  most  typical  signature  of  anisotropy,  is  observed 
in  an  increasing  number  of  regional  and  local  studies.  Of  special  interest  to  us 
is  the  analysis  of  Brooks  et  al  ( 1987)  who  show  evidence  of  crustal  anisotropy 
in  Scandinavia.  Anderson  et  al  (1974)  and  Hudson  (1981)  showed  how  a 
uniform  distribution  of  aligned  cracks  in  an  otherwise  isotropic  matrix  can  be 
viewed  in  the  long  wavelength  approximation  as  a  homogeneous  anisotropic 
material.  Crampin  et  al  (19841  proposed  that  widespread  crustal  anisotropy 
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coherent  over  large  areas  proceeds  from  opening  of  microcracks  in  preferred 
directions  related  to  the  regional  stress  field.  The  polarization  directions  of 
observed  split  S-waves  always  support  this  hypothesis,  which  has  become  the 
most  commonly  accepted  one  to  explain  upper  crustal  anisotropy. 

The  nature  and  depth  of  the  anisotropy  are  very  seldom  and  poorly  con¬ 
strained  by  S-wave  splitting  data.  AnisoUopy  confined  to  the  first  10  km 
of  the  crust  is  the  best  model  for  Japan  after  Kaneshima  et  al  (1988),  and 
anisotropy  to  at  least  6  km  depth  has  been  found  in  the  crystalline  base¬ 
ment  of  the  Urals  by  Koshubin  et  al  (1984).  When  the  depth  is  not  directly 
constrained,  an  anisotropy  in  the  first  few  kilometers  of  the  crust  is  usually 
consistent  with  the  observations.  Observed  shear  wave  delays  are  often  around 
0.2  s  (Kaneshima  et  al,  1988:  Zollo  k  Bernard.  1989:  Brooks  et  al,  1987)  and 
suggest  an  open  crack  density  of  0.1  over  5  km  depth  (or  any  combination 
giving  the  same  delay  such  as  0.05  crack  density  over  10  km).  More  directly 
constrained  velocities  support  a  crack  density  averaging  0.05  (Babuska  k  Pros. 
1984;  Koshubin  et  al.  1984). 

Model  and  modelling  method 

In  accordance  with  the  foregoing,  we  assume  that  the  upper  crust  in  Scan¬ 
dinavia  has  an  apparent  anisotropy  due  to  a  uniform  distribution  of  circular 
vertical  cracks  with  small  aspect  ratio  and  preferred  orientation  related  to 
the  regional  stress  field,  as  indicated  schematically  in  Fig.  7.6.1.  The  formula 
for  calculating  the  effective  elastic  coefficients  of  a  medium  containing  dry  og 
liquid-filled  cracks  are  given  in  Crampin  (1984)  after  Hudson  (1981).  We  re¬ 
tain  only  the  purely  elastic  perturbation  due  to  the  cracks,  and  neglect  the 
anelastic  part  since  one  can  estimate  from  Crampin  (1984)  that  metric  cracks 
with  a  0.1  crack  density  lead  to  an  attenuation  factor  Q~l  of  the  order  10-n 
for  the  1  Hz  Lg  wave.  The  attenuation  depends  on  the  crack  radius  in  the 
power  of  3  and  is  therefore  even  smaller  for  microcracks. 

Our  basic  model  for  the  Scandinavian  crust  is  taken  from  Cassell  et  al 
( 1983)  who  interpreted  P-wave  arrival  times  and  amplitudes  from  the  CANOBE 
refraction  experiment  along  a  profile  which  exactly  fits  the  propagation  path 
of  the  Lg  wave  from  CANOBE  shot  H2  to  NORSAR  (Fig.  7.6.1).  We  have 
replaced  the  upper  5  km  of  their  model,  unconstrained  by  their  data,  by  a 
layer  having  P-wave  velocities  increasing  from  6  to  6.25  km/s.  This  is  similar 
to  what  was  found  by  Gundem  (1984)  in  southern  Norway,  and  accounts  bet¬ 
ter  for  the  group  velocity  and  short  duration  of  the  observed  Lg  wavetrains 
in  this  area.  We  derive  the  S-wave  velocity  model  from  the  P-wave  velocity 
model  assuming  a  Poisson  ratio  equal  to  0.25.  Quality  factors  are  introduced 
in  the  model  with  frequency  dependence  of  \J~j  after  Campillo  et  al  (1985). 
The  model  is  displayed  in  Fig.  7.6.2. 
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The  Lg  modes  in  the  frequency  band  0.5  to  2  Hz  are  calculated  in  the 
anisotropic  models  at  first  order  in  deviation  from  isotropy  using  the  quasi- 
degenerate  perturbation  method  (Luh,  1973;  Maupin,  1989).  This  method  is 
capable  of  accounting  for  the  full  3-D  polarization  of  the  Lg  wave  generalized 
modes,  though  retaining  the  advantages  of  treating  anisotropy  as  a  small  per¬ 
turbation  of  an  isotropic  model.  The  source  for  the  explosions  is  taken  as  a 
unit  pressure  dirac  at  a  depth  of  150  m,  chosen  to  mimic  the  refraction  shots 
of  the  CANOBE  profile,  though  very  imperfectly  since  we  cannot  account  for 
the  water  layer  in  which  the  explosions  were  fired  or  the  sedimentary  layer 
underneath.  No  absolute  amplitude  is  tied  to  the  source,  and  amplitudes  of 
the  synthetic  seismograms  are  reported  relative  to  the  maximum  amplitude  of 
the  synthetic  seismogram  in  the  isotropic  reference  structure.  The  ground  dis¬ 
placement  is  convolved  with  the  NORSAR  short-period  instrument  response. 

Explosion  Lg  waves  and  upper  crustal  anisotropy 

Synthetic  Lg  waves  calculated  in  models  having  various  crack  densities 
at  various  depths  are  presented  in  Fig.  7.6.3.  These  cases  all  correspond  to  a 
wave  propagation  direction  at  15°  from  the  crack  normal  preferred  orientation, 
similar  to  the  angle  which  can  be  inferred  from  Brooks  et  al  (1987)  between 
crack  orientation  in  Scandinavia  and  propagation  path  to  NORSAR  of  the  Lg 
wave  from  the  CANOBE  shot  H2  and  FENNOLORA  shot  E2  (Fig.  7.6.1).  The 
source-station  distance  is  450  km.  The  reference  synthetic  Lg  wave  calculated 
in  the  isotropic  structure  has  of  course  no  transverse  component  (Fig.  7.6.3a). 
For  liquid-filled  cracks  between  0  and  10  km  depth  and  a  crack  density  of 
0.02,  a  small  transverse  component  appears  (Fig.  7.6.3b).  Increasing  the  crack 
density  to  0.05  (Fig.  7.6.3c),  the  transverse  component  grows  larger  than  the 
radial  one.  Note  that  this  anisotropy,  the  amplitude  and  depth  extension  of 
which  are  compatible  with  observations  of  split  shear  wave  delays  (Brooks  et 
al,  1987),  yields  relative  amplitudes  of  the  three  Lg  components  that  closely 
resemble  the  observations  (Fig.  7.6.1).  In  Fig.  7.6.3d  is  shown  the  synthetic 
Lg  wave  when  the  cracks  are  confined  to  the  upper  5  km  of  the  crust  with  0.1 
crack  density.  The  Lg  wavetrain  is  lengthened  on  the  three  components.  This 
lengthening  appears  also  for  other  propagation  directions.  The  last  example 
(Fig.  7.6.3e)  is  for  a  uniform  distribution  of  dry  cracks  in  the  upper  10  km  of 
the  crust  with  0.05  crack  density.  The  transverse  components  are  not  strongly 
developed  in  that  case. 

Synthetic  Lg  waves  for  wave  propagation  at  other  angles  from  the  crack- 
normal  direction  have  also  been  calculated.  For  propagation  in  the  crack 
orientation  symmetry  planes  at  0  and  90°,  no  transverse  component  of  course 
appears.  For  angles  ranging  from  0  to  45°,  the  Lg  waves  are  like  those  for  a  15° 
angle  and  displayed  in  Fig.  7.6.3.  For  angles  approaching  90°.  the  amplitude 
of  the  transverse  component  decreases  gradually  and  its  maximum  is  shifted 
towards  later  arrival  times.  However,  the  onset  time  of  the  Lg  wave  or  group 
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velocity  of  its  first  wave  packet  does  not  depend  on  azimuth.  The  effect  of 
upper  crustal  anisotropy  on  the  Lg  wave  is  basically  a  polarization  effect  and 
not  a  velocity  effect. 

Assuming  a  uniform  anisotropy  along  the  whole  propagation  path,  we  have 
calculated  Lg  synthetic  seismograms  for  different  epicentral  distances.  Be¬ 
sides  fluctuations  due  to  variations  in  modal  interferences,  which  occur  also 
in  isotropic  structures,  it  appears  that  the  ratio  between  different  components 
is  stable  with  epicentral  distance.  For  an  epicentral  distance  of  150  km,  at 
which  the  Lg  wavetrain  becomes  an  “individualized”  wavetrain,  the  ratio  be¬ 
tween  the  transverse  and  vertical  components  in  our  synthetic  data  is  already 
the  ratio  which  is  measured  at  750  km.  whether  large  or  small.  This  is  in 
agreement  with  observations  made  bv  Blandford  (1981). 

The  extensive  dilatancy  anisotropy  hypothesis  of  Crampin  et  al  (.984)  on 
which  we  have  based  our  models,  and  the  widespread  character  of  transverse 
component  observations,  rather  suggest  a  regional  crustal  anisotropy  than  a  lo¬ 
cal  one.  We  can,  however,  test  the  minimum  lateral  extension  of  an  anisotropic 
zone  necessary  to  explain  the  data.  As  for  body  waves,  if  anisotropy  along 
only  part  of  the  wavepath  introduces  polarization  anomalies,  they  are  propa¬ 
gated  to  the  station  even  if  the  remainder  of  the  propagation  path  is  isotropic. 
The  analysis  of  Lg  wave  character  as  a  function  of  epicentral  distance  implies 
that  anisotropy  over  150  km  is  sufficient  to  produce  the  observed  polarization 
anomalies. 

Conclusion 

We  have  shown  that  an  upper  crust  with  a  uniform  distribution  of  microc¬ 
racks  aligned  in  the  regional  stress  field  can  account  for  transverse  components 
of  Lg  waves  induced  by  explosions.  Our  preferred  model  to  explain  with  upper 
crustal  anisotropy  the  observations  in  Scandinavia  is  a  distribution  of  liquid- 
filled  cracks  in  the  upper  10  km  of  the  crust  with  a  crack  density  larger  than 
0.02  and  most  favorably  around  0.05.  This  model  is  consistent  with  observa¬ 
tions  of  shear  wave  splitting  by  Brooks  et  al  (1987)  in  northern  Scandinavia. 


Valerie  Maupin,  Inst,  de  Physique  du  Globe.  Strasbourg 
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Fig.  7.6.1  3-component  Lg  recordings  at  NORSAR  short-period  station  01A 
of  CANOBE  shot  H2  and  FENNOLORA  shots  B3  and  E2,  band-pass  filtered 
between  0.5  and  2  Hz.  The  arrival  time  of  the  wave  with  group  velocity 
3.6  km/s  is  indicated.  The  amplitude  scale  is  identical  on  the  three  components 
for  each  recording.  The  direction  of  regional  maximum  horizontal  compression 
and  a  schematic  vertical  crack  opened  by  the  regional  stress  field  are  indicated 
in  the  upper  left  corner  of  the  map  of  southern  Scandinavia. 
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Fig.  7.6.3  3-component  synthetic  Lg  waves  produced  bv  an  explosion  at  150  m 
depth  and  450  km  epicentral  distance  in  different  models  of  southern  Scandi¬ 
navia.  (a)  Isotropic  model;  (b)  model  with  liquid-filled  circular  vertical  cracks 
between  0  and  10  km  depth,  0.02  crack  density  and  orientation  of  the  crack 
normal  at  15°  from  the  wave  propagation  direction;  (c)  the  same  as  for  case 
(b)  for  0.05  crack  density;  (d)  the  same  as  for  case  (b)  for  0.1  crack  density  and 
cracks  between  0  and  5  km  depth;  (e)  the  same  as  case  (c)  for  dry  cracks.  The 
arrow  on  the  vertical  traces  indicates  the  arrival  time  of  the  3.6  km/s  group 
velocity  wave.  The  three  components  of  each  plot  have  the  same  amplitude 
scale.  The  ratio  between  the  maximum  amplitude  of  the  largest  component 
and  the  maximum  amplitude  of  the  vertical  component  in  the  isotropic  case 
is  indicated  in  each  3-component  plot  at  the  end  of  the  vertical  trace. 
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7.7  Statistically  optimal  event  detection  using  small  array 
data 

A  generalization  of  Capon's  maximum  likelihood  technique  for  detection  and 
estimation  of  seismic  signals  is  introduced.  By  using  a  multidimensional  au¬ 
toregressive  approximation  of  seismic  array  noise,  we  have  developed  a  tech¬ 
nique  to  use  Capon’s  group  filter  for  on-line  processing.  Such  autoregressive 
adaptation  to  the  current  noise  matrix  power  spectrum  yields  good  suppres¬ 
sion  of  mutually  correlated  array  noise  processes.  An  example  is  shown  of  this 
technique  as  applied  to  detection  of  the  a  small  Semipalatinsk  underground 
explosion  recorded  at  the  ARCESS  array. 

Nuclear  explosion  monitoring  using  seismic  data  is  faced  with  the  problem 
that  signals  of  small  explosions  are  masked  by  noise,  and  thus  have  to  be 
extracted  using  features  of  both  the  noise  and  the  signal.  Small  arrays  appear 
to  be  especially  suited  for  that  purpose.  This  is  due  to  the  strong  correlation  of 
noise  between  different  closely  located  receivers  which  gives  us  an  opportunity 
to  obtain  significant  noise  suppression.  To  realize  this  opportunity,  special 
software  is  needed. 

We  have  to  solve  two  main  problems:  (1)  to  detect  event  signals,  and  (2) 
to  classify  detected  signals  as  originating  from  either  an  explosion  or  an  earth¬ 
quake.  The  second  task  is  very  complex.  The  proper  identification  needs  first 
of  all  estimation  of  signal  waveform  and  then  estimation  of  signal  parameters 
such  as  onset  times  of  different  phases,  power  and  spectral  features,  and  so  on. 
While  the  detection  of  a  signal  can  be  done  in  a  relatively  narrow  frequency 
band  (for  example  using  a  high  frequency  band  only)  the  classification  must 
in  principle  be  based  on  wide  band  methods.  This  is  so  because  bandpass  fil¬ 
tering  distorts  not  only  the  noise  but  also  the  signal,  thus  possibly  eliminating 
useful  classification  features. 

The  problems  mentioned  can  be  formulated  in  terms  of  mathematical 
statistics,  and  optimal  decision  rules  for  these  mathematical  problems  have 
to  be  found.  The  first  task  is  to  detect  the  signed.  The  time  series  received  by 
the  array  has  the  following  structure:  A(  =  St  +  £(,  t  6  ...  -  1,0,1,...,  where 
5(  is  signal  and  is  noise.  At  is  a  vector  of  different  receiver  outputs.  These 
outputs  are  observed  through  a  moving  window.  Using  data  in  the  window 
we  must  make  a  decision:  does  it  contain  a  signal  or  not.  This  problem  can 
be  solved  in  terms  of  statistical  hypothesis  testing  theory.  It  is  necessary  to 
test  hypothesis  Ho'  observations  in  the  moving  window  are  pure  noise,  versus 
hypothesis  H\\  they  comprise  signal  plus  noise.  We  consider  hypothesis  Ho 
to  be  simple,  but  hypothesis  H\  to  be  complex.  This  is  because  the  statistical 
characteristics  of  noise  can  be  measured  before  the  signal  arrival  in  a  stage  of 
adaptation,  but  the  statistical  features  of  seismic  signals  are  almost  completely 
unknown. 
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The  problem  is  now  to  choose  a  decision  function  which  would  provide  the 
smallest  average  error  probability  for  all  possible  signals.  This  problem  can 
be  solved  using  a  Bayesian  approach  and  linite  dimensional  parameterization 
of  the  signal  (Kushnir  and  Lapshin,  1981)  (see  Appendix).  The  decision  has 
the  form  of  an  algorithm  as  shown  in  Fig.  7.7.1.  It  consists  of  a  group  filter 
followed  by  autocorrelators,  calculation  of  a  quadratic  form  and  a  trigger  which 
compares  the  quadratic  form  with  a  threshold.  The  transfer  function  of  the 
group  filter  is  a  vector  described  as  a  product  of  the  inverse  matrix  power 
spectrum  F~{  of  the  noise  and  a  vector  G  of  phase  shift  factors  due  to 
signal  delays  The  algorithm  shown  has  a  form  which  can  be  easily  realized 
as  an  on-line  procedure.  This  is  mostly  due  to  the  use  of  multidimensional 
autoregressive  (AR)  estimation  of  the  matrix  power  spectrum  density  F  of 
the  noise.  This  allows  us  to  avoid  direct  inversion  of  spectral  power  matrices 
and  is  very  convenient  in  an  adaptive  procedure  (Havkin.  1979;  Kushnir  et  at. 
1980). 

The  second  problem  mentioned  above  is  to  estimate  properly  the  signal 
waveform.  The  main  purpose  of  this  estimation  is  to  make  a  correct  decision 
of  whether  we  have  an  explosion  or  an  earthquake.  We  must  do  it  as  precisely 
as  possible.  The  model  of  observation  here  has  the  following  form:  Xt  =  Gt  * 
Ut  +  where  *  is  a  sign  of  convolution.  ut  is  a  signal  waveform  to  be  estimated 
—  (the  particle  motion  along  the  seismic  ray),  Gt  is  the  transfer  function  on 
the  path  from  the  seismic  source  to  the  receiver.  For  plane  waves,  this  vector  is 
defined  by  time  delays  only.  For  solving  the  estimation  problem,  a  conditional 
Wiener  filter  can  be  constructed.  This  filter  minimizes  the  variance  of  the 
estimate  E{(ut  —  ut  )2}  under  the  condition  that  the  mathematical  expectation 
of  the  estimate  coincides  with  the  real  signal:  E{ut}  =  ut. 

It  so  happens  that  the  Wiener  filter  we  are  looking  for  consists  of  the  same 
group  filter  which  is  used  for  the  detection  followed  bv  a  restitution  filter 
(which  makes  it  possible  to  obtain  the  signal  undistorted).  Capon  was  the 
first  to  propose  this  filter  for  seismic  signal  extraction  from  array  data  (Capon. 
1970).  The  complete  array  procedure  for  detection  and  classification  is  shown 
in  Fig. 7. 7. 2.  We  adapt  to  the  noise  matrix  power  spectrum  by  estimating  its 
AR  parameters  and  computing  vector  coefficients  of  the  group  filter.  Then 
we  perform  group  filtration  in  the  moving  window  and  detect  the  signal.  The 
first  two  operations  are  made  periodically  according  to  the  interval  of  noise 
stationarity.  The  third  —  group  filtration  —  can  be  devised  as  an  on-line 
procedure.  After  the  signal  is  detected  it  must  be  filtered  by  the  restitution 
filter,  which  refines  its  shape.  And  finally  classification  may  be  done  using  the 
estimates  of  signal  parameters. 

The  software  designed  was  tested  using  simulated  data  with  the  aim  of 
comparing  its  actual  efficiency  with  the  theoretical  one.  The  results  of  these 
tests  are  shown  in  Fig.  7. 7. .'5.  The  depicted  curves  are  the  gains  in  power  signal- 
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to-noise  ratio  of  undisturbed  optimal  group  filtration  relative  to  a  conventional 
beam  versus  the  coherence  coefficient  of  the  noise.  The  latter  is  defined  as 
the  ratio  between  the  largest  and  the  smallest  eigenvalues  of  the  noise  power 
spectrum  matrix.  We  see  that  the  mentioned  gain  may  be  very  significant  if 
the  array  noise  is  coherent  enough.  This  happens  in  practice  at  small  aperture 
arrays.  The  .z'culations  were  made  for  the  central  subarray  of  the  NORSAR 
array. 

Highly  promising  results  were  obtained  by  the  use  of  ARCESS  data  for 
signals  from  one  of  the  smallest  nuclear  tests  known  to  have  been  conducted 
at  the  Semipalatinsk  test  site.  In  Fig.  7.7.4  we  display  the  records  for  four 
ARCESS  channels  and  note  that  the  signal  is  obscured  by  the  noise.  For  the 
conventional  beam  trace  ( Fig.  7.7.5)  the  signal  is  likewise  not  seen,  but  inspec' 
ing  the  output  of  the  undisturbed  group  filter  we  can  see  the  signal  clearly. 
The  power  signal-to-noise  ratio  gain  relative  to  the  beam  is  approximately  a 
nctor  of  70-80  and  it  is  140-160  when  compared  with  a  single  channel.  The 
trace  shown  is  calculated  usi  g  6  matrix  AR  parameters  of  noise.  We  also 
used  other  numbers  of  AR  parameters  and  the  results  seemed  to  be  stable. 
Such  high  suppression  of  noise  is  achieved  mainly  due  to  the  high  correlation 
of  noise  records  in  the  inner  ARCESS  stations  (see  Fig.  7.7.4.  traces  AO. A3). 

Fig.  7.7.6  shows  that  the  6roup  filter  used  does  in  fact  retain  the  shape  of 
the  original  wave  rni.  The  first  trace  is  a  wide  frequency  band  waveform  used 
for  simulating  the  plane  wave  arriving  at  ARCESS.  These  simulated  data  were 
processed  by  the  gr^-n  filter  used  for  the  previously  shown  signed  extraction. 
The  resulting  (second)  trace  practically  coincides  with  the  first.  So.  if  the  real 
signal  is  plane  wave,  it  will  be  undisturbed  by  the  group  filter  in  the  frequency 
band  from  0.5  to  5  Hz. 

The  conventional  method  used  for  the  detection  of  weak  signals  is  the  fil¬ 
tration  of  the  array  beam  in  a  band  of  optimum  signal-to-noise  ratio  ( Kvaerna. 
1989).  The  two  traces  at  the  bottom  of  Fig.  7.7.5  show  the  signal  filtered  in 
the  frequency  band  2.5-4  Hz  after  beamforming  and  after  undisturbed  group 
filtration.  The  gain  here  is  not  as  large  as  in  broad  band,  but  still  exceeds  a 
factor  of  5  in  power  SNR.  Fig.  7.7.7  (at  the  top)  shows  the  same  signals,  but 
in  another  time  scale.  The  chosen  frequency  band  seems  to  be  the  best  for 
filtering  the  signal  after  beamforming.  For  comparison  we  have  plotted  two 
traces  at  the  bottom  of  Fig.  7.7.7  presenting  the  same  signals  filtered  in  the 
frequency  band  3-5  Hz. 

For  the  detection  of  signals  in  our  experiments.  4  different  variants  of  the 
optimal  detector  previously  described  were  used.  All  of  these  detectors  are 
sensitive  not  only  to  the  increase  in  trace  power  due  to  signal  arrival,  but  also 
to  changes  of  the  trace  spectrum  ( Kushnir  et  al,  1983).  The  first  is  optimal  in  a 
statistical  sense,  the  second  is  a  modification  of  STA/LTA  using  prewhitening 
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of  the  noise,  and  the  last  two  are  components  of  the  first.  Fig.  7.7.8  shows 
how  the  detectors  work  when  applied  to  the  beam  and  group  filtei  outputs. 
The  gain  due  to  optimal  group  filtration  is  evident. 

The  final  picture  (Fig.  7.7.9)  shows  the  results  of  the  signal  onset  time 
estimation.  Estimation  is  performed  by  an  algorithm  based  on  the  maximum 
likelihood  method  applied  to  the  problem  of  estimating  the  moment  in  time 
when  parameters  of  the  AR  process  are  abruptly  changed  (Pisarenko  et  aL 
1987).  One  can  see  that  the  likelihood  function  maximum  coincides  exactly 
with  the  beginning  of  the  signal. 

Conclusions 


1.  Application  of  an  adaptive  optimal  group  filtration  technique  to  small 
aperture  arrays  can  provide  large  gains  in  SNR  in  comparison  with  con¬ 
ventional  beamforrning  due  to  high  mutual  correlation  of  array  noise. 

2.  By  using  AR  estimation  of  the  power  noise  spectrum  for  group  filter 
adaptation,  we  greatly  reduce  time  and  memory  needed  for  the  adapta¬ 
tion  procedure  while  providing  high  quality  of  noise  suppression. 

Optimal  group  filtration  does  not  distort  the  signal.  Thus,  when  seen 
in  connection  with  1.  and  2.  above,  it  is  clear  that  optimum  filtration 
has  great  advantages  as  a  preprocessor  to  be  applied  prior  to  subsequent 
broad  band  operations  such  as  source  classification. 


In  future  work,  it  is  recommended  that  the  following  studies  be  undertaken: 

1.  To  perform  a  large-scale  experiment  concerning  detection  and  classifi¬ 
cation  capability  of  optimal  group  filtration  applied  to  seismic  signals 
recorded  at  N'ORESS  and  ARCESS.  Comparison  of  false  alarm  rates  of 
the  optimum  detector  versus  conventional  beam  detectors  should  form 
part  of  this  investigation. 

2.  To  develop  algorithms  for  the  compensation  of  signal  frequency  depen¬ 
dent  waveform  distortions  due  to  propagation  in  real  media  under  the 
array.  This  will  give  us  an  opportunity  to  equalize  signal  waveforms  in 
different  array  receivers  and  improve  signal  extraction  at  high  frequen¬ 
cies. 

•'!.  To  implement  these  algorithms  in  an  operational  environment  at  the 
NORSAR  data  center. 


i  it; 


4.  To  develop  algorithms  for  designing  optimal  array  geometry  on  the  basis 
of  the  optimal  group  filtration  features. 

5.  To  use  optimal  group  filtration  for  holographic  investigations  of  the 
earth’s  interior  using  NORSAR  data. 

We  plan  to  continue  further  work  along  these  lines  in  future  co-operative 
projects  between  NORSAR  and  the  Institute  of  Physics  of  the  Earth.  Moscow. 


A.F.  Kushnir,  Inst,  of  Physics  of  the  Earth.  Moscow 
V.I.  Pinsky,  Inst,  of  Physics  of  the  Earth.  Moscow 
J.  Fyen,  NORSAR 
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GROUP  FILTER  F  1  (X)G(X)  DETECTOR 


G (X ) =  (expf-iTjX) ,  1=1, m) 

P 

ha(X)=  2  e  cos(kX) 

O  It  ■  O  k 


Fig.  7.7.1  Optimal  group  detector  flowchart. 
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opCc/na.£  OLyrrciy 


Fig.  7.7.2 
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Fig.  7.7.3 


Fig.  7.7.4 
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UNDISTURBED  GROUP 


Fig.  7.7.5 
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BEAM,  DET.  1, 


BEAM,  DET.  3,  0.5-5  Hz 


BEAM,  DET.  4,  0.5-5  Hz 


UNDIST.  GR.  FILTR.,  Dl,  0.5-5  Hz 


UNDIST.  GR.  FILTR.,  D2,  0.5-5  Hz 


UNDIST.  GR.  FILTR.,  D3,  0.5-5  Hz 


A. 
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UNDIST.  GR.  FILTR.,  D4,  0.5-5  Hz 
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APPENDIX 

Statistically  optimal  event  detection  using  array  data 


Let  us  assume  that  the  noise  =  (£h,  •  ••• ,  £mi  )T  at  the  m  receivers  of  the 
array  is  a  multidimensional  Gaussian  time  series  with  zero  mean  and  matrix 
power  spectral  density  F{  A),  A  G  [0,  2k\  and  that  the  signal  is  represented 
as  a  Gaussian  scalar  process  put(Q)  at  the  seismic  event  source.  The  power 
spectral  density  of  ut(0)  is  denoted  59(A), A  E  [0, 2zr ] ,  where  0  is  an  un¬ 
known  vector  parameter.  The  medium  transfer  functions  are  assumed  to  be 
linear.  Then  the  observations  xt  =  (xu, . .  . .  xmt)T  become  multidimensional 
Gaussian  time  series  and  have  the  form  Gt  *  put  +  and  it  is  easy  to  write 
down  the  likelihood  function  u>(A'/v  )  p,0)  for  the  moving  window  observa¬ 
tions  A'v  =  (xf . iv)T.  Here  G,  =  (G1(, . . . ,  Gm,)T  is  the  vector  impulse 

response  function  of  the  media  along  the  paths  from  the  seismic  source  to  the 
receivers;  *  is  the  sign  of  convolution.  Then  the  hypothesis  H 0  is  that  p  =  0 
and  hypothesis  H\  is  that  p  ^  0.  We  consider  the  function  w(  A',v  |  0)  to  be 
known  through  adapative  esitmation.  If  the  unknown  signal  parameters  have 
an  a  priori  distribution  P(0),  then  the  best  (Bayesian)  test  for  testing  these 
hypotheses  ( providing  the  least  average  signal  miss  probability  for  a  given  false 
alarm  probability)  has  the  form: 


<?(X,v) 


1  (signal  is  present)  if  p(  A',v)  >  ka 
0  (signal  is  absent)  if  p(X/v)  <  ka 


(1) 


where 


MXn  1  M2, 0) 
ur( A'/v  |  0) 


djP(  0 ) 


(2) 


Here,  ka  is  the  detection  threshold,  determined  on  the  basis  of  the  given  false 
alarm  probability  a. 


It  is  practically  impossible  to  devise  an  on-line  algorithm  on  the  basis  of 
the  statistic  ( 1)  when  the  distribution  T(0)  is  arbitrary.  But  in  the  important 
particular  case  of  weak  signal  detection,  where  the  SNR,  p2 / <r2 .  is  sufficiently 
small  and  the  moving  window  size  N  is  sufficiently  large,  the  statistic  p(A',v) 
can  be  simplified.  As  it  is  shown  in  Kushnir  and  Lapshin  (1984),  if  the  like¬ 
lihood  ratio  tu(A’.v  I  7\/7v,0)/ui(AV.O)  in  (2)  is  replaced  by  its  exponential 
approximation,  i.e..  the  statistic  p(X\)  in  (2)  is  replaced  by  the  more  compu¬ 
tationally  convenient  statistic 


r(A'.v)  =  [  exp{7A(A,v.0)-  ^r.v(0)}dP(0)  (3) 

Je>  2 

the  asymptotic  error  probability  limits  for  the  test  (3).  p2  =  -/N.N  —  00, 
stay  constant. 
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The  functions  A(A'\.0)  and  F.v(0)  in  (3)  have  the  following  forms 


i  v 

A( A'.v. 0)  =  — ^  ^^(0)  (  W’x,  |*  -g}(Q)V} 

;=i 


r.v(0) 


4.V 


j  =  i 


where  A(  A'.v,  O)  is  an  asymptotically  sufficient  statistic  of  the  observations 
A'.v,  and 

WJ  =  G*(AJ)/~1(A7),  V'  =  VT;g-(A7).  g}(&)  =  9e(\j) 


r.v(0)  is  the  Fisher  information  quantity,  divided  by  .V. 

G(  A )  is  the  Fourier  transform  of  G<. 

Xj  =  4-  xte.ip{i\,t)  is  the  discrete  Fourier  transform  of  the  obser¬ 

vations  ,Y.\ . 

A,  =  j-: rj. 

F(  A)  is  the  matrix  power  spectral  density  of  the  noise 

For  calculation  of  the  integral  (3)  in  analytic  form  we  will  assume  that: 

P 

5e ( A)  =  1  +  ^  CkCoskX,  9  =  {ci . cP)T  (4) 

k-\ 

where  c^  =  2E {utut+k} / F {u^}  is  the  autocorrelation  of  the  source  signal. 
We  further  assume  that  the  a  prion  distribution  P{9)  in  (3)  is  Gaussian: 

dP{9)  =  {(2xf  dctB}~]/2  -  b)TB~'(9  -  b)}  (5) 


Then  we  have 


log{  r{  x.v) }  =  A'  +  -r(A0-  XT  Ad)+ ]-Sr  AX+ -  XT  AB~' b  =  A'  +  Mx.v)  (6) 


where  A'  is  a  constant  independent  of  A'.v.  A  =  (A] . Ap)1. 


An 


-57?  Ei  "V  o  i* -v. 


.V 


;=i 


Afc  =  — —  yt  \  \v:Vj  i2  -r,)cosa-A,i 
2vTy  fr\  1 


IIS 


6  =  (61,...Jp)t 

6k  =  V?  cos(kXj) 

J  =  l 

A  =  (r  +  7-25-')  =  (BT  Ai~2I)-xB 

j  v  -  _ 

T  =  [—  y  V/  cos(mA;)cos(nAj);  m.n  6  l.P] 
j=i 

Thereby,  the  asymptotically  optimal  Bayes  test  for  array  detection  of  seis¬ 
mic  signals  with  unknown  spectrum  is: 

1  (signal  is  present)  if  ry(x^ )  >  ka 
|  0  (signal  is  absent)  if  r^{xtv)  <  ka 

where  the  threshold  ka  is  determined  on  the  basis  of  the  given  false  alarm 
error  probability,  and  the  fitting  parameter  7  is  determined  so  as  to  provide 
the  highest  asymptotical  efficiency  for  the  test  (Kushnir  et  al,  1983). 

Let  us  consider  in  particular  two  important  (diametrically  opposite)  cases, 
where: 


1.  The  power  spectrum  parameters  of  the  source  signal  are  almost  known: 

II  b  ||«||  r-1  || 

2.  The  a  priori  information  about  these  parameters  is  negligible: 

II  b  ||»  j|  r-1  || 

In  these  cases,  the  statistic  ry(x^)  of  the  test  (7)  is  simplified  and  looks 
like,  respectively: 

1)  rx(xN)  =  A0  +  bT  A 

2 ) r2,(zw)  =  7(^0- ATo)+  ^  ATr-'A  (8) 

Calculation  of  the  test  statistic  (8)  can  be  realized  in  the  time  domain  as 
shown  in  Fig.  7.7.1. 
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